Changeset 3102 for trunk/LMDZ.GENERIC
- Timestamp:
- Oct 25, 2023, 1:46:14 PM (14 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.GENERIC/libf/phystd/rain_generic.F90
r3043 r3102 298 298 ! JL22. Accounts for gravitational energy of falling precipitations (probably not to be used in the GCM 299 299 ! where the counterpart is not included in the dynamics.) 300 ! only for mesoscale simulations !!! 300 301 #endif 301 302 endif ! of if (precip_rate(i) .GT.0.) … … 307 308 308 309 if(precip_scheme_generic.eq.1)then 310 ! for now this is the only precipitation scheme working here and is therefore the default one (NC2023) 309 311 310 312 do i = 1, ngrid … … 322 324 323 325 elseif (precip_scheme_generic.ge.2) then 324 326 ! all the precipitation schemes below have not been tested and still have water-dependant routines 327 ! these schemes should be tested, improved ... 328 write(*,*) "Be careful, you have chosen a precip_scheme_generic equal to or greater than 2. For now it has not been tested." 325 329 do i = 1, ngrid 326 330 if (rneb(i,k).GT.0.0) then … … 412 416 enddo 413 417 414 endif ! if precip_scheme_generic =1418 endif ! if precip_scheme_generic 415 419 416 420 enddo ! of do k = nlayer, 1, -1 … … 422 426 call abort 423 427 endif 424 if (t(i,1) .LT. T_h2O_ice_liq) then 425 dqssnow_generic(i,i_ice_generic) = precip_rate(i) 426 dqsrain_generic(i,i_ice_generic) = 0.0 427 else 428 dqssnow_generic(i,i_ice_generic) = 0.0 429 dqsrain_generic(i,i_ice_generic) = precip_rate(i) ! liquid water = ice for now 430 endif 431 432 ! For now we force : 428 ! lines below should be uncommented and customized if you want to make a difference between rain and snow 429 !if (t(i,1) .LT. T_h2O_ice_liq) then 430 ! dqssnow_generic(i,i_ice_generic) = precip_rate(i) 431 ! dqsrain_generic(i,i_ice_generic) = 0.0 432 !else 433 ! dqssnow_generic(i,i_ice_generic) = 0.0 434 ! dqsrain_generic(i,i_ice_generic) = precip_rate(i) ! liquid water = ice for now 435 !endif 436 437 ! For now we make no distinction between rain and snow, and consider only rain: 433 438 dqsrain_generic(i,i_ice_generic) = precip_rate(i) 434 439 dqssnow_generic(i,i_ice_generic) = 0.0
Note: See TracChangeset
for help on using the changeset viewer.