Changeset 3674

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

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• hazecloud.F90 (modified) (3 diffs)
• physiq_mod.F90 (modified) (5 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

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

@@ -518 +518 @@
          zdtsw(:,:) = 0.0
          zdtlw(:,:) = 0.0
+         zdqc(:,:,:)=0.
+         zdqsc(:,:)=0.
 
 !        Initialize tracer names, indexes and properties.
@@ -1270 +1272 @@
                                          zdqch4fast(1:ngrid)
             zdtsurf(1:ngrid)=zdtsurf(1:ngrid)+lw_ch4*zdqsch4fast(1:ngrid)/capcal(1:ngrid)
-            end if
+         end if
 !        ------------------------------------------------------------------
 !        CO surface sublimation and condensation in fast model (nogcm)
@@ -1361 +1363 @@
          endif
 
-         zdqc(:,:,:)=0.
-         zdqsc(:,:)=0.
          call condense_n2(ngrid,nlayer,nq,ptimestep,                    &
                            capcal,pplay,pplev,tsurf,pt,                 &
@@ -1497 +1497 @@
                zdqphot_ch4(:,:)    = 0.
                pdqmufi_prod(:,:,:) = 0.
-               call hazecloud(ngrid,nlayer,nq,ptimestep,                             &
+               call hazecloud(ngrid,nlayer,nq,ptimestep,zday,                        &
                               pplay,pplev,pq,pdq,dist_star,mu0,zfluxuv,pdqmufi_prod, &
                               zdqphot_prec,zdqphot_ch4,zdqconv_prec,declin)
@@ -1520 +1520 @@
                   zdqhaze(:,:,i_haze)=(profmmr(:,:)-pq(:,:,igcm_haze))/ptimestep
                else
-                  call hazecloud(ngrid,nlayer,nq,ptimestep,            &
+                  call hazecloud(ngrid,nlayer,nq,ptimestep,zday,       &
                      pplay,pplev,pq,pdq,dist_star,mu0,zfluxuv,zdqhaze, &
                      zdqphot_prec,zdqphot_ch4,zdqconv_prec,declin)