Changeset 3275 for trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90

Timestamp:

Mar 20, 2024, 3:05:14 PM (8 months ago)

Author:

afalco

Message:

Pluto PCM:
Changed _vap to _gas;
Included surfprop.F90;
callcorrk includes methane
AF

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• . (modified) (1 prop)
• physiq_mod.F90 (modified) (21 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90

@@ -300 +300 @@
       real omega(ngrid,nlayer)            ! omega velocity (Pa/s, >0 when downward)
 
-      integer i,l,ig,ierr,iq,nw,isoil,iesp, igcm_generic_vap, igcm_generic_ice
-      logical call_ice_vap_generic ! to call only one time the ice/vap pair of a tracer
+      integer i,l,ig,ierr,iq,nw,isoil,iesp, igcm_generic_gas, igcm_generic_ice
+      logical call_ice_gas_generic ! to call only one time the ice/vap pair of a tracer
 
       real zls                       ! Solar longitude (radians).
@@ -638 +638 @@
          qsurf_hist(:,:)=qsurf(:,:)
 
-!!         call WriteField_phy("post_qsurf_hist_qsurf",qsurf(1:ngrid,igcm_h2o_vap),1)
+!!         call WriteField_phy("post_qsurf_hist_qsurf",qsurf(1:ngrid,igcm_h2o_gas),1)
 
 !        Initialize variable for dynamical heating and zonal wind tendency diagnostic
@@ -681 +681 @@
          endif
 
-!!         call WriteField_phy("post_physdem_qsurf",qsurf(1:ngrid,igcm_h2o_vap),1)
+!!         call WriteField_phy("post_physdem_qsurf",qsurf(1:ngrid,igcm_h2o_gas),1)
 #endif
          if (corrk) then
@@ -702 +702 @@
          endif
 
-!!         call WriteField_phy("post_corrk_firstcall_qsurf",qsurf(1:ngrid,igcm_h2o_vap),1)
+!!         call WriteField_phy("post_corrk_firstcall_qsurf",qsurf(1:ngrid,igcm_h2o_gas),1)
          ! XIOS outputs
 #ifdef CPP_XIOS
@@ -711 +711 @@
 #endif
 
-!!         call WriteField_phy("post_xios_qsurf",qsurf(1:ngrid,igcm_h2o_vap),1)
+!!         call WriteField_phy("post_xios_qsurf",qsurf(1:ngrid,igcm_h2o_gas),1)
 
          write(*,*) "physiq: end of firstcall"
       endif ! end of 'firstcall'
 
-!!      call WriteField_phy("post_firstcall_qsurf",qsurf(1:ngrid,igcm_h2o_vap),1)
-!!      call writediagfi(ngrid,"firstcall_post_qsurf"," "," ",2,qsurf(1:ngrid,igcm_h2o_vap))
+!!      call WriteField_phy("post_firstcall_qsurf",qsurf(1:ngrid,igcm_h2o_gas),1)
+!!      call writediagfi(ngrid,"firstcall_post_qsurf"," "," ",2,qsurf(1:ngrid,igcm_h2o_gas))
 
       if (check_physics_inputs) then
@@ -882 +882 @@
             ! Eclipse incoming sunlight !AF24: removed
 
-!!            call writediagfi(ngrid,"corrk_pre_dqsurf"," "," ",2,dqsurf(1:ngrid,igcm_h2o_vap))
-!!            call writediagfi(ngrid,"corrk_pre_qsurf"," "," ",2,qsurf(1:ngrid,igcm_h2o_vap))
+!!            call writediagfi(ngrid,"corrk_pre_dqsurf"," "," ",2,dqsurf(1:ngrid,igcm_h2o_gas))
+!!            call writediagfi(ngrid,"corrk_pre_qsurf"," "," ",2,qsurf(1:ngrid,igcm_h2o_gas))
 
 
@@ -900 +900 @@
                   do iq=1,nq
 
-                     call generic_tracer_index(nq,iq,igcm_generic_vap,igcm_generic_ice,call_ice_vap_generic)
-
-                     if (call_ice_vap_generic) then ! to call only one time the ice/vap pair of a tracer
+                     call generic_tracer_index(nq,iq,igcm_generic_gas,igcm_generic_ice,call_ice_gas_generic)
+
+                     if (call_ice_gas_generic) then ! to call only one time the ice/vap pair of a tracer
 
                         epsi_generic=constants_epsi_generic(iq)
 
-                        muvar(1:ngrid,1:nlayer)=mugaz/(1.e0+(1.e0/epsi_generic-1.e0)*pq(1:ngrid,1:nlayer,igcm_generic_vap))
-                        muvar(1:ngrid,nlayer+1)=mugaz/(1.e0+(1.e0/epsi_generic-1.e0)*pq(1:ngrid,nlayer,igcm_generic_vap))
+                        muvar(1:ngrid,1:nlayer)=mugaz/(1.e0+(1.e0/epsi_generic-1.e0)*pq(1:ngrid,1:nlayer,igcm_generic_gas))
+                        muvar(1:ngrid,nlayer+1)=mugaz/(1.e0+(1.e0/epsi_generic-1.e0)*pq(1:ngrid,nlayer,igcm_generic_gas))
 
                      endif
@@ -931 +931 @@
                 call callcorrk(ngrid,nlayer,pq,nq,qsurf,                          &
                               albedo,albedo_equivalent,emis,mu0,pplev,pplay,pt,   &
-                              tsurf,fract,dist_star,aerosol,muvar,                &
+                              zzlay,tsurf,fract,dist_star,aerosol,muvar,                &
                               zdtlw,zdtsw,fluxsurf_lw,fluxsurf_sw,                &
                               fluxsurfabs_sw,fluxtop_lw,                          &
@@ -1039 +1039 @@
       endif ! of if (callrad)
 
-!!      call writediagfi(ngrid,"vdifc_pre_dqsurf"," "," ",2,dqsurf(1:ngrid,igcm_h2o_vap))
-!!      call writediagfi(ngrid,"vdifc_pre_qsurf"," "," ",2,qsurf(1:ngrid,igcm_h2o_vap))
+!!      call writediagfi(ngrid,"vdifc_pre_dqsurf"," "," ",2,dqsurf(1:ngrid,igcm_h2o_gas))
+!!      call writediagfi(ngrid,"vdifc_pre_qsurf"," "," ",2,qsurf(1:ngrid,igcm_h2o_gas))
 
 
@@ -1137 +1137 @@
          ! endif
 
-!!         call writediagfi(ngrid,"vdifc_post_zdqsdif"," "," ",2,zdqsdif(1:ngrid,igcm_h2o_vap))
+!!         call writediagfi(ngrid,"vdifc_post_zdqsdif"," "," ",2,zdqsdif(1:ngrid,igcm_h2o_gas))
 
          if (tracer) then
@@ -1144 +1144 @@
          end if ! of if (tracer)
 
-!!         call writediagfi(ngrid,"vdifc_post_dqsurf"," "," ",2,dqsurf(1:ngrid,igcm_h2o_vap))
-!!         call writediagfi(ngrid,"vdifc_post_qsurf"," "," ",2,qsurf(1:ngrid,igcm_h2o_vap))
+!!         call writediagfi(ngrid,"vdifc_post_dqsurf"," "," ",2,dqsurf(1:ngrid,igcm_h2o_gas))
+!!         call writediagfi(ngrid,"vdifc_post_qsurf"," "," ",2,qsurf(1:ngrid,igcm_h2o_gas))
 
          ! test energy conservation
@@ -1264 +1264 @@
          dqsurf(1:ngrid,igcm_n2) = dqsurf(1:ngrid,igcm_n2) + zdqsc(1:ngrid,igcm_n2)
 
-!!         call writediagfi(ngrid,"condense_n2_post_dqsurf"," "," ",2,dqsurf(1:ngrid,igcm_h2o_vap))
-!!         call writediagfi(ngrid,"condense_n2_post_qsurf"," "," ",2,qsurf(1:ngrid,igcm_h2o_vap))
+!!         call writediagfi(ngrid,"condense_n2_post_dqsurf"," "," ",2,dqsurf(1:ngrid,igcm_h2o_gas))
+!!         call writediagfi(ngrid,"condense_n2_post_qsurf"," "," ",2,qsurf(1:ngrid,igcm_h2o_gas))
 
          ! test energy conservation
@@ -1458 +1458 @@
                cloudfrac(:,:)=0.0
                do iq=1,nq
-                  call generic_tracer_index(nq,iq,igcm_generic_vap,igcm_generic_ice,call_ice_vap_generic)
-                  if (call_ice_vap_generic) then ! to call only one time the ice/vap pair of a tracer
+                  call generic_tracer_index(nq,iq,igcm_generic_gas,igcm_generic_ice,call_ice_gas_generic)
+                  if (call_ice_gas_generic) then ! to call only one time the ice/vap pair of a tracer
                      do l = 1, nlayer
                         do ig=1,ngrid
@@ -1491 +1491 @@
             dqsurf(1:ngrid,1:nq) = dqsurf(1:ngrid,1:nq) + zdqssed(1:ngrid,1:nq)
 
-!!            call writediagfi(ngrid,"callsedim_post_dqsurf"," "," ",2,dqsurf(1:ngrid,igcm_h2o_vap))
+!!            call writediagfi(ngrid,"callsedim_post_dqsurf"," "," ",2,dqsurf(1:ngrid,igcm_h2o_gas))
 
             ! ! Test water conservation !AF24: removed
@@ -1497 +1497 @@
          end if ! end of 'sedimentation'
 
-!!         call writediagfi(ngrid,"mass_redist_pre_dqsurf"," "," ",2,dqsurf(1:ngrid,igcm_h2o_vap))
-!!         call writediagfi(ngrid,"mass_redist_pre_qsurf"," "," ",2,qsurf(1:ngrid,igcm_h2o_vap))
+!!         call writediagfi(ngrid,"mass_redist_pre_dqsurf"," "," ",2,dqsurf(1:ngrid,igcm_h2o_gas))
+!!         call writediagfi(ngrid,"mass_redist_pre_qsurf"," "," ",2,qsurf(1:ngrid,igcm_h2o_gas))
 
   ! ---------------
@@ -1509 +1509 @@
             zdmassmr(1:ngrid,1:nlayer) = mass(1:ngrid,1:nlayer) * 0
             !    (   zdqevap(1:ngrid,1:nlayer)                          &
-            !    !   + zdqrain(1:ngrid,1:nlayer,igcm_h2o_vap)             &
-            !    !   + dqmoist(1:ngrid,1:nlayer,igcm_h2o_vap)             &
+            !    !   + zdqrain(1:ngrid,1:nlayer,igcm_h2o_gas)             &
+            !    !   + dqmoist(1:ngrid,1:nlayer,igcm_h2o_gas)             &
             !      + dqvaplscale(1:ngrid,1:nlayer) )
 
@@ -1536 +1536 @@
          endif
 
-!         call writediagfi(ngrid,"mass_redistribution_post_dqsurf"," "," ",2,dqsurf(1:ngrid,igcm_h2o_vap))
-
-!!         call writediagfi(ngrid,"slab_pre_dqsurf"," "," ",2,dqsurf(1:ngrid,igcm_h2o_vap))
-!!         call writediagfi(ngrid,"slab_pre_qsurf"," "," ",2,qsurf(1:ngrid,igcm_h2o_vap))
+!         call writediagfi(ngrid,"mass_redistribution_post_dqsurf"," "," ",2,dqsurf(1:ngrid,igcm_h2o_gas))
+
+!!         call writediagfi(ngrid,"slab_pre_dqsurf"," "," ",2,dqsurf(1:ngrid,igcm_h2o_gas))
+!!         call writediagfi(ngrid,"slab_pre_qsurf"," "," ",2,qsurf(1:ngrid,igcm_h2o_gas))
 
 
@@ -1794 +1794 @@
          do iq=1,nq
 
-            call generic_tracer_index(nq,iq,igcm_generic_vap,igcm_generic_ice,call_ice_vap_generic)
-
-            if (call_ice_vap_generic) then ! to call only one time the ice/vap pair of a tracer
+            call generic_tracer_index(nq,iq,igcm_generic_gas,igcm_generic_ice,call_ice_gas_generic)
+
+            if (call_ice_gas_generic) then ! to call only one time the ice/vap pair of a tracer
 
                do l = 1, nlayer
                   do ig=1,ngrid
                      call Psat_generic(zt(ig,l),pplay(ig,l),metallicity,psat_tmp_generic,qsat_generic(ig,l,iq))
-                     RH_generic(ig,l,iq) = zq(ig,l,igcm_generic_vap) / qsat_generic(ig,l,iq)
+                     RH_generic(ig,l,iq) = zq(ig,l,igcm_generic_gas) / qsat_generic(ig,l,iq)
                   enddo
                enddo
@@ -2364 +2364 @@
 
       do iq=1,nq
-         call generic_tracer_index(nq,iq,igcm_generic_vap,igcm_generic_ice,call_ice_vap_generic)
-
-         if (call_ice_vap_generic) then ! to call only one time the ice/vap pair of a tracer
+         call generic_tracer_index(nq,iq,igcm_generic_gas,igcm_generic_ice,call_ice_gas_generic)
+
+         if (call_ice_gas_generic) then ! to call only one time the ice/vap pair of a tracer
 
             reffrad_generic_zeros_for_wrf(:,:) = 1.
@@ -2373 +2373 @@
             comm_TOTCLOUDFRAC(1:ngrid) = totcloudfrac(1:ngrid) !??????
             comm_SURFRAIN(1:ngrid) = zdqsrain_generic(1:ngrid,iq)
-            comm_DQVAP(1:ngrid,1:nlayer) = pdq(1:ngrid,1:nlayer,igcm_generic_vap)
+            comm_DQVAP(1:ngrid,1:nlayer) = pdq(1:ngrid,1:nlayer,igcm_generic_gas)
             comm_DQICE(1:ngrid,1:nlayer)=pdq(1:ngrid,1:nlayer,igcm_generic_ice)
             ! comm_H2OICE_REFF(1:ngrid,1:nlayer) = reffrad_generic_zeros_for_wrf(1:ngrid,1:nlayer) ! for now zeros !
@@ -2417 +2417 @@
 !        do ig=1,ngrid
 !          if ((tracer).and.(water)) then
-!           pdq(ig,:,igcm_h2o_vap) = pdq(ig,:,igcm_h2o_vap) + lsf_dq(:)
+!           pdq(ig,:,igcm_h2o_gas) = pdq(ig,:,igcm_h2o_gas) + lsf_dq(:)
 !          endif
 !          pdt(ig,:) = pdt(ig,:) + lsf_dt(:)