Changeset 2254 for LMDZ5/trunk/libf/phylmd/surf_ocean_mod.F90

Timestamp:

Apr 3, 2015, 12:16:36 PM (10 years ago)

Author:

fhourdin

Message:

Modification du calcul des flux air/mer
1) Introduction d'un facteur f_qsat_oce=0.98 devant qsat dans le calcul
de l'évaporation sur océan pour tenir compte de la moindre évaporation
de l'eau salée.
2) Introduction d'une différentiation entre z0 pour le sensible, z0h,
et le latent, z0q, imposé constant z0q=f_z0qh_oce*z0h

Modification of air/sea fluxes computation
1) Introduction of a correcting factor f_qsat_oce=0.98 on qsat
to account for the weaker evaporation of salty water.
2) Introduction of z0q=f_z0qh_oce*z0h

File:

• LMDZ5/trunk/libf/phylmd/surf_ocean_mod.F90 (modified) (14 diffs)

LMDZ5/trunk/libf/phylmd/surf_ocean_mod.F90

@@ -6 +6 @@
 CONTAINS
 !
-!****************************************************************************************
+!******************************************************************************
 !
   SUBROUTINE surf_ocean(rlon, rlat, swnet, lwnet, alb1, &
        windsp, rmu0, fder, tsurf_in, &
        itime, dtime, jour, knon, knindex, &
-       p1lay, cdragh, cdragm, precip_rain, precip_snow, temp_air, spechum, &
+       p1lay, z1lay, cdragh, cdragm, precip_rain, precip_snow, temp_air, spechum, &
        AcoefH, AcoefQ, BcoefH, BcoefQ, &
        AcoefU, AcoefV, BcoefU, BcoefV, &
@@ -38 +38 @@
 
 ! Input variables
-!****************************************************************************************
+!******************************************************************************
     INTEGER, INTENT(IN)                      :: itime, jour, knon
     INTEGER, DIMENSION(klon), INTENT(IN)     :: knindex
@@ -50 +50 @@
     REAL, DIMENSION(klon), INTENT(IN)        :: fder
     REAL, DIMENSION(klon), INTENT(IN)        :: tsurf_in
-    REAL, DIMENSION(klon), INTENT(IN)        :: p1lay
+    REAL, DIMENSION(klon), INTENT(IN)        :: p1lay,z1lay ! pression (Pa) et altitude (m) du premier niveau
     REAL, DIMENSION(klon), INTENT(IN)        :: cdragh
     REAL, DIMENSION(klon), INTENT(IN)        :: cdragm
@@ -63 +63 @@
 
 ! In/Output variables
-!****************************************************************************************
+!******************************************************************************
     REAL, DIMENSION(klon), INTENT(INOUT)     :: snow
     REAL, DIMENSION(klon), INTENT(INOUT)     :: qsurf
@@ -69 +69 @@
 
 ! Output variables
-!****************************************************************************************
+!******************************************************************************
     REAL, DIMENSION(klon), INTENT(OUT)       :: z0m, z0h
 !albedo SB >>>
@@ -84 +84 @@
 
 ! Local variables
-!****************************************************************************************
+!******************************************************************************
     INTEGER               :: i, k
     REAL                  :: tmp
@@ -90 +90 @@
     REAL, DIMENSION(klon) :: alb_eau
     REAL, DIMENSION(klon) :: radsol
+    REAL, DIMENSION(klon) :: cdragq ! Cdrag pour l'evaporation
 
 ! End definition
-!****************************************************************************************
-
-
-!****************************************************************************************
+!******************************************************************************
+
+
+!******************************************************************************
 ! Calculate total net radiance at surface
 !
-!****************************************************************************************
+!******************************************************************************
     radsol(:) = 0.0
     radsol(1:knon) = swnet(1:knon) + lwnet(1:knon)
 
-!****************************************************************************************
+!******************************************************************************
+! Cdragq computed from cdrag
+! The difference comes only from a factor (f_z0qh_oce) on z0, so that
+! it can be computed inside surf_ocean
+! More complicated appraches may require the propagation through
+! pbl_surface of an independant cdragq variable.
+!******************************************************************************
+
+    IF ( f_z0qh_oce .ne. 1.) THEN
+       cdragq(:)=cdragh(:)*                                      &
+       log(z1lay(:)/z0h(:))/log(z1lay(:)/(f_z0qh_oce*z0h(:)))
+    ELSE
+       cdragq(:)=cdragh(:)
+    ENDIF
+
+!******************************************************************************
 ! Switch according to type of ocean (couple, slab or forced)
-!****************************************************************************************
+!******************************************************************************
     SELECT CASE(type_ocean)
     CASE('couple')
@@ -111 +127 @@
             windsp, fder, & 
             itime, dtime, knon, knindex, &
-            p1lay, cdragh, cdragm, precip_rain, precip_snow,temp_air,spechum,& 
+            p1lay, cdragh, cdragq, cdragm, precip_rain, precip_snow,temp_air,spechum,& 
             AcoefH, AcoefQ, BcoefH, BcoefQ, &
             AcoefU, AcoefV, BcoefU, BcoefV, &
@@ -122 +138 @@
        CALL ocean_slab_noice( &
             itime, dtime, jour, knon, knindex, &
-            p1lay, cdragh, cdragm, precip_rain, precip_snow, temp_air, spechum,&
+            p1lay, cdragh, cdragq, cdragm, precip_rain, precip_snow, temp_air, spechum,&
             AcoefH, AcoefQ, BcoefH, BcoefQ, &
             AcoefU, AcoefV, BcoefU, BcoefV, &
@@ -133 +149 @@
        CALL ocean_forced_noice( &
             itime, dtime, jour, knon, knindex, &
-            p1lay, cdragh, cdragm, precip_rain, precip_snow, &
+            p1lay, cdragh, cdragq, cdragm, precip_rain, precip_snow, &
             temp_air, spechum, &
             AcoefH, AcoefQ, BcoefH, BcoefQ, &
@@ -143 +159 @@
     END SELECT
 
-!****************************************************************************************
+!******************************************************************************
 ! fcodron: compute lmt_bils  forced case (same as wfbils_oce / 1.-contfracatm)
-!****************************************************************************************
+!******************************************************************************
     IF (type_ocean.NE.'slab') THEN
         lmt_bils(:)=0.
@@ -154 +170 @@
     END IF
 
-!****************************************************************************************
+!******************************************************************************
 ! Calculate albedo
-!
-!****************************************************************************************
+!******************************************************************************
 !albedo SB >>>
-
-
   if(iflag_albedo==1)then
     call ocean_albedo(knon,rmu0,knindex,windsp,SFRWL,alb_dir_new,alb_dif_new)
@@ -180 +193 @@
 !albedo SB <<<
 
-!****************************************************************************************
+!******************************************************************************
 ! Calculate the rugosity
-!
-!****************************************************************************************
+!******************************************************************************
 IF (iflag_z0_oce==0) THEN
     DO i = 1, knon
@@ -197 +209 @@
 ENDIF
 !
-!****************************************************************************************
-!    
+!******************************************************************************
   END SUBROUTINE surf_ocean
-!
-!****************************************************************************************
+!******************************************************************************
 !
 END MODULE surf_ocean_mod