Changeset 3081 for trunk/LMDZ.GENERIC/libf

Timestamp:

Oct 11, 2023, 12:36:59 PM (13 months ago)

Author:

yjaziri

Message:

Generic PCM:

Move mu0 (cosinus of solar zenith angle) calculation out of callrad
This was a issue for photochemistry because mu0 was updated only iradia times

YJ

File:

• trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90 (modified) (1 diff)

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

@@ -903 +903 @@
 !---------------------------------
 
+      ! Compute local stellar zenith angles
+      ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+      if (tlocked) then
+      ! JL14 corrects tidally resonant (and inclined) cases. nres=omega_rot/omega_orb
+         ztim1=SIN(declin)
+         ztim2=COS(declin)*COS(zlss)
+         ztim3=COS(declin)*SIN(zlss)
+
+         call stelang(ngrid,sinlon,coslon,sinlat,coslat,    &
+                        ztim1,ztim2,ztim3,mu0,fract, flatten)
+
+      elseif (diurnal) then
+         ztim1=SIN(declin)
+         ztim2=COS(declin)*COS(2.*pi*(zday-.5))
+         ztim3=-COS(declin)*SIN(2.*pi*(zday-.5))
+
+         call stelang(ngrid,sinlon,coslon,sinlat,coslat,    &
+                        ztim1,ztim2,ztim3,mu0,fract, flatten)
+      else if(diurnal .eqv. .false.) then
+
+         call mucorr(ngrid,declin,latitude,mu0,fract,10000.,rad,flatten)
+         ! WARNING: this function appears not to work in 1D
+
+         if ((ngrid.eq.1).and.(global1d)) then ! Fixed zenith angle 'szangle' in 1D simulations w/ globally-averaged sunlight.
+            mu0 = cos(pi*szangle/180.0)
+         endif
+
+      endif 
+
       if (callrad) then
          if( mod(icount-1,iradia).eq.0.or.lastcall) then
-
-          ! Compute local stellar zenith angles
-          ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-            if (tlocked) then
-            ! JL14 corrects tidally resonant (and inclined) cases. nres=omega_rot/omega_orb
-               ztim1=SIN(declin)
-               ztim2=COS(declin)*COS(zlss)
-               ztim3=COS(declin)*SIN(zlss)
-
-               call stelang(ngrid,sinlon,coslon,sinlat,coslat,    &
-                            ztim1,ztim2,ztim3,mu0,fract, flatten)
-
-            elseif (diurnal) then
-               ztim1=SIN(declin)
-               ztim2=COS(declin)*COS(2.*pi*(zday-.5))
-               ztim3=-COS(declin)*SIN(2.*pi*(zday-.5))
-
-               call stelang(ngrid,sinlon,coslon,sinlat,coslat,    &
-                            ztim1,ztim2,ztim3,mu0,fract, flatten)
-            else if(diurnal .eqv. .false.) then
-
-               call mucorr(ngrid,declin,latitude,mu0,fract,10000.,rad,flatten)
-               ! WARNING: this function appears not to work in 1D
-
-               if ((ngrid.eq.1).and.(global1d)) then ! Fixed zenith angle 'szangle' in 1D simulations w/ globally-averaged sunlight.
-                  mu0 = cos(pi*szangle/180.0)
-                  !print*,'acosz=',mu0,', szangle=',szangle
-               endif
-
-            endif 
            
             ! Eclipse incoming sunlight (e.g. Saturn ring shadowing).