Changeset 3198 for trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90

Timestamp:

Feb 2, 2024, 3:08:36 PM (10 months ago)

Author:

tbertrand

Message:

Cleaning newstart.F and adapting it to Pluto + small adjustments in the physics
TB

File:

• trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90 (modified) (2 diffs)

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

@@ -920 +920 @@
                print*,'and the surface albedo is taken equal to the first visible spectral value'
 
+               albedo_equivalent(1:ngrid)=albedo(1:ngrid,1)
                fluxsurfabs_sw(1:ngrid) = fluxtop_dn(1:ngrid)*(1.-albedo(1:ngrid,1))
+               ! TB24:
+               fluxabs_sw(1:ngrid)=fluxsurfabs_sw(1:ngrid)
                fluxrad_sky(1:ngrid)    = fluxsurfabs_sw(1:ngrid)
                fluxtop_lw(1:ngrid)  = emis(1:ngrid)*sigma*tsurf(1:ngrid)**4
@@ -1662 +1665 @@
       if(callrad)then
 
-         !call writediagfi(ngrid,"ALB","Surface albedo"," ",2,albedo_equivalent)
+         call writediagfi(ngrid,"ALB","Surface albedo"," ",2,albedo_equivalent)
          !call writediagfi(ngrid,"ALB_1st","First Band Surface albedo"," ",2,albedo(:,1))
          call writediagfi(ngrid,"ISR","incoming stellar rad.","W m-2",2,fluxtop_dn)