Changeset 963


Ignore:
Timestamp:
Jun 9, 2008, 11:54:29 AM (17 years ago)
Author:
Laurent Fairhead
Message:

On passe les parametres lus dans parameter_mix.data et supcrit.data dans le
physiq.def
LF

Location:
LMDZ4/trunk/libf/phylmd
Files:
5 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/trunk/libf/phylmd/YOMCST2.h

    r879 r963  
    11
    2       INTEGER choice
     2      INTEGER choice, iflag_mix
    33      REAL  gammas, alphas, betas, Fmax, qqa1, qqa2, qqa3, scut
    44      REAL  Qcoef1max,Qcoef2max,Supcrit1,Supcrit2
    5 C
    6       COMMON/YOMCST2/gammas,    alphas, betas, Fmax, choice, scut,
    7      $               qqa1, qqa2, qqa3,
    8      $               Qcoef1max,Qcoef2max,
    9      $               Supcrit1, Supcrit2
    10 C    --------------------------------------------------------------------
     5!
     6      COMMON/YOMCST2/gammas,    alphas, betas, Fmax, scut,              &
     7     &               qqa1, qqa2, qqa3,                                  &
     8     &               Qcoef1max,Qcoef2max,                               &
     9     &               Supcrit1, Supcrit2,                                &
     10     &               choice,iflag_mix
     11!    --------------------------------------------------------------------
    1112

  • LMDZ4/trunk/libf/phylmd/concvl.F

    r940 r963  
    118118       REAL zx_t,zdelta,zx_qs,zcor
    119119c
    120        INTEGER iflag_mix
    121        SAVE iflag_mix
     120!       INTEGER iflag_mix
     121!       SAVE iflag_mix
    122122       INTEGER noff, minorig
    123123       INTEGER i,k,itra
    … …  
    215215C
    216216      if (iflag_con.eq.3) then
    217 c     CALL cv3_inicp(iflag_clos,iflag_mix)
    218       CALL cv3_inip(iflag_mix)
     217c     CALL cv3_inicp()
     218      CALL cv3_inip()
    219219      endif
    220220c
    … …  
    223223C===========================================================================
    224224C
    225          open (56,file='supcrit.data')
    226          read (56,*) Supcrit1, Supcrit2
    227          close (56)
     225c$$$         open (56,file='supcrit.data')
     226c$$$         read (56,*) Supcrit1, Supcrit2
     227c$$$         close (56)
    228228c
    229229         print*, 'supcrit1, supcrit2' ,supcrit1, supcrit2

  • LMDZ4/trunk/libf/phylmd/conf_phys.F90

    r956 r963  
    2222 include "nuage.h"
    2323 include "YOMCST.h"
     24 include "YOMCST2.h"
    2425!IM : on inclut/initialise les taux de CH4, N2O, CFC11 et CFC12
    2526include "clesphys.h"
    … …  
    8384  integer,SAVE :: iflag_coupl_omp,iflag_clos_omp,iflag_wake_omp
    8485  integer,SAVE :: iflag_cvl_sigd_omp
     86  REAL, SAVE :: supcrit1_omp, supcrit2_omp
     87  INTEGER, SAVE :: iflag_mix_omp
     88  real, save :: scut_omp, qqa1_omp, qqa2_omp, gammas_omp, Fmax_omp, alphas_omp
    8589
    8690  REAL,SAVE :: R_ecc_omp,R_peri_omp,R_incl_omp,solaire_omp,co2_ppm_omp
    … …  
    991995  ok_slab_sicOBS_omp = .true.
    992996  call getin('ok_slab_sicOBS', ok_slab_sicOBS_omp)
     997!
     998!
     999! PARAMETERS FOR CONVECTIVE INHIBITION BY TROPOS. DRYNESS
     1000!
     1001!Config Key  = supcrit1
     1002!Config Desc =
     1003!Config Def  = .540
     1004!Config Help =
     1005!
     1006  supcrit1_omp = .540
     1007  call getin('supcrit1',supcrit1_omp)
     1008
     1009!
     1010!Config Key  = supcrit2
     1011!Config Desc =
     1012!Config Def  = .600
     1013!Config Help =
     1014!
     1015  supcrit2_omp = .600
     1016  call getin('supcrit2',supcrit2_omp)
     1017
     1018!
     1019! PARAMETERS FOR THE MIXING DISTRIBUTION
     1020!
     1021!
     1022!Config Key  = iflag_mix
     1023!Config Desc =
     1024!Config Def  = 1
     1025!Config Help =
     1026!
     1027  iflag_mix_omp = 1
     1028  call getin('iflag_mix',iflag_mix_omp)
     1029
     1030!
     1031!Config Key  = scut
     1032!Config Desc =
     1033!Config Def  = 0.95
     1034!Config Help =
     1035!
     1036  scut_omp = 0.95
     1037  call getin('scut',scut_omp)
     1038
     1039!
     1040!Config Key  = qqa1
     1041!Config Desc =
     1042!Config Def  = 1.0
     1043!Config Help =
     1044!
     1045  qqa1_omp = 1.0
     1046  call getin('qqa1',qqa1_omp)
     1047
     1048!
     1049!Config Key  = qqa2
     1050!Config Desc =
     1051!Config Def  = 0.0
     1052!Config Help =
     1053!
     1054  qqa2_omp = 0.0
     1055  call getin('qqa2',qqa2_omp)
     1056
     1057!
     1058!Config Key  = gammas
     1059!Config Desc =
     1060!Config Def  = 0.05
     1061!Config Help =
     1062!
     1063  gammas_omp = 0.05
     1064  call getin('gammas',gammas_omp)
     1065
     1066!
     1067!Config Key  = Fmax
     1068!Config Desc =
     1069!Config Def  = 0.65
     1070!Config Help =
     1071!
     1072  Fmax_omp = 0.65
     1073  call getin('Fmax',Fmax_omp)
     1074
     1075!
     1076!Config Key  = alphas 
     1077!Config Desc =
     1078!Config Def  = -5.
     1079!Config Help =
     1080!
     1081  alphas_omp = -5.
     1082  call getin('alphas',alphas_omp)
     1083
    9931084
    9941085
    … …  
    10971188    cvl_corr = cvl_corr_omp
    10981189    ok_lic_melt = ok_lic_melt_omp
     1190    supcrit1 = supcrit1_omp
     1191    supcrit2 = supcrit2_omp
     1192    iflag_mix = iflag_mix_omp
     1193    scut = scut_omp
     1194    qqa1 = qqa1_omp
     1195    qqa2 = qqa2_omp
     1196    gammas = gammas_omp
     1197    Fmax = Fmax_omp
     1198    alphas = alphas_omp
    10991199
    11001200!$OMP MASTER
    … …  
    11781278  write(numout,*)' inertie_ice = ', inertie_ice
    11791279  write(numout,*)' inertie_sno = ', inertie_sno
     1280  write(numout,*)' supcrit1 = ', supcrit1
     1281  write(numout,*)' supcrit2 = ', supcrit2
     1282  write(numout,*)' iflag_mix = ', iflag_mix
     1283  write(numout,*)' scut = ', scut
     1284  write(numout,*)' qqa1 = ', qqa1
     1285  write(numout,*)' qqa2 = ', qqa2
     1286  write(numout,*)' gammas = ', gammas
     1287  write(numout,*)' Fmax = ', Fmax
     1288  write(numout,*)' alphas = ', alphas
    11801289
    11811290  write(numout,*)' lonmin lonmax latmin latmax bilKP_ins =',&

  • LMDZ4/trunk/libf/phylmd/cv3_inicp.F

    r879 r963  
    1         SUBROUTINE cv3_inicp(iflag_clos,iflag_mix)
     1        SUBROUTINE cv3_inicp()
    22*
    33***************************************************************

  • LMDZ4/trunk/libf/phylmd/cv3_inip.F

    r879 r963  
    1         SUBROUTINE cv3_inip(iflag_mix)
     1        SUBROUTINE cv3_inip()
    22***************************************************************
    33*                                                             *
    … …  
    1111#include "YOMCST2.h"
    1212c
    13       INTEGER iflag_mix
     13c      INTEGER iflag_mix
    1414c
    1515c --   Mixing probability distribution functions
    … …  
    5252C
    5353
    54         open(57,file='parameter_mix.data')
     54c$$$        open(57,file='parameter_mix.data')
     55c$$$
     56c$$$        read(57,*) iflag_mix, scut
     57c$$$        read(57,*)
     58c$$$        if(iflag_mix .gt. 0) then
     59c$$$          read(57,*) qqa1, qqa2
     60c$$$              read(57,*)
     61c$$$              read(57,*) gammas, Fmax
     62c$$$              read(57,*)
     63c$$$              read(57,*) alphas
     64c$$$         endif
     65c$$$     close(57)
    5566
    56         read(57,*) iflag_mix, scut
    57         read(57,*)
    58         if(iflag_mix .gt. 0) then
    59               read(57,*) qqa1, qqa2
    60               read(57,*)
    61               read(57,*) gammas, Fmax
    62               read(57,*)
    63               read(57,*) alphas
    64          endif
    65          close(57)
    6667c
    6768      if(iflag_mix .gt. 0) then
Note: See TracChangeset for help on using the changeset viewer.