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

Timestamp:

Mar 24, 2015, 2:28:51 PM (9 years ago)

Author:

fhourdin

Message:

Revisite de la formule des flux de surface
(en priorité sur l'océan) en tenant compte des bourrasques de
vent et de la différence entre les hauteurs de rugosités pour
la quantité de mouvement, l'enthalpie et éventuellement l'humidité.

Etape 2 :

• Séparation des z0 pour la quantité de mouvement et l'enthalpie.

rugs (ou frugs, rugos, yrugos ...) disparait au profit de z0m, z0h.
Les variables qui étaient à la fois dans pbl_surface_init et

• dans l'interface de pbl_surface sont suprimées de pbl_surface_init.

On travaille directement pour ces variables (evap, z0, qsol, agesno)
avec les versions de phys_state_var_mod (qui étaient
précédemment dans phys_local_var_mod

• Nouveaux paramètres de contrôle :
  • iflag_z0_oce (par défaut 0, et seule option active jusque là)
  • z0m_seaice_omp, z0h_seaice_omp, comme leur nom l'indique (utilisées dans surf_landice
  • z0min appliqué sur z0m et z0h dans pbl_surface
• Introduction des fonction phyeta0_get et phyetat0_srf pour lire

les conditions de initiales dans startphy.
Du coup une seule ligne suffit pour lire et contrôler d'éventuels
problèmes.

• Pour la variable fxrugs, elle est remplacée par z0m(:,nbsrf+1)

Ce choix déjà utilisé pour d'autres variables pourrait être
systématiser pour alléger l'interface de pbl_surface_mod.

• Dans les sorties, les variables rugs* ont été remplacées par

des z0m* et z0h*

• Nettoyage des anciens alb1/alb2 dans les lectures/écritures

des états de redémarrage (et dans pbl_surface_mod.F90).

File:

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

LMDZ5/trunk/libf/phylmd/surf_ocean_mod.F90

@@ -9 +9 @@
 !
   SUBROUTINE surf_ocean(rlon, rlat, swnet, lwnet, alb1, &
-       rugos, windsp, rmu0, fder, tsurf_in, &
+       windsp, rmu0, fder, tsurf_in, &
        itime, dtime, jour, knon, knindex, &
        p1lay, cdragh, cdragm, precip_rain, precip_snow, temp_air, spechum, &
@@ -16 +16 @@
        ps, u1, v1, gustiness, rugoro, pctsrf, &
        snow, qsurf, agesno, &
-!albedo SB >>>
-!      z0_new, alb1_new, alb2_new, evap, fluxsens, fluxlat, &
-       z0_new, SFRWL, alb_dir_new, alb_dif_new, evap, fluxsens, fluxlat, & 
-!albedo SB <<<
+       z0m, z0h, SFRWL, alb_dir_new, alb_dif_new, evap, fluxsens, fluxlat, & 
        tsurf_new, dflux_s, dflux_l, lmt_bils, &
        flux_u1, flux_v1)
@@ -49 +46 @@
     REAL, DIMENSION(klon), INTENT(IN)        :: lwnet  ! net longwave radiation at surface  
     REAL, DIMENSION(klon), INTENT(IN)        :: alb1   ! albedo in visible SW interval
-    REAL, DIMENSION(klon), INTENT(IN)        :: rugos
     REAL, DIMENSION(klon), INTENT(IN)        :: windsp
     REAL, DIMENSION(klon), INTENT(IN)        :: rmu0  
@@ -74 +70 @@
 ! Output variables
 !****************************************************************************************
-    REAL, DIMENSION(klon), INTENT(OUT)       :: z0_new
+    REAL, DIMENSION(klon), INTENT(OUT)       :: z0m, z0h
 !albedo SB >>>
 !    REAL, DIMENSION(klon), INTENT(OUT)       :: alb1_new  ! new albedo in visible SW interval
@@ -188 +184 @@
 !
 !****************************************************************************************
+IF (iflag_z0_oce==0) THEN
     DO i = 1, knon
        tmp = MAX(cepdu2,u1(i)**2+v1(i)**2)
-       z0_new(i) = 0.018*cdragm(i) * (u1(i)**2+v1(i)**2)/RG  &
+       z0m(i) = 0.018*cdragm(i) * (u1(i)**2+v1(i)**2)/RG  &
             +  0.11*14e-6 / SQRT(cdragm(i) * tmp)
-       z0_new(i) = MAX(1.5e-05,z0_new(i))
+       z0m(i) = MAX(1.5e-05,z0m(i))
     ENDDO   
+    z0h(1:knon)=z0m(1:knon) ! En attendant mieux
+
+ELSE
+       STOP'Alina, au boulot :)'
+ENDIF
 !
 !****************************************************************************************