Changeset 2891 for trunk/LMDZ.GENERIC


Ignore:
Timestamp:
Feb 7, 2023, 4:13:22 PM (23 months ago)
Author:
jleconte
Message:

potential energy in rain_generic added

correction : generic tracer is not always CH4

Location:
trunk/LMDZ.GENERIC/libf/phystd
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90

    r2890 r2891  
    16711671            zdqssnow_generic(1:ngrid,1:nq)    = 0.0
    16721672
    1673             call rain_generic(ngrid,nlayer,nq,ptimestep,pplev,pplay,pt,pdt,pq,pdq,            &
     1673            call rain_generic(ngrid,nlayer,nq,ptimestep,pplev,pplay,pphi,pt,pdt,pq,pdq,            &
    16741674                        zdtrain_generic,dq_rain_generic_vap,dq_rain_generic_cld,              &
    16751675                        zdqsrain_generic,zdqssnow_generic,reevap_precip_generic,cloudfrac)
     
    26922692         CALL send_xios_field('h2o_ice',zq(:,:,igcm_h2o_ice))
    26932693      ENDIF
    2694      
     2694
    26952695      CALL send_xios_field("ISR",fluxtop_dn)
    26962696      CALL send_xios_field("OLR",fluxtop_lw)
  • trunk/LMDZ.GENERIC/libf/phystd/rain_generic.F90

    r2782 r2891  
    1 subroutine rain_generic(ngrid,nlayer,nq,ptimestep,pplev,pplay,t,pdt,pq,pdq,d_t,dq_rain_generic_vap,dq_rain_generic_cld,dqsrain_generic,dqssnow_generic,reevap_precip,rneb)
     1subroutine rain_generic(ngrid,nlayer,nq,ptimestep,pplev,pplay,pphi,t,pdt,pq,pdq,d_t,dq_rain_generic_vap,dq_rain_generic_cld,dqsrain_generic,dqssnow_generic,reevap_precip,rneb)
    22
    33
     
    3333   real,intent(in) :: pplev(ngrid,nlayer+1) ! inter-layer pressure (Pa)
    3434   real,intent(in) :: pplay(ngrid,nlayer)   ! mid-layer pressure (Pa)
     35   real,intent(in) :: pphi(ngrid,nlayer)   ! mid-layer geopotential
    3536   real,intent(in) :: t(ngrid,nlayer) ! input temperature (K)
    3637   real,intent(in) :: pdt(ngrid,nlayer) ! input tendency on temperature (K/s)     
     
    283284                        precip_rate(i)  = precip_rate_tmp(i)
    284285                     end if
     286#ifdef MESOSCALE
     287                     d_t(i,k) = d_t(i,k)+(pphi(i,k+1)-pphi(i,k))*precip_rate(i)*ptimestep/(RCPD*dmass(i,k))
     288                     ! JL22. Accounts for gravitational energy of falling precipitations (probably not to be used in the GCM
     289                     !   where the counterpart is not included in the dynamics.)
     290#endif
    285291                  endif ! of if (precip_rate(i) .GT.0.)
    286292               enddo
Note: See TracChangeset for help on using the changeset viewer.