Changeset 3269 for trunk/LMDZ.GENERIC

Timestamp:

Mar 14, 2024, 2:03:31 PM (9 months ago)

Author:

mturbet

Message:

add net surface longwave heat flux diagnostics

Location:

trunk/LMDZ.GENERIC

Files:

• changelog.txt (modified) (1 diff)
• libf/phystd/physiq_mod.F90 (modified) (4 diffs)

trunk/LMDZ.GENERIC/changelog.txt

@@ -1886 +1886 @@
 - Preliminary cleanup+commenting in hydrol.F90
 - Changes in computation of snow albedo when ok_slab_ocean is activated
+- add diagnostics in physiq_mod.F90 to compute net surface longwave flux (required for intercomparison projects)

trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90

@@ -452 +452 @@
       real int_dtaui1(ngrid,nlayer,L_NSPECTI),int_dtauv1(ngrid,nlayer,L_NSPECTV) ! For optical thickness diagnostics.
       real tf, ntf
+      
+      real net_fluxsurf_lw(ngrid) ! net longwave flux at the surface (for diagnostics only)
 
       real nconsMAX, vdifcncons(ngrid), cadjncons(ngrid) ! Vdfic water conservation test. By RW
@@ -1122 +1124 @@
                GSR=0.0
                fluxrad_sky(1:ngrid)=emis(1:ngrid)*fluxsurf_lw(1:ngrid)+fluxsurfabs_sw(1:ngrid)
+               net_fluxsurf_lw(1:ngrid)=emis(1:ngrid)*fluxsurf_lw(1:ngrid)
 
                             !if(noradsurf)then ! no lower surface; SW flux just disappears
@@ -1176 +1179 @@
          zplanck(1:ngrid)=tsurf(1:ngrid)*tsurf(1:ngrid)
          zplanck(1:ngrid)=emis(1:ngrid)*sigma*zplanck(1:ngrid)*zplanck(1:ngrid)
+         net_fluxsurf_lw(1:ngrid)=net_fluxsurf_lw(1:ngrid)-zplanck(1:ngrid)
          fluxrad(1:ngrid)=fluxrad_sky(1:ngrid)-zplanck(1:ngrid)
          pdt(1:ngrid,1:nlayer)=pdt(1:ngrid,1:nlayer)+dtrad(1:ngrid,1:nlayer)
@@ -2534 +2538 @@
 !           call writediagfi(ngrid,"fluxsurfswcs","sw surface flux (cs).","W m-2",2,fluxsurf_sw1)
 !           call writediagfi(ngrid,"fluxsurflwcs","lw back radiation (cs).","W m-2",2,fluxsurf_lw1)
+!           call writediagfi(ngrid,"netfluxsurflw","lw net surface flux","W m-2",2,net_fluxsurf_lw)
 
          if(ok_slab_ocean) then