Changeset 3684

Timestamp:

Mar 14, 2025, 3:01:04 PM (3 months ago)

Author:

debatzbr

Message:

Cleaning for new diagnostics of optical thickness

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• callcorrk.F90 (modified) (2 diffs)
• callkeys_mod.F90 (modified) (1 diff)
• inifis_mod.F90 (modified) (1 diff)
• physiq_mod.F90 (modified) (2 diffs)

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

@@ -32 +32 @@
       use comcstfi_mod, only: pi, mugaz, cpp, r, g
       use callkeys_mod, only: diurnal,tracer,varfixed, &
-                              diagdtau,kastprof,strictboundcorrk,specOLR, &
+                              kastprof,strictboundcorrk,specOLR, &
                               tplanckmin,tplanckmax,global1d, &
                               optichaze,haze_radproffix,&
@@ -1061 +1061 @@
          end do
 
-      !  Optical thickness diagnostics (added by JVO)
-      if (diagdtau) then
-        do l=1,L_NLAYRAD
-          do nw=1,L_NSPECTV
-            int_dtauv(ig,l,nw) = 0.0d0
-             DO k=1,L_NGAUSS
-              ! Output exp(-tau) because gweight ponderates exp and not tau itself
-              int_dtauv(ig,l,nw)= int_dtauv(ig,l,nw) + exp(-dtauv(l,nw,k))*gweight(k)
-             ENDDO
-          enddo
-          do nw=1,L_NSPECTI
-           int_dtaui(ig,l,nw) = 0.0d0
-             DO k=1,L_NGAUSS
-              ! Output exp(-tau) because gweight ponderates exp and not tau itself
-              int_dtaui(ig,l,nw)= int_dtaui(ig,l,nw) + exp(-dtaui(l,nw,k))*gweight(k)
-             ENDDO
-          enddo
-        enddo
-      endif
+         ! Optical thickness diagnostics
+         ! Output exp(-tau) because gweight ponderates exp and not tau itself
+         int_dtauv(ig,:,:) = 0.0d0
+         int_dtaui(ig,:,:) = 0.0d0
+         do l=1,L_NLAYRAD
+            do nw=1,L_NSPECTV
+               do k=1,L_NGAUSS
+                  int_dtauv(ig,l,nw)= int_dtauv(ig,l,nw) + exp(-dtauv(l,nw,k))*gweight(k)
+               enddo
+            enddo
+            do nw=1,L_NSPECTI
+               do k=1,L_NGAUSS
+                  int_dtaui(ig,l,nw)= int_dtaui(ig,l,nw) + exp(-dtaui(l,nw,k))*gweight(k)
+               enddo
+            enddo
+         enddo
 
 ! **********************************************************

trunk/LMDZ.PLUTO/libf/phypluto/callkeys_mod.F90 ¶

@@ -22 +22 @@
       logical,save :: specOLR
       logical,save :: kastprof
-      logical,save :: diagdtau
-!$OMP THREADPRIVATE(enertest,nonideal,meanOLR,kastprof,diagdtau)
+!$OMP THREADPRIVATE(enertest,nonideal,meanOLR,kastprof)
       logical,save :: newtonian
       logical,save :: force_cpp

trunk/LMDZ.PLUTO/libf/phypluto/inifis_mod.F90 ¶

@@ -375 +375 @@
      if (is_master) write(*,*)trim(rname)//": specOLR = ",specOLR
 
-     if (is_master) write(*,*)trim(rname)//&
-       ": Output diagnostic optical thickness attenuation exp(-tau)?"
-     diagdtau=.false.
-     call getin_p("diagdtau",diagdtau)
-     if (is_master) write(*,*)trim(rname)//": diagdtau = ",diagdtau
-
      if (is_master) write(*,*)trim(rname)//": Operate in kastprof mode?"
      kastprof=.false.

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

@@ -48 +48 @@
                               callrad, callsoil, nosurf,                      &
                               callconduct,callmolvis,callmoldiff,             &
-                              corrk, diagdtau,                                &
+                              corrk,                                          &
                               diurnal, enertest, fat1au,                      &
                               icetstep, intheat, iradia, kastprof,            &
@@ -2190 +2190 @@
 
       ! Diagnostics of optical thickness
-      ! Warning this is exp(-tau), I let you postproc with -log to have tau itself - JVO 19
-      if (diagdtau.and..not.oldplutocorrk) then
-            do nw=1,L_NSPECTV
-               write(str2,'(i2.2)') nw
-               call write_output('dtauv'//str2,'Layer optical thickness attenuation in VI band '//str2,'',int_dtauv(:,nlayer:1:-1,nw))
-            enddo
-            do nw=1,L_NSPECTI
-               write(str2,'(i2.2)') nw
-               call write_output('dtaui'//str2,'Layer optical thickness attenuation in IR band '//str2,'',int_dtaui(:,nlayer:1:-1,nw))
-            enddo
-      endif
-
+      ! Warning this is exp(-tau), I let you postproc with -log to have tau itself
+      !do nw=1,L_NSPECTV
+      !   write(str2,'(i2.2)') nw
+      !   call write_output('dtauv'//str2,'Layer optical thickness attenuation in VI band '//str2,'',int_dtauv(:,nlayer:1:-1,nw))
+      !enddo
+      !do nw=1,L_NSPECTI
+      !   write(str2,'(i2.2)') nw
+      !   call write_output('dtaui'//str2,'Layer optical thickness attenuation in IR band '//str2,'',int_dtaui(:,nlayer:1:-1,nw))
+      !enddo
+      
       ! Temporary inclusions for heating diagnostics.
       if (.not.fast) then