Changeset 3309 for trunk/LMDZ.GENERIC/libf

Timestamp:

Apr 22, 2024, 3:03:22 AM (6 months ago)

Author:

yjaziri

Message:

GENERIC PCM:

• Cosmetic + clarifying some variables and comments in chemistry
• add controle variable for actinique UV flux in photochemistry and output surface UV flux in diagspecUV.nc

YJ

Location:

trunk/LMDZ.GENERIC/libf

Files:

• aeronostd/calchim_asis.F90 (modified) (5 diffs)
• aeronostd/chimiedata_h.F90 (modified) (1 diff)
• aeronostd/iniwrite_spec.F (added)
• aeronostd/photochemistry_asis.F90 (modified) (1 diff)
• aeronostd/photolysis_mod.F90 (modified) (5 diffs)
• aeronostd/photolysis_online.F (modified) (4 diffs)
• aeronostd/writediagspecUV.F (added)
• phystd/callkeys_mod.F90 (modified) (1 diff)
• phystd/inifis_mod.F90 (modified) (1 diff)

trunk/LMDZ.GENERIC/libf/aeronostd/calchim_asis.F90

@@ -5 +5 @@
                          pu,pdu,pv,pdv,surfdust,surfice,icount,zdtchim)
 
-      use tracer_h, only: mmol, noms, nesp, is_chim
-
-      use conc_mod, only: mmean ! mean molecular mass of the atmosphere
+      use tracer_h,         only: mmol, noms, nesp, is_chim
+
+      use conc_mod,         only: mmean ! mean molecular mass of the atmosphere
       USE comcstfi_mod
       use callkeys_mod
@@ -14 +14 @@
       use chimiedata_h
       use radcommon_h,      only: glat
-      use wstats_mod, only: wstats
+      use wstats_mod,       only: wstats
 
       implicit none
@@ -256 +256 @@
 !        chemistry in lower atmosphere
 
-!         if (photochem) then
-
             call photochemistry_asis(nlayer,ngrid,                         &
                                 ig,lswitch,zycol,szacol,fractcol,ptimestep,&
@@ -338 +336 @@
 ! value of parameter 'output' to trigger writting of outputs
 ! is set above at the declaration of the variable.
-      if (photochem .and. output) then
+      if (output) then
 
             ! photoloysis
@@ -417 +415 @@
             endif
          endif ! end depos
+
+         if (jonline .and. output_writediagspecUV) then
+
+            !flux at the surface
+            call writediagspecUV(ngrid,"flux","flux(lon,lat,band)","photon.s-1.nm-1.cm-2",3,fluxUV(:,:,1))
+
+         endif ! end jonline
 
       end if ! of if (output)

trunk/LMDZ.GENERIC/libf/aeronostd/chimiedata_h.F90

@@ -56 +56 @@
       real, save, allocatable :: albedo_snow_chim(:)    ! albedo snow on chemistry wavelenght grid
       real, save, allocatable :: albedo_co2_ice_chim(:) ! albedo co2 ice on chemistry wavelenght grid
+      real, save, allocatable :: fluxUV(:,:,:)          ! total actinic flux at each ngrid, wavelength and altitude (photon.s-1.nm-1.cm-2)
 
 !--------------------------------------------

trunk/LMDZ.GENERIC/libf/aeronostd/photochemistry_asis.F90

@@ -137 +137 @@
      allocate(albedo_snow_chim(nw))
      allocate(albedo_co2_ice_chim(nw))
+     allocate(fluxUV(ngrid,nw,nlayer))
+     fluxUV(:,:,:) = 0.
     ! Step 1 : Initialisation of the Spectral Albedos.
      DO iw=1,nw

trunk/LMDZ.GENERIC/libf/aeronostd/photolysis_mod.F90

@@ -32 +32 @@
     ! solar flux
   
-    real, allocatable, save :: f(:)            !  solar flux (photon.s-1.nm-1.cm-2) at 1 au
+    real, allocatable, save :: fstar1AU(:)     ! stellar flux (photon.s-1.nm-1.cm-2) at 1 au
   
     ! cross-sections and quantum yields
@@ -69 +69 @@
     call gridw(nw,wl,wc,wu,mopt)
   
-    allocate(f(nw))
+    allocate(fstar1AU(nw))
   
     ! read and grid solar flux data
    
-    call rdsolarflux(nw,wl,wc,f)
+    call rdsolarflux(nw,wl,wc,fstar1AU)
   
     ! read maximum wavelength until photochemical heating is calculated
@@ -474 +474 @@
     !==============================================================================
    
-    subroutine rdsolarflux(nw,wl,wc,f)
+    subroutine rdsolarflux(nw,wl,wc,fstar1AU)
   
         !     Read and re-grid solar flux data.            
@@ -490 +490 @@
         !     output
   
-        real, dimension(nw) :: f  ! solar flux (photon.s-1.nm-1.cm-2)
+        real, dimension(nw) :: fstar1AU  ! stellar flux (photon.s-1.nm-1.cm-2)
   
         !     local 
@@ -568 +568 @@
         ! 5.039e11 = 1.e-4*1e-9/(hc = 6.62e-34*2.998e8)
   
-        ! f need to be in photon.s-1.nm-1.cm-2
+        ! fstar1AU need to be in photon.s-1.nm-1.cm-2
         DO iw = 1, nw-1
            !f(iw) = yg1(iw)*wc(iw)*5.039e11  ! If yg1 in w.m-2.nm-1
-           f(iw) = yg1(iw)
+           fstar1AU(iw) = yg1(iw)
         ENDDO
   

trunk/LMDZ.GENERIC/libf/aeronostd/photolysis_online.F

@@ -20 +20 @@
       use photolysis_mod
       use tracer_h
-      use chimiedata_h, only: indexchim
+      use chimiedata_h, only: indexchim, fluxUV
       use types_asis,   only: nb_phot_hv_max, nb_phot_max, jlabel,
      $                        reactions
@@ -45 +45 @@
       real (kind = 8), dimension(nlayer,nb_phot_max) :: e_phot      ! photolysis rates by energie (J.mol-1.s-1)
  
-!     solar flux at mars
-
-      real, dimension(nw) :: fmars                                  ! solar flux (w.m-2.nm-1)
+!     stellar flux at planet
+
+      real, dimension(nw) :: fplanet                                ! stellar flux (photon.s-1.nm-1.cm-2) at planet
       real :: factor
 
@@ -198 +198 @@
       call sphers(nlev,zalt,sza,dsdh,nid)
 
-!==== solar flux at mars
+!==== stellar flux at planet
 
       factor = (1./dist_sol)**2.
       do iw = 1,nw-1
-         fmars(iw) = f(iw)*factor
+         fplanet(iw) = fstar1AU(iw)*factor
       end do
 
@@ -229 +229 @@
 
          do ilay = 1,nlayer
-           saflux(ilay) = fmars(iw)*(fdir(ilay) + fdn(ilay) + fup(ilay))
+           saflux(ilay) = fplanet(iw)*
+     $                     (fdir(ilay) + fdn(ilay) + fup(ilay))
+           fluxUV(ig,iw,ilay) = saflux(ilay)
          end do
 

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

@@ -87 +87 @@
       logical,save :: depos
       logical,save :: haze
-!$OMP THREADPRIVATE(photochem,photoheat,jonline,depos)
+      logical,save :: output_writediagspecUV
+!$OMP THREADPRIVATE(photochem,photoheat,jonline,depos,output_writediagspecUV)
       logical,save :: calllott_nonoro
 !$OMP THREADPRIVATE(calllott_nonoro)

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

@@ -543 +543 @@
      if (is_master) write(*,*)trim(rname)//": haze = ",haze
 
+     if (is_master) write(*,*)trim(rname)//": Output photochemical surface UV flux ?"
+     output_writediagspecUV=.false. ! default value
+     call getin_p("output_writediagspecUV",output_writediagspecUV)
+     if (is_master) write(*,*)trim(rname)//": output_writediagspecUV = ",output_writediagspecUV
+
 ! Global1D mean and solar zenith angle