Changeset 2560

Timestamp:

Sep 7, 2021, 3:28:26 PM (3 years ago)

Author:

slebonnois

Message:

SL: Implementation of SW computation based on generic model. Switch between this new SW module or old module that reads R. Haus tables implemented with a key (solarchoice)

Location:

trunk/LMDZ.VENUS

Files:

• deftank/context_lmdz_physics-78niv.xml (modified) (1 diff)
• deftank/context_lmdz_physics-notrac.xml (modified) (1 diff)
• deftank/context_lmdz_physics.xml (modified) (1 diff)
• deftank/physiq-96x96x50-chemistry.def (modified) (1 diff)
• deftank/physiq-96x96x50-notrac.def (modified) (1 diff)
• deftank/physiq-96x96x78-chemistry.def (modified) (1 diff)
• deftank/physiq-96x96x90-chemistry.def (modified) (1 diff)
• deftank/physiq.def (modified) (1 diff)
• libf/phyvenus/aerave_new.F (added)
• libf/phyvenus/aeropacity.F90 (added)
• libf/phyvenus/aeroptproperties.F90 (added)
• libf/phyvenus/aerosol_mod.F90 (added)
• libf/phyvenus/ave_stelspec.F90 (added)
• libf/phyvenus/bilinearbig.F90 (added)
• libf/phyvenus/blackl.F (added)
• libf/phyvenus/calc_rayleigh.F90 (added)
• libf/phyvenus/clesphys.h (modified) (2 diffs)
• libf/phyvenus/conf_phys.F90 (modified) (2 diffs)
• libf/phyvenus/datafile_mod.F90 (added)
• libf/phyvenus/dsolver.F (added)
• libf/phyvenus/dtridgl.F (added)
• libf/phyvenus/gases_h.F90 (added)
• libf/phyvenus/gfluxv.F (added)
• libf/phyvenus/iniaerosol.F (added)
• libf/phyvenus/interpolateH2Ocont_CKD.F90 (added)
• libf/phyvenus/optcv.F90 (added)
• libf/phyvenus/radcommon_h.F90 (added)
• libf/phyvenus/radinc_h.F90 (added)
• libf/phyvenus/radlwsw.F (modified) (7 diffs)
• libf/phyvenus/setspv.F90 (added)
• libf/phyvenus/sfluxv.F (added)
• libf/phyvenus/su_aer_radii.F90 (added)
• libf/phyvenus/su_gases.F90 (added)
• libf/phyvenus/suaer_corrk.F90 (added)
• libf/phyvenus/sugas_corrk.F90 (added)
• libf/phyvenus/sw_venus_corrk.F90 (added)
• libf/phyvenus/tpindex.F (added)
• libf/phyvenus/write_histins.h (modified) (2 diffs)
• libf/phyvenus/write_histmth.h (modified) (2 diffs)

trunk/LMDZ.VENUS/deftank/context_lmdz_physics-78niv.xml

@@ -101 +101 @@
                    long_name="Vertical wind" 
                    unit="Pa/s" />
+           <field id="vitwz"
+                   long_name="Vertical wind (m/s)" 
+                   unit="m/s" />
             <field id="Kz"
                    long_name="vertical diffusion coef" 

trunk/LMDZ.VENUS/deftank/context_lmdz_physics-notrac.xml

@@ -101 +101 @@
                    long_name="Vertical wind" 
                    unit="Pa/s" />
+           <field id="vitwz"
+                   long_name="Vertical wind (m/s)" 
+                   unit="m/s" />
             <field id="Kz"
                    long_name="vertical diffusion coef" 

trunk/LMDZ.VENUS/deftank/context_lmdz_physics.xml

@@ -101 +101 @@
                    long_name="Vertical wind" 
                    unit="Pa/s" />
+           <field id="vitwz"
+                   long_name="Vertical wind (m/s)" 
+                   unit="m/s" />
             <field id="Kz"
                    long_name="vertical diffusion coef" 

trunk/LMDZ.VENUS/deftank/physiq-96x96x50-chemistry.def

@@ -73 +73 @@
 #
 iflag_ajs = 1
+
+# Solar radiation module
+##~~~~~~~~~~~~~~~~~~~~~~
+## 1 = Rainer Haus Tables
+## 2 = Generic solar module
 #
+solarchoice=1
+
 # Thermosphere and nlte parameters
 ##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

trunk/LMDZ.VENUS/deftank/physiq-96x96x50-notrac.def

@@ -73 +73 @@
 #
 iflag_ajs = 1
+
+# Solar radiation module
+##~~~~~~~~~~~~~~~~~~~~~~
+## 1 = Rainer Haus Tables
+## 2 = Generic solar module
 #
+solarchoice=1
+
 # Thermosphere and nlte parameters
 ##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

trunk/LMDZ.VENUS/deftank/physiq-96x96x78-chemistry.def

@@ -73 +73 @@
 #
 iflag_ajs = 1
+
+# Solar radiation module
+##~~~~~~~~~~~~~~~~~~~~~~
+## 1 = Rainer Haus Tables
+## 2 = Generic solar module
 #
+solarchoice=1
+
 # Thermosphere and nlte parameters
 ##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

trunk/LMDZ.VENUS/deftank/physiq-96x96x90-chemistry.def

@@ -73 +73 @@
 #
 iflag_ajs = 1
+
+# Solar radiation module
+##~~~~~~~~~~~~~~~~~~~~~~
+## 1 = Rainer Haus Tables
+## 2 = Generic solar module
 #
+solarchoice=1
+
 # Thermosphere and nlte parameters
 ##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

trunk/LMDZ.VENUS/deftank/physiq.def

@@ -73 +73 @@
 #
 iflag_ajs = 1
+
+# Solar radiation module
+##~~~~~~~~~~~~~~~~~~~~~~
+## 1 = Rainer Haus Tables
+## 2 = Generic solar module
 #
+solarchoice=1
+
 # Thermosphere and nlte parameters
 ##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

trunk/LMDZ.VENUS/libf/phyvenus/clesphys.h

@@ -19 +19 @@
        INTEGER nircorr, nltemodel, solvarmod
        INTEGER nb_mode
+       INTEGER solarchoice
        REAL    ecriphy
        REAL    z0, lmixmin
@@ -33 +34 @@
      &     iflag_con, iflag_ajs,                                        &
      &     lev_histins, lev_histday, lev_histmth, tr_scheme,            &
-     &     cl_scheme, nircorr, nltemodel, solvarmod, nb_mode
+     &     cl_scheme, nircorr, nltemodel, solvarmod, nb_mode,           &
+     &     solarchoice
 
        COMMON/clesphys_r/ ecriphy, z0, lmixmin,                         &

trunk/LMDZ.VENUS/libf/phyvenus/conf_phys.F90

@@ -388 +388 @@
   call getin('nb_mode',nb_mode)
   
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+! PARAMETER FOR SOLAR RADIATION
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+
+!
+!Config Key  = solarchoice
+!Config Desc = 
+!Config Def  = 1
+!Config Help = 
+!
+! 1 = RH Tables
+! 2 = Generic module
+  solarchoice = 1
+  call getin('solarchoice',solarchoice)
+
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 ! PARAMETER FOR NLTE PHYSICS
@@ -514 +529 @@
   write(lunout,*)' ok_sedim = ',ok_sedim
   write(lunout,*)' nb_mode = ',nb_mode
+  write(lunout,*)' solarchoice = ',solarchoice
   write(lunout,*)' callnlte = ',callnlte
   write(lunout,*)' nltemodel = ',nltemodel

trunk/LMDZ.VENUS/libf/phyvenus/radlwsw.F

@@ -33 +33 @@
      .      topsw,toplw,solsw,sollw,sollwdown,lwnet,swnet
       use write_field_phy
+      use radinc_h, only: ini_radinc_h
 
 #ifdef CPP_XIOS      
@@ -60 +61 @@
 !===========
       INTEGER k, kk, i, j, band
+      integer,save :: i_sw
+      integer,parameter :: subloop=100
 
       REAL   PPB(klev+1)
@@ -103 +106 @@
         call load_ksi(ksive)
 
+        i_sw = subloop
+c for sw_venus_corrk
+        if (solarchoice.eq.2) call ini_radinc_h(klev)
+        
       endif ! firstcall
 c-------------------------------------------
@@ -108 +115 @@
       DO k = 1, klev
        DO i = 1, klon
-         heat(i,k)=0.
+! Solar heating rates from generic => directly on klon
+! only once every subloop calls, because of high frequency needed by lw...
+         if (i_sw.eq.subloop) heat(i,k)=0.
          cool(i,k)=0.
        ENDDO
@@ -600 +609 @@
 c      CALL SW_venus_rh_1Dglobave(zrmu0,zfract,   ! pour moy globale
 c     S        PPB,temp, 
-      CALL SW_venus_rh(zrmu0,zfract,latdeg,
+c     S        zheat, 
+c     S        ztopsw,zsolsw,ZFSNET)
+
+cc SOLAR : RH TABLES
+       if (solarchoice.eq.1) then
+        CALL SW_venus_rh(zrmu0,zfract,latdeg,
      S        PPA,PPB,temp, 
      S        zheat, 
      S        ztopsw,zsolsw,ZFSNET)
+       endif ! solarchoice.eq.1
+
 c======================================================================
-         radsol(j) = zsolsw - zsollw  ! + vers bas
-         topsw(j) = ztopsw            ! + vers bas
+c lw into klon grid:
+         radsol(j) =  - zsollw        ! + vers bas
          toplw(j) = ztoplw            ! + vers haut
-         solsw(j) = zsolsw            ! + vers bas
          sollw(j) = -zsollw           ! + vers bas
          sollwdown(j) = zsollwdown    ! + vers bas
 
          DO k = 1, klev+1
-         lwnet  (j,k)   = ZFLNET(k)
-         swnet  (j,k)   = ZFSNET(k)
+           lwnet  (j,k)   = ZFLNET(k)    ! + vers haut
          ENDDO
 
+cc SOLAR : RH TABLES
+c sw into klon grid (if not using the generic routine) 
+         if (solarchoice.eq.1) then
+           radsol(j)= radsol(j)+zsolsw  ! + vers bas
+           topsw(j) = ztopsw            ! + vers bas
+           solsw(j) = zsolsw            ! + vers bas
+           DO k = 1, klev+1
+            swnet  (j,k)   = ZFSNET(k)   ! + vers bas
+           ENDDO
+         endif ! solarchoice.eq.1
+
 c
-C heat/cool with upper atmosphere 
+C cool with upper atmosphere 
 C
       IF(callnlte) THEN
          DO k = 1,nlaylte
-           heat(j,k) = zheat(k)
            cool(j,k) = zcool(k)
          ENDDO
@@ -628 +652 @@
        if (klev.gt.nlaylte) then
          do k = nlaylte+1,  klev
-           heat(j,k) = 0.
            cool(j,k) = 0.
          enddo
@@ -634 +657 @@
       ELSE
          DO k = 1, klev
-           heat(j,k) = zheat(k)
            cool(j,k) = zcool(k)
          ENDDO
       ENDIF ! callnlte
+
+cc SOLAR : RH TABLES
+C heat with upper atmosphere 
+      IF(solarchoice.eq.1) THEN
+       IF(callnlte) THEN
+         DO k = 1,nlaylte
+           heat(j,k) = zheat(k)
+         ENDDO
+c     Zero tendencies for any remaining layers between nlaylte and klev
+        if (klev.gt.nlaylte) then
+         do k = nlaylte+1,  klev
+           heat(j,k) = 0.
+         enddo
+        endif    
+       ELSE
+         DO k = 1, klev
+           heat(j,k) = zheat(k)
+         ENDDO
+       ENDIF ! callnlte
+      ENDIF ! solarchoice.eq.1
+
       ENDDO ! of DO j = 1, klon
 c+++++++ FIN BOUCLE SUR LA GRILLE +++++++++++++++++++++++++
+
+cc SOLAR : GENERICS
+! Solar heating rates from generic => directly on klon
+! only once every subloop calls, because of high frequency needed by lw...
+      if (solarchoice.eq.2) then
+
+       if (i_sw.eq.subloop) then
+
+        call sw_venus_corrk(klon,klev,rmu0,paprs,pplay,t,tsol,
+     .       fract,dist,heat,solsw,topsw,swnet,firstcall)
+
+        IF(callnlte) THEN
+c     Zero tendencies for any remaining layers between nlaylte and klev
+         if (klev.gt.nlaylte) then
+          do k = nlaylte+1,klev
+           heat(:,k) = 0.
+          enddo
+         endif    
+        ENDIF ! callnlte
+       
+        i_sw=0
+       endif ! i_sw=subloop
+
+       radsol(:)=radsol(:)+solsw(:)
+       i_sw=i_sw+1
+      endif ! solarchoice.eq.2
+
+!----------------------------------------------------------
 
 ! for tests: write output fields...

trunk/LMDZ.VENUS/libf/phyvenus/write_histins.h

@@ -46 +46 @@
 c     call histwrite_phy(nid_ins,.false.,"mang",itau_w,mang)
 c     call histwrite_phy(nid_ins,.false.,"Kz",itau_w,ycoefh)
-      call histwrite_phy(nid_ins,.false.,"mmean",itau_w,mmean)
-      call histwrite_phy(nid_ins,.false.,"rho",itau_w,rho)
+c     call histwrite_phy(nid_ins,.false.,"mmean",itau_w,mmean)
+c     call histwrite_phy(nid_ins,.false.,"rho",itau_w,rho)
 
 c plusieurs traceurs  !!!outputs in [vmr]
@@ -104 +104 @@
        ENDIF
              
-      call histwrite_phy(nid_ins,.false.,"WH2SO4",itau_w,WH2SO4)
+c     call histwrite_phy(nid_ins,.false.,"WH2SO4",itau_w,WH2SO4)
       
-      call histwrite_phy(nid_ins,.false.,"rho_droplet",itau_w,
-     & rho_droplet)
+c     call histwrite_phy(nid_ins,.false.,"rho_droplet",itau_w,
+c    & rho_droplet)
                 ENDIF
 

trunk/LMDZ.VENUS/libf/phyvenus/write_histmth.h

@@ -51 +51 @@
 c plusieurs traceurs  !!!outputs in [vmr]
        IF (iflag_trac.eq.1) THEN
-         DO iq=1,nqmax
+         DO iq=1,nqmax-nmicro
        call histwrite_phy(nid_mth,.false.,tname(iq),itau_w,qx(:,:,iq)
      &                   *mmean(:,:)/M_tr(iq))
+         ENDDO
+         DO iq=nqmax-nmicro+1,nqmax
+       call histwrite_phy(nid_mth,.false.,tname(iq),itau_w,qx(:,:,iq))
          ENDDO
        ENDIF
@@ -105 +108 @@
 c en K/s      
       call histwrite_phy(nid_mth,.false.,"dtswr",itau_w,dtsw)
-      call histwrite_phy(nid_mth,.false.,"dtswrNLTE",itau_w,d_t_nirco2)
-      call histwrite_phy(nid_mth,.false.,"dtswrDCrisp",itau_w,heat)
+c     call histwrite_phy(nid_mth,.false.,"dtswrNLTE",itau_w,d_t_nirco2)
+c     call histwrite_phy(nid_mth,.false.,"dtswrDCrisp",itau_w,heat)
 c en K/s      
       call histwrite_phy(nid_mth,.false.,"dtlwr",itau_w,dtlw)