Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phylmd/phys_state_var_mod.F90

    r2240 r2243  
    2424      REAL, ALLOCATABLE, SAVE :: ftsol(:,:)
    2525!$OMP THREADPRIVATE(ftsol)
     26      REAL,ALLOCATABLE,SAVE :: qsol(:),fevap(:,:),z0m(:,:),z0h(:,:),agesno(:,:)
     27!$OMP THREADPRIVATE(qsol,fevap,z0m,z0h,agesno)
    2628!      character(len=6), SAVE :: ocean
    2729!!!!!!$OMP THREADPRIVATE(ocean)
     
    418420      ALLOCATE(pctsrf(klon,nbsrf))
    419421      ALLOCATE(ftsol(klon,nbsrf))
     422      ALLOCATE(qsol(klon),fevap(klon,nbsrf))
     423      ALLOCATE(z0m(klon,nbsrf+1),z0h(klon,nbsrf+1),agesno(klon,nbsrf))
    420424      ALLOCATE(falb1(klon,nbsrf))
    421425      ALLOCATE(falb2(klon,nbsrf))
     
    589593
    590594      deallocate(rlat, rlon, pctsrf, ftsol, falb1, falb2)
     595      deallocate(qsol,fevap,z0m,z0h,agesno)
    591596      deallocate(rain_fall, snow_fall, solsw, sollw, radsol)
    592597      deallocate(zmea, zstd, zsig, zgam)
Note: See TracChangeset for help on using the changeset viewer.