Changeset 3312

Timestamp:

Apr 22, 2024, 6:44:48 PM (8 months ago)

Author:

yjaziri

Message:

Generic PCM:

Photochemistry: correction for using diurnal equal true in 1D
No correction factor with diurnal equal true for photolysis rate
+ writediagspecUV change name to flux_surf for clarity

YJ

Location:

trunk/LMDZ.GENERIC

Files:

• changelog.txt (modified) (1 diff)
• libf/aeronostd/calchim_asis.F90 (modified) (1 diff)
• libf/aeronostd/photochemistry_asis.F90 (modified) (1 diff)

trunk/LMDZ.GENERIC/changelog.txt

@@ -1916 +1916 @@
 will before open_startphy was doing both by create_restartphy
 at firstcall
+
+== 22/04/2024 == YJ
+Photochemistry: correction for using diurnal equal true in 1D
+No correction factor with diurnal equal true for photolysis rate
++ writediagspecUV change name to flux_surf for clarity

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

@@ -419 +419 @@
 
             !flux at the surface
-            call writediagspecUV(ngrid,"flux","flux(lon,lat,band)","photon.s-1.nm-1.cm-2",3,fluxUV(:,:,1))
+            call writediagspecUV(ngrid,"flux_surf","Surface flux(lon,lat,band)","photon.s-1.nm-1.cm-2",3,fluxUV(:,:,1))
 
          endif ! end jonline

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

@@ -169 +169 @@
       call photolysis_online(nlayer, alt, press, temp, zmmean, rm,   &
                              tau, sza, dist_sol, v_phot, e_phot, ig, ngrid, nreact)
-      if (ngrid.eq.1) then
-        do iphot = 1,nb_phot_hv_max
-          v_phot(:,iphot) = v_phot(:,iphot)* 0.25 / cos(sza*pi/180.) ! globally averaged = divide by 4
-          e_phot(:,iphot) = e_phot(:,iphot)* 0.25 / cos(sza*pi/180.) ! globally averaged = divide by 4
-        end do
-      elseif(diurnal .eqv. .false.) then
-        do iphot = 1,nb_phot_hv_max
-          v_phot(:,iphot) = v_phot(:,iphot) * fractcol
-          e_phot(:,iphot) = e_phot(:,iphot) * fractcol
-        end do
+      if (diurnal .eqv. .false.) then
+        if (ngrid.eq.1) then
+          do iphot = 1,nb_phot_hv_max
+            v_phot(:,iphot) = v_phot(:,iphot)* 0.25 / cos(sza*pi/180.) ! globally averaged = divide by 4
+            e_phot(:,iphot) = e_phot(:,iphot)* 0.25 / cos(sza*pi/180.) ! globally averaged = divide by 4
+          end do
+        else
+          do iphot = 1,nb_phot_hv_max
+            v_phot(:,iphot) = v_phot(:,iphot) * fractcol
+            e_phot(:,iphot) = e_phot(:,iphot) * fractcol
+          end do
+        endif
       endif
    else ! night