Changeset 3982

Timestamp:

Dec 2, 2025, 11:54:50 AM (18 hours ago)

Author:

debatzbr

Message:

Pluto PCM: Add optical/HR diagnostics.

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• callcorrk.F90 (modified) (3 diffs)
• physiq_mod.F90 (modified) (1 diff)

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

@@ -44 +44 @@
       use recombin_corrk_mod, only: corrk_recombin, call_recombin
       use planetwide_mod, only: planetwide_maxval, planetwide_minval
-      use radcommon_h, only: wavev,wavei
+      use radcommon_h, only: wavev,wavei,bwnv,bwni
       use mp2m_diagnostics
       implicit none
@@ -219 +219 @@
       real :: vmrch4(ngrid,nlayer)              ! vmr ch4 from vmrch4_proffix
 
+      LOGICAL diagrad_rates ! Diagnostics
       REAL dtlw_nu(nlayer,L_NSPECTI) ! heating rate (K/s) due to LW in spectral bands
       REAL dtsw_nu(nlayer,L_NSPECTV) ! heating rate (K/s) due to SW in spectral bands
@@ -1233 +1234 @@
          endif
 
-      endif
+         !  Heating rates vs altitude in a .txt file
+         diagrad_rates=.true.
+         if(diagrad_rates)then
+            open(120,file='heating_rates.out')
+            write(120,*) "Pressure - Alt - HR tot - Rates (wavel SW)"
+            do l=1,nlayer
+              write(120,*) pplay(1,l),zzlay(1,l),dtsw(1,l),dtsw_nu(l,:)
+            enddo
+            close(120)
+
+            open(121,file='cooling_rates.out')
+            write(121,*) "Pressure - Alt - CR tot - Rates (wavel LW)"
+            do l=1,nlayer
+              write(121,*) pplay(1,l),zzlay(1,l),dtlw(1,l),dtlw_nu(l,:)
+            enddo
+            close(121)
+
+            open(122,file='bands.out')
+            write(122,*) "wavel - bands boundaries (microns)"
+            do nw=1,L_NSPECTV
+              write(122,*) wavev(nw),1.e4/bwnv(nw+1),1.e4/bwnv(nw)
+            enddo
+            do nw=1,L_NSPECTI
+              write(122,*) wavei(nw),1.e4/bwni(nw+1),1.e4/bwni(nw)
+            enddo
+            close(122)
+
+            open(123,file='c2h2_rates.out')
+            write(123,*) "Pressure - Alt - CR c2h2"
+            do l=1,nlayer
+              write(123,*) pplay(1,l),zzlay(1,l),dtlw_hcn_c2h2(1,l)
+            enddo
+            close(123)
+
+            open(124,file='opt_depth_vi.out')
+            write(124,*) "Pressure - Alt - dtau vi (wavel SW)"
+            do l=1,nlayer
+              write(124,*) pplay(1,l),zzlay(1,l),int_dtauv(1,l,:)
+            enddo
+            close(124)
+
+            open(125,file='opt_depth_ir.out')
+            write(125,*) "Pressure - Alt - dtau vi (wavel SW)"
+            do l=1,nlayer
+              write(125,*) pplay(1,l),zzlay(1,l),int_dtaui(1,l,:)
+            enddo
+            close(125)
+            endif
+      endif ! end of (lastcall) and (ngrid.eq.1)
 
       ! See physiq.F for explanations about CLFvarying. This is temporary.

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

@@ -2420 +2420 @@
       ! Warning this is exp(-dtau), I let you postproc with -log to have tau and k itself
       ! VI
-      call write_output('dtauv_4656nm','Layer optical thickness attenuation in VI band','',int_dtauv(:,nlayer:1:-1,2))  ! 4.656 um (17x28)
-      call write_output('dtauv_1181nm','Layer optical thickness attenuation in VI band','',int_dtauv(:,nlayer:1:-1,21)) ! 1.181 um (17x28)
-      call write_output('dtauv_700nm','Layer optical thickness attenuation in VI band','',int_dtauv(:,nlayer:1:-1,24))  ! 0.700 um (17x28)
-      call write_output('dtauv_185nm','Layer optical thickness attenuation in VI band','',int_dtauv(:,nlayer:1:-1,27))  ! 0.185 um (17x28)
-      call write_output('dtauv_118nm','Layer optical thickness attenuation in VI band','',int_dtauv(:,nlayer:1:-1,28))  ! 0.118 um (17x28)
+      call write_output('dtauv_4656nm','Layer optical thickness attenuation in VI band','',int_dtauv(:,nlayer:1:-1,2))  ! 4.656 um (28 VIS Bands)
+      call write_output('dtauv_1181nm','Layer optical thickness attenuation in VI band','',int_dtauv(:,nlayer:1:-1,21)) ! 1.181 um (28 VIS Bands)
+      call write_output('dtauv_700nm','Layer optical thickness attenuation in VI band','',int_dtauv(:,nlayer:1:-1,24))  ! 0.700 um (28 VIS Bands)
+      call write_output('dtauv_185nm','Layer optical thickness attenuation in VI band','',int_dtauv(:,nlayer:1:-1,27))  ! 0.185 um (28 VIS Bands)
+      call write_output('dtauv_118nm','Layer optical thickness attenuation in VI band','',int_dtauv(:,nlayer:1:-1,28))  ! 0.118 um (28 VIS Bands)
       ! IR
-      call write_output('dtaui_81250nm','Layer optical thickness attenuation in IR band','',int_dtaui(:,nlayer:1:-1,2)) ! 81.250 um (17x27)
-      call write_output('dtaui_3859nm','Layer optical thickness attenuation in IR band','',int_dtaui(:,nlayer:1:-1,16)) ! 3.859 um (17x27)
+      call write_output('dtaui_81250nm','Layer optical thickness attenuation in IR band','',int_dtaui(:,nlayer:1:-1,2)) ! 81.250 um (17 IR Bands)
+      call write_output('dtaui_3859nm','Layer optical thickness attenuation in IR band','',int_dtaui(:,nlayer:1:-1,16)) ! 3.859 um (17 IR Bands)
+      
+      !call write_output('dtaui_25250nm','Layer optical thickness attenuation in IR band','',int_dtaui(:,nlayer:1:-1,4)) ! 25.250 um (25 IR Bands)
+      !call write_output('dtaui_20800nm','Layer optical thickness attenuation in IR band','',int_dtaui(:,nlayer:1:-1,6)) ! 20.800 um (25 IR Bands)
+      !call write_output('dtaui_18000nm','Layer optical thickness attenuation in IR band','',int_dtaui(:,nlayer:1:-1,8)) ! 18.000 um (25 IR Bands)
+      !call write_output('dtaui_15050nm','Layer optical thickness attenuation in IR band','',int_dtaui(:,nlayer:1:-1,10)) ! 15.050 um (25 IR Bands)
+      
       if (callmufi) then
          ! Aerosol optical thickness