Changeset 2898 for trunk/LMDZ.GENERIC/libf/phystd/dyn1d

Timestamp:

Feb 21, 2023, 9:08:49 AM (21 months ago)

Author:

emillour

Message:

Generic PCM:
Some cleanup in aerosol_mod.F90 : integrate "iniaerosol" in the module.
Also adapted kcm1d to follow up on recent code updates.
EM

File:

• trunk/LMDZ.GENERIC/libf/phystd/dyn1d/kcm1d.F90 (modified) (10 diffs)

trunk/LMDZ.GENERIC/libf/phystd/dyn1d/kcm1d.F90

@@ -10 +10 @@
   use callkeys_mod, only: pceil, tstrat, tracer, global1d
   use inifis_mod, only: inifis
+  use aerosol_mod, only: iniaerosol
+  use callcorrk_mod, only: callcorrk
   use comcstfi_mod
   use mod_grid_phy_lmdz, only : regular_lonlat
@@ -56 +58 @@
   real,allocatable :: vmr(:,:) ! tracer mixing ratio [mol/mol]
   real,allocatable :: qsurf(:)        ! tracer surface budget [kg/kg] ????
-  real psurf,psurf_n,tsurf      
-
-  real emis, albedo, albedo_equivalent
+  real psurf,psurf_n,tsurf(1)      
+
+  real emis(1), albedo(1), albedo_equivalent(1)
   real albedo_wv(1,L_NSPECTV)
 
@@ -65 +67 @@
   real dtsw(llm) ! heating rate (K/s) due to SW
   real dtlw(llm) ! heating rate (K/s) due to LW
-  real fluxsurf_lw   ! incident LW flux to surf (W/m2)
-  real fluxtop_lw    ! outgoing LW flux to space (W/m2)
-  real fluxsurf_sw   ! incident SW flux to surf (W/m2)
-  real fluxsurfabs_sw ! absorbed SW flux by the surf (W/m2)
-  real fluxabs_sw    ! SW flux absorbed by planet (W/m2)
-  real fluxtop_dn    ! incident top of atmosphere SW flux (W/m2)
+  real fluxsurf_lw(1)   ! incident LW flux to surf (W/m2)
+  real fluxtop_lw(1)    ! outgoing LW flux to space (W/m2)
+  real fluxsurf_sw(1)   ! incident SW flux to surf (W/m2)
+  real fluxsurfabs_sw(1) ! absorbed SW flux by the surf (W/m2)
+  real fluxabs_sw(1)    ! SW flux absorbed by planet (W/m2)
+  real fluxtop_dn(1)    ! incident top of atmosphere SW flux (W/m2)
 
   ! not used
-  real cloudfrac(llm)
-  real totcloudfrac
-  real tau_col
+  real cloudfrac(1,llm)
+  real totcloudfrac(1)
+  real tau_col(1)
 
   real dTstrat
@@ -81 +83 @@
   real OLR_nu(1,L_NSPECTI)
   real OSR_nu(1,L_NSPECTV)
+  real GSR_nu(1,L_NSPECTV)
+  real int_dtaui(1,llm,L_NSPECTI)
+  real int_dtauv(1,llm,L_NSPECTV)
   real Eatmtot
 
@@ -100 +105 @@
   pi=2.E+0*asin(1.E+0)
 
-  cloudfrac(:)  = 0.0
-  totcloudfrac  = 0.0
+  cloudfrac(1,:)  = 0.0
+  totcloudfrac(1)  = 0.0
 
 
@@ -254 +259 @@
   write(*,*) " albedo = ",albedo
   do iw=1,L_NSPECTV
-     albedo_wv(1,iw)=albedo
+     albedo_wv(1,iw)=albedo(1)
   enddo
 
@@ -360 +365 @@
          
                
-     call iniaerosol()
+     call iniaerosol
      
 
@@ -372 +377 @@
           dtlw,dtsw,fluxsurf_lw,fluxsurf_sw,              & 
           fluxsurfabs_sw,fluxtop_lw,                      &
-          fluxabs_sw,fluxtop_dn,OLR_nu,OSR_nu,tau_col,    &
-          cloudfrac,totcloudfrac,.false.,firstcall,lastcall)
+          fluxabs_sw,fluxtop_dn,OLR_nu,OSR_nu,GSR_nu,     &
+          int_dtaui,int_dtauv,                            &
+          tau_col,cloudfrac,totcloudfrac,                 &
+          .false.,firstcall,lastcall)
 
      !write(*,*) 'BASE 3'
@@ -397 +404 @@
      ! skin temperature (gray approx.) using analytic pure H2 expression
      !Tstrat  = (fluxabs_sw/(2*sigma))**0.25 ! skin temperature (gray approx.)
-     Tstrat  = (fluxtop_lw/(2*sigma))**0.25 ! skin temperature (gray approx.)
+     Tstrat  = (fluxtop_lw(1)/(2*sigma))**0.25 ! skin temperature (gray approx.)
      dTstrat = dTstrat-Tstrat
 
@@ -422 +429 @@
        tsurf,fract,dist_star,aerosol,muvar,                    &
        dtlw,dtsw,fluxsurf_lw,fluxsurf_sw,fluxsurfabs_sw,       &
-       fluxtop_lw, fluxabs_sw,fluxtop_dn,OLR_nu,OSR_nu,        &
-       tau_col,                                                &
-       cloudfrac,totcloudfrac,.false.,firstcall,lastcall)
+       fluxtop_lw, fluxabs_sw,fluxtop_dn,OLR_nu,OSR_nu,GSR_nu, &
+       int_dtaui,int_dtauv,                                    &
+       tau_col,cloudfrac,totcloudfrac,                         &
+       .false.,firstcall,lastcall)