Changeset 2019 for trunk/LMDZ.GENERIC

Timestamp:

Oct 19, 2018, 11:38:16 AM (6 years ago)

Author:

mlefevre

Message:

Added outputs and large scale forcing (heating rate and water) for mesoscale generic model

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• comm_wrf.F90 (modified) (3 diffs)
• physiq_mod.F90 (modified) (3 diffs)

trunk/LMDZ.GENERIC/libf/phystd/comm_wrf.F90

@@ -32 +32 @@
   REAL,SAVE,ALLOCATABLE :: comm_H2OICE_REFF(:,:)
   REAL,SAVE,ALLOCATABLE :: comm_LATENT_HF(:)
+  REAL,SAVE,ALLOCATABLE :: comm_REEVAP(:)
+  REAL,SAVE,ALLOCATABLE :: comm_SURFRAIN(:)
 
 contains
@@ -66 +68 @@
   allocate(comm_H2OICE_REFF(ngrid,nlayer))
   allocate(comm_LATENT_HF(ngrid))
+  allocate(comm_REEVAP(ngrid))
+  allocate(comm_SURFRAIN(ngrid)
 
   end subroutine allocate_comm_wrf
@@ -98 +102 @@
   deallocate(comm_H2OICE_REFF)
   deallocate(comm_LATENT_HF)
+  deallocate(comm_REEVAP)
+  deallocate(comm_SURFRAIN)
 
   end subroutine deallocate_comm_wrf

trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90

@@ -392 +392 @@
       real :: flux_sens_lat(ngrid)
       real :: qsurfint(ngrid,nq)
-
+#ifdef MESOSCALE 
+      REAL :: lsf_dt(nlayer)
+      REAL :: lsf_dq(nlayer)
+#endif
 
 !==================================================================================================
@@ -2157 +2160 @@
       comm_HR_SW(1:ngrid,1:nlayer) = zdtsw(1:ngrid,1:nlayer)
       comm_HR_LW(1:ngrid,1:nlayer) = zdtlw(1:ngrid,1:nlayer)
-      comm_CLOUDFRAC(1:ngrid,1:nlayer)=cloudfrac(1:ngrid,1:nlayer)
-      comm_TOTCLOUDFRAC(1:ngrid)=totcloudfrac(1:ngrid)
-      comm_RAIN(1:ngrid,1:nlayer)=zdqrain(1:ngrid,1:nlayer,igcm_h2o_vap)
-      comm_SNOW(1:ngrid,1:nlayer)=zdqrain(1:ngrid,1:nlayer,igcm_h2o_ice)
-      comm_ALBEQ(1:ngrid)=albedo_equivalent(1:ngrid)
+      if ((tracer).and.(water)) then
+        comm_CLOUDFRAC(1:ngrid,1:nlayer)=cloudfrac(1:ngrid,1:nlayer)
+        comm_TOTCLOUDFRAC(1:ngrid)=totcloudfrac(1:ngrid)
+        comm_RAIN(1:ngrid,1:nlayer)=zdqrain(1:ngrid,1:nlayer,igcm_h2o_vap)
+        comm_SURFRAIN(1:ngrid)=zdqsrain(1:ngrid)
+        comm_SNOW(1:ngrid,1:nlayer)=zdqrain(1:ngrid,1:nlayer,igcm_h2o_ice)
+        comm_ALBEQ(1:ngrid)=albedo_equivalent(1:ngrid)
+        comm_H2OICE_REFF(1:ngrid,1:nlayer)=reffrad(1:ngrid,1:nlayer,iaero_h2o)
+        comm_REEVAP(1:ngrid)=reevap_precip(1:ngrid)
+      endif
       comm_FLUXTOP_DN(1:ngrid)=fluxtop_dn(1:ngrid)
       comm_FLUXABS_SW(1:ngrid)=fluxabs_sw(1:ngrid)
@@ -2169 +2177 @@
       comm_FLXGRD(1:ngrid)=fluxgrd(1:ngrid)
       comm_LSCEZ(1:ngrid,1:nlayer)=lscaledEz(1:ngrid,1:nlayer)
-      comm_H2OICE_REFF(1:ngrid,1:nlayer)=reffrad(1:ngrid,1:nlayer,iaero_h2o)
       sensibFlux(1:ngrid) = zflubid(1:ngrid) - capcal(1:ngrid)*zdtsdif(1:ngrid) !!! ????
+
+      if (turb_resolved) then
+        open(17,file='lsf.txt',form='formatted',status='old')
+        rewind(17)
+        DO l=1,nlayer
+          read(17,*) lsf_dt(l),lsf_dq(l)
+        ENDDO
+        close(17)
+        do ig=1,ngrid
+          if ((tracer).and.(water)) then
+           pdq(ig,:,igcm_h2o_vap) = pdq(ig,:,igcm_h2o_vap) + lsf_dq(:)
+          endif
+          pdt(ig,:) = pdt(ig,:) + lsf_dt(:)
+          comm_HR_DYN(ig,:) = lsf_dt(:)
+        enddo
+      endif
 #endif