Changeset 3674 for trunk/LMDZ.PLUTO/libf/phypluto/hazecloud.F90

Timestamp:

Mar 10, 2025, 2:49:16 PM (3 months ago)

Author:

tbertrand

Message:

Pluto: small fix in hazecloud to get the anti-solar longitude when running in parallel
TB

File:

• trunk/LMDZ.PLUTO/libf/phypluto/hazecloud.F90 (modified) (3 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/hazecloud.F90

@@ -1 @@
-      subroutine hazecloud(ngrid,nlayer,nq,ptimestep, &
+      subroutine hazecloud(ngrid,nlayer,nq,ptimestep,zday, &
           pplay,pplev,pq,pdq,pdist_sol,mu0,pfluxuv,zdqhaze, &
           zdqphot_prec,zdqphot_ch4,zdqconv_prec,declin)
@@ -49 +49 @@
       DATA firstcall/.true./
       REAL ptimestep
+      REAL zday
       REAL pplay(ngrid,nlayer),pplev(ngrid,nlayer+1)
       REAL,INTENT(IN) :: pq(ngrid,nlayer,nq)
@@ -124 +125 @@
       ELSE
 !     1)  get longitude/latitude (rad) of anti-subsolar point (max de mu0 - 180)
-        longit=longitude((MAXLOC(mu0,DIM=1,MASK=mu0.GT.0.9)))     ! rad
-        IF (longit.GE.0) THEN
-           longit=longit-pi
-        ELSE
-           longit=longit+pi
+        longit=-(zday-int(zday))*2.*pi 
+        IF (longit.LE.-pi) THEN
+           longit=longit+2.*pi
         ENDIF
         latit=-declin                                        ! rad
-
 !     2) Define input parameter for the fit
         valmin=48.74e10    ! minimum value of ipm flux in ph/m2/s