Changeset 3940 for LMDZ6/trunk/libf/phylmdiso/surf_seaice_mod.F90

Timestamp:

Jun 15, 2021, 1:18:14 PM (3 years ago)

Author:

crisi

Message:

replace files by symbloic liks from phylmdiso towards phylmd.
Many files at once

File:

• LMDZ6/trunk/libf/phylmdiso/surf_seaice_mod.F90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmdiso/surf_seaice_mod.F90

@@ -45 +45 @@
     INCLUDE "dimsoil.h"
     INCLUDE "clesphys.h"
+
+    INCLUDE "YOMCST.h"
+    ! for rd and retv
 
 ! Input arguments
@@ -113 +116 @@
     REAL, DIMENSION(klon) :: alb1_new,alb2_new
 !albedo SB <<<
-!
+
+    real rhoa(knon) ! density of moist air  (kg / m3)
+
 ! End definitions
 !****************************************************************************************
@@ -124 +129 @@
     radsol(:) = 0.0
     radsol(1:knon) = swnet(1:knon) + lwnet(1:knon)
+
+    rhoa = PS(:KNON) / (Rd * temp_air(:knon) * (1. + retv * spechum(:knon)))
 
 !****************************************************************************************
@@ -142 +149 @@
             radsol, snow, qsurf, &
             alb1_new, alb2_new, evap, fluxsens, fluxlat, flux_u1, flux_v1, &
-            tsurf_new, dflux_s, dflux_l)
+            tsurf_new, dflux_s, dflux_l, rhoa)
        
     ELSE IF (type_ocean == 'slab'.AND.version_ocean=='sicINT') THEN
@@ -164 +171 @@
             radsol, snow, qsol, agesno, tsoil, &
             qsurf, alb1_new, alb2_new, evap, fluxsens, fluxlat, flux_u1, flux_v1, &
-            tsurf_new, dflux_s, dflux_l &
+            tsurf_new, dflux_s, dflux_l, rhoa &
 #ifdef ISO
             ,xtprecip_rain, xtprecip_snow, xtspechum,Roce, &