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

Timestamp:

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

Author:

tbertrand

Message:

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

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• physiq_mod.F90 (modified) (2 diffs)
• surfini.F (modified) (1 diff)
• turbdiff_mod.F90 (modified) (1 diff)

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)

trunk/LMDZ.PLUTO/libf/phypluto/surfini.F

@@ -44 +44 @@
       ! Step 2 : We get the bare ground albedo from the start files.
       DO ig=1,ngrid
-         albedo_bareground(ig)=albedodat(ig)
+         albedo_bareground(ig)=0.1 ! TB24 albedodat(ig)
          DO nw=1,L_NSPECTV
-            albedo(ig,nw)=albedo_bareground(ig)
+            albedo(ig,nw)=0.1 !albedo_bareground(ig)
          ENDDO
       ENDDO

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

@@ -138 +138 @@
 
       IF (firstcall) THEN
+         ivap=1 !TB24
+         iliq=0
+         iliq_surf=0
+         iice_surf=0 ! simply to make the code legible                    
          sensibFlux(:)=0.
          firstcall=.false.