Changeset 4352 for LMDZ6/trunk


Ignore:
Timestamp:
Nov 20, 2022, 12:01:18 PM (3 years ago)
Author:
evignon
Message:

option pour conditionner l'activation des params SSO a un critere
sur la "planitude" de la maille afin de ne pas freiner de facon irrealiste
les vents sur les calottes
Travail de Valentin Wiener

Location:
LMDZ6/trunk/libf/phylmd
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmd/clesphys.h

    r4245 r4352  
    77!
    88!..include cles_phys.h
    9 !
     9
     10       ! threshold on zrel_mount to activate SSO schemes
     11       REAL zrel_mount_t
    1012       INTEGER iflag_cycle_diurne
    1113       LOGICAL soil_model,new_oliq,ok_orodr,ok_orolf
     
    122124     &     , co2_ppm0                                                   &
    123125!FC
    124      &     , Cd_frein                                                   &
     126     &     , Cd_frein, zrel_mount_t                                     &
    125127     &     , ecrit_LES                                                  &
    126128     &     , ecrit_ins, ecrit_hf, ecrit_day                             &
  • LMDZ6/trunk/libf/phylmd/conf_phys_m.F90

    r4250 r4352  
    226226    LOGICAL, SAVE :: ok_lic_cond_omp
    227227    !
     228    REAL, SAVE    :: zrel_mount_t_omp
    228229    INTEGER, SAVE :: iflag_cycle_diurne_omp
    229230    LOGICAL, SAVE :: soil_model_omp,new_oliq_omp
     
    900901    CALL getin('ok_orolf', ok_orolf_omp)
    901902
     903
     904    !Config  Key  =  zrel_mount_t
     905    !Config  Desc = zrel_mount_t
     906    !Config  Def  = 0.
     907    !Config  Help = Connais pas !
     908    zrel_mount_t_omp = 0.
     909    CALL getin('zrel_mount_t', zrel_mount_t_omp)
     910
     911
    902912    !Config  Key  = ok_limitvrai
    903913    !Config  Desc = Force la lecture de la bonne annee
     
    24262436    ok_orodr = ok_orodr_omp
    24272437    ok_orolf = ok_orolf_omp
     2438    zrel_mount_t=zrel_mount_t_omp
    24282439    ok_limitvrai = ok_limitvrai_omp
    24292440    nbapp_rad = nbapp_rad_omp
     
    28512862    WRITE(lunout,*) ' ok_orodr=',ok_orodr
    28522863    WRITE(lunout,*) ' ok_orolf=',ok_orolf
     2864    WRITE(lunout,*) ' zrel_mount_t=',zrel_mount_t
    28532865    WRITE(lunout,*) ' ok_limitvrai=',ok_limitvrai
    28542866    WRITE(lunout,*) ' nbapp_rad=',nbapp_rad
  • LMDZ6/trunk/libf/phylmd/physiq_mod.F90

    r4298 r4352  
    972972
    973973    REAL picefra(klon,klev)
     974    REAL zrel_mount(klon)
    974975    !IM cf. AM 081204 END
    975976    !
     
    45474548       DO i=1,klon
    45484549          itest(i)=0
    4549           !        IF ((zstd(i).gt.10.0)) THEN
    4550           IF (((zpic(i)-zmea(i)).GT.100.).AND.(zstd(i).GT.10.0)) THEN
     4550          zrel_mount(i)=zstd(i)/(max(zsig(i),1.E-8)*sqrt(cell_area(i)))
     4551          !zrel_mount: relative mountain height wrt relief explained by mean slope
     4552          ! -> condition on zrel_mount can deactivate the drag on tilted planar terrains
     4553          !    such as ice sheets (work by V. Wiener)
     4554          IF (((zpic(i)-zmea(i)).GT.100.).AND.(zstd(i).GT.10.0).AND.(zrel_mount(i).GE.zrel_mount_t)) THEN
    45514555             itest(i)=1
    45524556             igwd=igwd+1
     
    45974601       DO i=1,klon
    45984602          itest(i)=0
    4599           IF ((zpic(i)-zmea(i)).GT.100.) THEN
     4603          !zrel_mount: relative mountain height wrt relief explained by mean slope
     4604          ! -> condition on zrel_mount can deactivate the lifting on tilted planar terrains
     4605          !    such as ice sheets (work by V. Wiener)
     4606          zrel_mount(i)=zstd(i)/(max(zsig(i),1.E-8)*sqrt(cell_area(i)))
     4607          IF (((zpic(i)-zmea(i)).GT.100.).AND.(zrel_mount(i).GE.zrel_mount_t)) THEN
    46004608             itest(i)=1
    46014609             igwd=igwd+1
Note: See TracChangeset for help on using the changeset viewer.