Changeset 2727 for trunk/LMDZ.GENERIC

Timestamp:

Jun 21, 2022, 11:05:47 AM (3 years ago)

Author:

aslmd

Message:

added variable Generic tracer in callcorrk

File:

• trunk/LMDZ.GENERIC/libf/phystd/callcorrk.F90 (modified) (5 diffs)

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

@@ -33 +33 @@
       use callkeys_mod, only: varactive,diurnal,tracer,water,varfixed,satval,diagdtau,    &
                               kastprof,strictboundcorrk,specOLR,CLFvarying,               &
-                              tplanckmin,tplanckmax,global1d
+                              tplanckmin,tplanckmax,global1d,generic_condensation
       use optcv_mod, only: optcv
       use optci_mod, only: optci
       use recombin_corrk_mod, only: corrk_recombin, call_recombin
+      use generic_cloud_common_h, only: Psat_generic, epsi_generic
+      use generic_tracer_index_mod, only: generic_tracer_index
       implicit none
 
@@ -154 +156 @@
       REAL*8 albv(L_NSPECTV)                         ! Spectral Visible Albedo.
 
-      INTEGER ig,l,k,nw,iaer
+      INTEGER ig,l,k,nw,iaer,iq
 
       real*8,allocatable,save :: taugsurf(:,:)
       real*8,allocatable,save :: taugsurfi(:,:)
 !$OMP THREADPRIVATE(taugsurf,taugsurfi)
-      real*8 qvar(L_LEVELS)   ! Mixing ratio of variable component (mol/mol).
+      real*8 qvar(L_LEVELS)   ! Mixing ratio of variable component (mol/mol). index 1 is the top of the atmosphere, index L_LEVELS is the bottom
 
       ! Local aerosol optical properties for each column on RADIATIVE grid.
@@ -201 +203 @@
       ! local variable
       integer ok ! status (returned by NetCDF functions)
+
+      integer igcm_generic_vap, igcm_generic_ice! index of the vap and ice of generic_tracer
+      logical call_ice_vap_generic ! to call only one time the ice/vap pair of a tracer
+      real, save :: metallicity ! metallicity of planet --- is not used here, but necessary to call function Psat_generic
+!$OMP THREADPRIVATE(metallicity)
+      REAL, SAVE :: qvap_deep   ! deep mixing ratio of water vapor when simulating bottom less planets
+!$OMP THREADPRIVATE(qvap_deep)
 
 !===============================================================
@@ -394 +403 @@
          endif
 
-         if((igcm_h2o_vap.eq.0) .and. varactive)then
+         if((igcm_h2o_vap.eq.0) .and. varactive .and. water)then
             message='varactive in callcorrk but no h2o_vap tracer.'
             call abort_physic(subname,message,1)
+         endif
+
+         if(varfixed .and. generic_condensation .and. .not. water)then
+            write(*,*) "Deep water vapor mixing ratio ? (no effect if negative) "
+            qvap_deep=-1. ! default value
+            call getin_p("qvap_deep",qvap_deep)
+            write(*,*) " qvap_deep = ",qvap_deep
+
+            metallicity=0.0 ! default value --- is not used here but necessary to call function Psat_generic
+            call getin_p("metallicity",metallicity) ! --- is not used here but necessary to call function Psat_generic
          endif
 
@@ -648 +667 @@
 
 
-!-----------------------------------------------------------------------
-!     Water vapour (to be generalised for other gases eventually ...)
-!-----------------------------------------------------------------------
+      !-----------------------------------------------------------------------
+      !     Water vapour (to be generalised for other gases eventually ...)
+      !-----------------------------------------------------------------------
+            
+      if (water) then
+         if(varactive)then
+
+            i_var=igcm_h2o_vap
+            do l=1,nlayer
+               qvar(2*l)   = pq(ig,nlayer+1-l,i_var)
+               qvar(2*l+1) = pq(ig,nlayer+1-l,i_var)    
+               !JL13index   qvar(2*l+1) = (pq(ig,nlayer+1-l,i_var)+pq(ig,max(nlayer-l,1),i_var))/2    
+               !JL13index   Average approximation as for temperature...
+            end do
+            qvar(1)=qvar(2)
+
+         elseif(varfixed)then
+
+            do l=1,nlayer ! Here we will assign fixed water vapour profiles globally.
+               RH = satval * ((pplay(ig,l)/pplev(ig,1) - 0.02) / 0.98)
+               if(RH.lt.0.0) RH=0.0
+               
+               call Psat_water(pt(ig,l),pplay(ig,l),psat,qsat)
+
+               !pq_temp(l) = qsat      ! fully saturated everywhere
+               pq_temp(l) = RH * qsat ! ~realistic profile (e.g. 80% saturation at ground)
+            end do
+            
+            do l=1,nlayer
+               qvar(2*l)   = pq_temp(nlayer+1-l)
+               qvar(2*l+1) = (pq_temp(nlayer+1-l)+pq_temp(max(nlayer-l,1)))/2
+            end do
+            
+            qvar(1)=qvar(2)
+
+            ! Lowest layer of atmosphere
+            RH = satval * (1 - 0.02) / 0.98
+            if(RH.lt.0.0) RH=0.0
+
+            qvar(2*nlayer+1)= RH * qsat ! ~realistic profile (e.g. 80% saturation at ground)
+   
+         else
+            do k=1,L_LEVELS
+               qvar(k) = 1.0D-7
+            end do
+         end if ! varactive/varfixed
       
-      if(varactive)then
-
-         i_var=igcm_h2o_vap
-         do l=1,nlayer
-            qvar(2*l)   = pq(ig,nlayer+1-l,i_var)
-            qvar(2*l+1) = pq(ig,nlayer+1-l,i_var)    
-!JL13index            qvar(2*l+1) = (pq(ig,nlayer+1-l,i_var)+pq(ig,max(nlayer-l,1),i_var))/2    
-!JL13index            ! Average approximation as for temperature...
-         end do
-         qvar(1)=qvar(2)
-
-      elseif(varfixed)then
-
-         do l=1,nlayer ! Here we will assign fixed water vapour profiles globally.
-            RH = satval * ((pplay(ig,l)/pplev(ig,1) - 0.02) / 0.98)
-            if(RH.lt.0.0) RH=0.0
+      endif ! if (water)
+
+      !-----------------------------------------------------------------------
+      !  GCS (Generic Condensable Specie) Vapor
+      !  If you have GCS tracers and they are : variable & radiatively active
+      !
+      !  NC22
+      !-----------------------------------------------------------------------
+
+      if (generic_condensation .and. .not. water ) then
+
+         ! For now, only one GCS tracer can be both variable and radiatively active
+         ! If you set two GCS tracers, that are variable and radiatively active,
+         ! the last one in tracer.def will be chosen as the one that will be vadiatively active
+
+         do iq=1,nq
+
+            call generic_tracer_index(nq,iq,igcm_generic_vap,igcm_generic_ice,call_ice_vap_generic)
             
-            call Psat_water(pt(ig,l),pplay(ig,l),psat,qsat)
-
-            !pq_temp(l) = qsat      ! fully saturated everywhere
-            pq_temp(l) = RH * qsat ! ~realistic profile (e.g. 80% saturation at ground)
-         end do
+            if (call_ice_vap_generic) then ! to call only one time the ice/vap pair of a tracer
+
+               if(varactive)then
+
+                  i_var=igcm_generic_vap
+                  do l=1,nlayer
+                     qvar(2*l)   = pq(ig,nlayer+1-l,i_var)
+                     qvar(2*l+1) = pq(ig,nlayer+1-l,i_var)    
+                     !JL13index            qvar(2*l+1) = (pq(ig,nlayer+1-l,i_var)+pq(ig,max(nlayer-l,1),i_var))/2    
+                     !JL13index            ! Average approximation as for temperature...
+                  end do
+                  qvar(1)=qvar(2)
+
+               elseif(varfixed .and. (qvap_deep .ge. 0))then
+
+                  do l=1,nlayer ! Here we will assign fixed water vapour profiles globally.
+                                          
+                     call Psat_generic(pt(ig,l),pplay(ig,l),metallicity,psat,qsat)
+
+                     if (qsat .lt. qvap_deep) then
+                        pq_temp(l) = qsat      ! fully saturated everywhere
+                     else 
+                        pq_temp(l) = qvap_deep
+                     end if
+
+                  end do
+                  
+                  do l=1,nlayer
+                     qvar(2*l)   = pq_temp(nlayer+1-l)
+                     qvar(2*l+1) = (pq_temp(nlayer+1-l)+pq_temp(max(nlayer-l,1)))/2
+                  end do
+                  
+                  qvar(1)=qvar(2)
          
-         do l=1,nlayer
-            qvar(2*l)   = pq_temp(nlayer+1-l)
-            qvar(2*l+1) = (pq_temp(nlayer+1-l)+pq_temp(max(nlayer-l,1)))/2
-         end do
-         
-         qvar(1)=qvar(2)
-
-         ! Lowest layer of atmosphere
-         RH = satval * (1 - 0.02) / 0.98
-         if(RH.lt.0.0) RH=0.0
-
-         qvar(2*nlayer+1)= RH * qsat ! ~realistic profile (e.g. 80% saturation at ground)
- 
-      else
+               else
+                  do k=1,L_LEVELS
+                     qvar(k) = 1.0D-7
+                  end do
+               end if ! varactive/varfixed
+
+            endif
+
+         end do ! do iq=1,nq loop on tracers 
+
+      end if ! if (generic_condensation .and. .not. water )
+
+      !-----------------------------------------------------------------------
+      !  No Water vapor and No GCS (Generic Condensable Specie) vapor
+      !-----------------------------------------------------------------------
+
+      if (.not. generic_condensation .and. .not. water ) then
          do k=1,L_LEVELS
             qvar(k) = 1.0D-7
          end do
-      end if ! varactive/varfixed
+      end if ! if (.not. generic_condensation .and. .not. water )
+
 
       if(.not.kastprof)then
          ! IMPORTANT: Now convert from kg/kg to mol/mol.
          do k=1,L_LEVELS
-            qvar(k) = qvar(k)/(epsi+qvar(k)*(1.-epsi))
+            if (water) then
+               qvar(k) = qvar(k)/(epsi+qvar(k)*(1.-epsi))
+            endif
+            if (generic_condensation .and. .not. water) then
+               qvar(k) = qvar(k)/(epsi_generic+qvar(k)*(1.-epsi_generic))
+            endif
          end do
       end if