Possibility to insert initial clouds into cold-start #1388

gthompsnWRF · 2021-01-28T03:05:08Z

TYPE: new feature

KEYWORDS: cloud fraction, cloud initialization

SOURCE: Greg Thompson (JCSDA)

DESCRIPTION OF CHANGES:
Problem:
Cold start simulations without microphysics species take significant time to 'spin-up' the clouds causing the radiation
to be way off balance at the start of most model simulations. This greatly affects the first few hours up to 12 hours
or more.

Solution:
The setting of initial cloud fraction and associated cloud water/ice mixing ratio gives a big advantage to 'spin up'
the clouds by assuming that clouds can be created from conditions of very high relative humidity. The cloud fraction
scheme of icloud=3 is a viable solution to making initial clouds along with the setting of the water vapor mixing ratio
to the point of saturation at initialization in association with the inserted cloud fields.

LIST OF MODIFIED FILES:
M Registry/Registry.EM_COMMON
M dyn_em/module_initialize_real.F
M dyn_em/depend.dyn_em
M run/README.namelist

TESTS CONDUCTED:

Numerous simulations in case studies plus real-time simulations with 3km spacing (nearly HRRR-CONUS) in
support of FAA's ICICLE field project (2019Jan28-2019Mar08).
Simulations by P. Jimenez in testing WRF-Solar-V2
Jenkins is pass.

RELEASE NOTE: In combination with updates to the cloud fraction scheme (icloud=3), a new feature permits the usage of the cloud fraction scheme as part of real.exe, specifically designed to use with "cold start" simulations to reduce the spin-up problem of clouds and associated radiation.

…ption

gthompsnWRF · 2021-01-28T03:07:54Z

This PR is a replacement for #1323 which I decided to close because it had the complexity of being based from an earlier code version. I already reported that it had negligible additional run time in real.exe compared to skipping the option. I am just trying to get this approval moving forward with a clean starting point of v4.2.2 as the base in the develop branch. Contact me for any questions. I am trying to dis-entangle multiple PRs into more tiny chunks from here onwards.

davegill · 2021-01-28T04:09:03Z

jenkins:

Please find result of the WRF regression test cases in the attachment. This build is for Commit ID: 19790ec3e234ebb296e6992a8f5f4a19ea5061a4, requested by: gthompsnWRF for PR: https://github.com/wrf-model/WRF/pull/1388. For any query please send e-mail to David Gill.

    Test Type              | Expected  | Received |  Failed
    = = = = = = = = = = = = = = = = = = = = = = = =  = = = =
    Number of Tests        : 19           18
    Number of Builds       : 48           46
    Number of Simulations  : 163           161        0
    Number of Comparisons  : 103           102        0

    Failed Simulations are: 
    None
    Which comparisons are not bit-for-bit: 
    None

New lines in dyn_em/module_initialize_real: ``` USE module_radiation_driver USE module_dm, ONLY : wrf_dm_max_real ``` These require two additional dependencies for module_initialize_real.o to be added in the dyn_em Makefile dependency file. modified: dyn_em/depend.dyn_em

davegill · 2021-01-28T04:25:58Z

@gthompsnWRF
Greg,
Could you fix your affiliation for me please?

davegill · 2021-01-28T04:28:06Z

@gthompsnWRF

I am trying to dis-entangle multiple PRs into more tiny chunks from here onwards.

davegill · 2021-01-28T08:09:05Z

jenkins

Please find result of the WRF regression test cases in the attachment. This build is for Commit ID: 0b491619162d9d4a338712196809a7c65da55f83, requested by: gthompsnWRF for PR: https://github.com/wrf-model/WRF/pull/1388. For any query please send e-mail to David Gill.

    Test Type              | Expected  | Received |  Failed
    = = = = = = = = = = = = = = = = = = = = = = = =  = = = =
    Number of Tests        : 19           18
    Number of Builds       : 48           46
    Number of Simulations  : 163           161        0
    Number of Comparisons  : 103           102        0

    Failed Simulations are: 
    None
    Which comparisons are not bit-for-bit: 
    None

gthompsnWRF · 2021-01-28T13:56:22Z

@gthompsnWRF
Greg,
Could you fix your affiliation for me please?

Got it.

weiwangncar · 2021-01-28T17:38:37Z

I'm ok with this PR.

dyn_em/module_initialize_real.F

davegill · 2021-01-28T20:12:15Z

@gthompsnWRF @dudhia @weiwangncar
Folks,
Does an initialized CLDFRA from real actually get used?

The first time cldfra is used in the first rk step part 1 is in the radiation driver. In the radiation driver, the first time that the CLDFRA field is used is an initialization.

! Choose how to compute cloud fraction (since 3.6)
! Initialize to zero
     DO j=jts,jte
     DO k=kts,kte
     DO i=its,ite
        CLDFRA(i,k,j) = 0.
     END DO
     END DO
     END DO

dudhia · 2021-01-28T20:58:31Z

It only works in conjunction with initial cloud and humidity adjustments that go with it, so that cldfra recomputed in radiation should be the same as the input one in the first time step.

gthompsnWRF · 2021-01-28T21:18:33Z

@gthompsnWRF @dudhia @weiwangncar
Folks,
Does an initialized CLDFRA from real actually get used?

The first time cldfra is used in the first rk step part 1 is in the radiation driver. In the radiation driver, the first time that the CLDFRA field is used is an initialization.
! Choose how to compute cloud fraction (since 3.6)
! Initialize to zero
     DO j=jts,jte
     DO k=kts,kte
     DO i=its,ite
        CLDFRA(i,k,j) = 0.
     END DO
     END DO
     END DO

Dave: The advantage of seeing it done by real.exe is that initial values can be seen in both wrfinput and the very first history file. By waiting until radiation_driver is called, neither happens. The history file is written ahead of the first radiation time step if I recall. I think my own copy of radiation_driver had disabled the zero setting in radiation_driver. I suppose Jimy is right that it just gets made again on first time step - although, this time, if modifying the water vapor, then it will be more cloudy by default.

davegill · 2021-01-28T21:51:46Z

@dudhia @gthompsnWRF

The advantage of seeing it done by real.exe is that initial values can be seen in both wrfinput and the very first history file.

I think I see. This is basically diagnostically providing an initial value for CLDFRA that could be available for plotting purposes. Since this has no real impact on the model solution, I am even less concerned.

gthompsnWRF · 2021-01-28T22:02:39Z

@dudhia @gthompsnWRF

The advantage of seeing it done by real.exe is that initial values can be seen in both wrfinput and the very first history file.

I think I see. This is basically diagnostically providing an initial value for CLDFRA that could be available for plotting purposes. Since this has no real impact on the model solution, I am even less concerned.

Actually, the statement of 'no real impact' is not correct. Since it is being done in real.exe, the fields of cloud water and ice are initialized by cal_cldfra3 as well as a change to the water vapor field (to ensure saturation). So cloud evolution is whoppingly different - that is the point. In cold starts without clouds, it can take 3-15 hours to get clouds forming. With these changes, I get clouds immediately. It is true that some still deteriorate due to downward motions or mixing and dilution. But the change to the first 3 hours is dramatic - giving real-looking clouds at onset that stick around a while as the dynamics of WRF build other (or add onto existing) clouds.

davegill · 2021-01-28T22:45:07Z

@gthompsnWRF @weiwangncar @dudhia

But the change to the first 3 hours is dramatic - giving real-looking clouds at onset that stick around a while as the dynamics of WRF build other (or add onto existing) clouds.

Folks,
This is a namelist option. Do we want it on by default in the Registry? Most everyone runs 5 microphysics species (or more), so we would have an improved cloud state during the early simulation for most users.

dyn_em/module_initialize_real.F

gthompsnWRF · 2021-01-28T22:49:46Z

@gthompsnWRF @weiwangncar @dudhia

But the change to the first 3 hours is dramatic - giving real-looking clouds at onset that stick around a while as the dynamics of WRF build other (or add onto existing) clouds.

Folks,
This is a namelist option. Do we want it on by default in the Registry? Most everyone runs 5 microphysics species (or more), so we would have an improved cloud state during the early simulation for most users.

Dave: I am pretty sure users will not want this enabled by default, since it will be a big new change they may not have expected. In the past, I have seen a lot of hesitancy to enable something new when it makes big changes to users' regular testing simulations.

…ldfra3

weiwangncar · 2021-01-28T22:51:25Z

I agree with Greg. Let's not make it default for now. We can revisit this in the future.

… into insert_init_cloud2

davegill · 2021-02-09T16:57:59Z

@gthompsnWRF @dudhia @weiwangncar
Folks,
Is this part of the modification that is interwoven with a proposal about MAD WRF?
Is this able to be merged?
Where is the input for this option from? Is this is a standard field from a GRIB file?

dudhia · 2021-02-09T17:01:37Z

Yes, this is one we want to resolve after the V4.3 branch since Pedro now has a renamed clfra3 for madwrf routine that we view as temporary.

…

On Tue, Feb 9, 2021 at 9:58 AM Dave Gill ***@***.***> wrote: @gthompsnWRF <https://github.com/gthompsnWRF> @dudhia <https://github.com/dudhia> @weiwangncar <https://github.com/weiwangncar> Folks, Is this part of the modification that is interwoven with a proposal about MAD WRF? Is this able to be merged? Where is the input for this option from? Is this is a standard field from a GRIB file? — You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub <#1388 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AEIZ77FEMJCM62MTWH4RILTS6FSS5ANCNFSM4WWKMFFA> .

gthompsnWRF · 2021-02-09T17:05:15Z

@gthompsnWRF @dudhia @weiwangncar
Folks,
Is this part of the modification that is interwoven with a proposal about MAD WRF?
Is this able to be merged?
Where is the input for this option from? Is this is a standard field from a GRIB file?

Hang on everyone. Let me ask for clarification. My proposal for inserting initial clouds does not "rely" on anything in MADWRF and I want my code to progress asap, earliest branch possible. Did I misinterpret Dave's question or Jimy's reply?

dudhia · 2021-02-09T17:07:29Z

The idea is to commit both yours before we create the V4.3 branch.

…

On Tue, Feb 9, 2021 at 10:05 AM gthompsnWRF ***@***.***> wrote: @gthompsnWRF <https://github.com/gthompsnWRF> @dudhia <https://github.com/dudhia> @weiwangncar <https://github.com/weiwangncar> Folks, Is this part of the modification that is interwoven with a proposal about MAD WRF? Is this able to be merged? Where is the input for this option from? Is this is a standard field from a GRIB file? Hang on everyone. Let me ask for clarification. My proposal for inserting initial clouds does not "rely" on anything in MADWRF and I want my code to progress asap, not in v4.3 branch. Did I misinterpret Dave's question or Jimy's reply? — You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub <#1388 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AEIZ77FKFHZIJ526HXBEO23S6FTN3ANCNFSM4WWKMFFA> .

gthompsnWRF · 2021-02-09T17:09:33Z

The idea is to commit both yours before we create the V4.3 branch.

OK, thanks for clarifying. I read your comment as having to wait until a future time when we unify codes in a new module_cloud_fraction (or similar). Just being sure that is not what you meant, thanks.

dudhia · 2021-02-18T17:01:23Z

I am going ahead and merging this one so that the MADWRF code can be compared for clashes. I see at least one area where they overlap in the Registry changes.

TYPE: new feature KEYWORDS: cloud fraction, cloud initialization SOURCE: Greg Thompson (JCSDA) DESCRIPTION OF CHANGES: Problem: Cold start simulations without microphysics species take significant time to 'spin-up' the clouds causing the radiation to be way off balance at the start of most model simulations. This greatly affects the first few hours up to 12 hours or more. Solution: The setting of initial cloud fraction and associated cloud water/ice mixing ratio gives a big advantage to 'spin up' the clouds by assuming that clouds can be created from conditions of very high relative humidity. The cloud fraction scheme of icloud=3 is a viable solution to making initial clouds along with the setting of the water vapor mixing ratio to the point of saturation at initialization in association with the inserted cloud fields. LIST OF MODIFIED FILES: M Registry/Registry.EM_COMMON M dyn_em/module_initialize_real.F M dyn_em/depend.dyn_em M run/README.namelist TESTS CONDUCTED: Numerous simulations in case studies plus real-time simulations with 3km spacing (nearly HRRR-CONUS) in support of FAA's ICICLE field project (2019Jan28-2019Mar08). Simulations by P. Jimenez in testing WRF-Solar-V2 Jenkins is pass. RELEASE NOTE: In combination with updates to the cloud fraction scheme (icloud=3), a new feature permits the usage of the cloud fraction scheme as part of real.exe, specifically designed to use with "cold start" simulations to reduce the spin-up problem of clouds and associated radiation.* add in possibility to initialize clouds in cold-start using icloud3 option Co-authored-by: Dave Gill <[email protected]>

add in possibility to initialize clouds in cold-start using icloud3 o…

19790ec

…ption

gthompsnWRF requested a review from a team as a code owner January 28, 2021 03:05

davegill added bug PHYSICS Develop Branch and removed Develop Branch PHYSICS bug labels Jan 28, 2021

davegill added PHYSICS New Feature Develop Branch labels Jan 28, 2021

davegill changed the title ~~Possibility to insert initial clouds into cold-start (PR replacement of #1323)~~ Possibility to insert initial clouds into cold-start Jan 28, 2021

davegill reviewed Jan 28, 2021

View reviewed changes

dyn_em/module_initialize_real.F Outdated Show resolved Hide resolved

remove redundant IF-test for namelist setting

7853fc9

Update Registry.EM_COMMON

0f9f7b7

dudhia reviewed Jan 28, 2021

View reviewed changes

dyn_em/module_initialize_real.F Outdated Show resolved Hide resolved

specify *only* for the required subroutine in radiation_driver, cal_c…

d6ed26c

…ldfra3

Merge branch 'insert_init_cloud2' of https://github.com/gthompsnWRF/WRF…

b930273

… into insert_init_cloud2

dudhia approved these changes Feb 11, 2021

View reviewed changes

davegill added the Infrastructure: OK label Feb 12, 2021

dudhia merged commit 6c4aef8 into wrf-model:develop Feb 18, 2021

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

Possibility to insert initial clouds into cold-start #1388

Possibility to insert initial clouds into cold-start #1388

gthompsnWRF commented Jan 28, 2021 •

edited by davegill

Loading

gthompsnWRF commented Jan 28, 2021

davegill commented Jan 28, 2021

davegill commented Jan 28, 2021

davegill commented Jan 28, 2021

davegill commented Jan 28, 2021

gthompsnWRF commented Jan 28, 2021

weiwangncar commented Jan 28, 2021

davegill commented Jan 28, 2021

dudhia commented Jan 28, 2021

gthompsnWRF commented Jan 28, 2021

davegill commented Jan 28, 2021

gthompsnWRF commented Jan 28, 2021

davegill commented Jan 28, 2021

gthompsnWRF commented Jan 28, 2021

weiwangncar commented Jan 28, 2021

davegill commented Feb 9, 2021

dudhia commented Feb 9, 2021 via email

gthompsnWRF commented Feb 9, 2021 •

edited

Loading

dudhia commented Feb 9, 2021 via email

gthompsnWRF commented Feb 9, 2021

dudhia commented Feb 18, 2021

Possibility to insert initial clouds into cold-start #1388

Possibility to insert initial clouds into cold-start #1388

Conversation

gthompsnWRF commented Jan 28, 2021 • edited by davegill Loading

gthompsnWRF commented Jan 28, 2021

davegill commented Jan 28, 2021

davegill commented Jan 28, 2021

davegill commented Jan 28, 2021

davegill commented Jan 28, 2021

gthompsnWRF commented Jan 28, 2021

weiwangncar commented Jan 28, 2021

davegill commented Jan 28, 2021

dudhia commented Jan 28, 2021

gthompsnWRF commented Jan 28, 2021

davegill commented Jan 28, 2021

gthompsnWRF commented Jan 28, 2021

davegill commented Jan 28, 2021

gthompsnWRF commented Jan 28, 2021

weiwangncar commented Jan 28, 2021

davegill commented Feb 9, 2021

dudhia commented Feb 9, 2021 via email

gthompsnWRF commented Feb 9, 2021 • edited Loading

dudhia commented Feb 9, 2021 via email

gthompsnWRF commented Feb 9, 2021

dudhia commented Feb 18, 2021

gthompsnWRF commented Jan 28, 2021 •

edited by davegill

Loading

gthompsnWRF commented Feb 9, 2021 •

edited

Loading