Changeset 3444

Timestamp:

Sep 27, 2024, 6:54:41 PM (3 months ago)

Author:

afalco

Message:

Pluto PCM: reffrad not written if not initialized.
AF

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• callcorrk.F90 (modified) (4 diffs)
• callcorrk_pluto_mod.F90 (modified) (2 diffs)
• physiq_mod.F90 (modified) (1 diff)

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

@@ -493 +493 @@
          end do !iaer=1,naerkind.
          if (haze_radproffix) then
+            if (haze_radproffix) then
+                call haze_reffrad_fix(ngrid,nlayer,zzlay,&
+                    reffrad,nueffrad)
+            endif
+
             print*, 'haze_radproffix=T : fixed profile for haze rad'
          else
@@ -507 +512 @@
 
       if (aerohaze) then
-         if (haze_radproffix) then
-            call haze_reffrad_fix(ngrid,nlayer,zzlay,&
-                reffrad,nueffrad)
-          endif
 
          ! Get 3D aerosol optical properties.
@@ -576 +577 @@
 
 
-      ! AF24: for now only consider one aerosol (=haze) 
-      if (aerohaze) then 
+      ! AF24: for now only consider one aerosol (=haze)
+      if (aerohaze) then
         do iaer=1,naerkind
             ! Shortwave.
@@ -689 +690 @@
 !     Aerosol optical depths
 !-----------------------------------------------------------------------
-      if (aerohaze) then 
+      if (aerohaze) then
         do iaer=1,naerkind     ! a bug was here
             do k=0,nlayer-1

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

@@ -256 +256 @@
            end do !iaer=1,naerkind.
            if (haze_radproffix) then
+                if (haze_radproffix) then
+                    call haze_reffrad_fix(ngrid,nlayer,zzlay, &
+                        reffrad,nueffrad)
+                endif
               print*, 'haze_radproffix=T : fixed profile for haze rad'
            else
@@ -298 +302 @@
         ! uniform. Called only if aerohaze=true.
         !--------------------------------------------------
-        if (haze_radproffix) then
-          call haze_reffrad_fix(ngrid,nlayer,zzlay, &
-              reffrad,nueffrad)
-        endif
 
         call aeroptproperties(ngrid,nlayer,reffrad,nueffrad,         &

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

@@ -2409 +2409 @@
          !Pluto specific
          call writediagfi(ngrid,'n2_iceflux','n2_iceflux',"kg m^-2 s^-1",2,flusurf(1,igcm_n2) )
-         call writediagfi(ngrid,'haze_reff','haze_reff','m',3,reffrad(1,1,1))
+         if (aerohaze)then
+            call writediagfi(ngrid,'haze_reff','haze_reff','m',3,reffrad(1,1,1))
+         end if
          if (methane) then
             call writediagfi(ngrid,'ch4_iceflux','ch4_iceflux',&