Changeset 1059


Ignore:
Timestamp:
2008-06-05T11:47:29+02:00 (12 years ago)
Author:
ctlod
Message:

trunk: follow DOCTOR rules for TKE namelist parameters names and zdftke.F90 style code review, see ticket: #188

Location:
trunk/CONFIG
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • trunk/CONFIG/GYRE/EXP00/namelist

    r1052 r1059  
    200200!                =1 observed ice-cover          , 
    201201!                =2 ice-model used                                      ("key_lim3" or "key_lim2) 
     202!  nn_ico_cpl    ice-ocean coupling : =0 LIM-3 old case  
     203!                                     =1 stresses computed using now ocean velocity 
     204!                                     =2 combination of 0 and 1 cases 
    202205!  ln_dm2dc     Daily mean to Diurnal Cycle short wave (qsr) 
    203206!  ln_ssr       Sea Surface Restoring on T and/or S (fill namsbc_ssr) 
     
    213216   ln_cpl      = .false. 
    214217   nn_ice      = 0 
     218   nn_ico_cpl  = 0 
    215219   ln_dm2dc    = .false. 
    216220   ln_rnf      = .false. 
     
    709713!       namtke   turbulent eddy kinetic dependent vertical diffusion    ("key_zdftke") 
    710714!----------------------------------------------------------------------- 
    711 !  ln_rstke flag to restart with tke from a run without tke (default F) 
    712 !  ediff    coef. to compute vertical eddy coef. (avt=ediff*mxl*sqrt(e) ) 
    713 !  ediss    coef. of the Kolmogoroff dissipation 
    714 !  ebb      coef. of the surface input of tke 
    715 !  efave    coef. to applied to the tke diffusion ( avtke=efave*avm ) 
    716 !  emin     minimum value of tke (m^2/s^2) 
    717 !  emin0    surface minimum value of tke (m^2/s^2) 
    718 !  nitke    number of restart iterative loops 
    719 !  ri_c     critic richardson number = 2/9 = 0.22222222 (hard coded) 
    720 !  nmxl     length used   = 0 bounded by the distance to surface and bottom 
    721 !           !             = 1 bounded by the local vertical scale factor 
    722 !           !             = 2 first vertical derivative of mixing length bounded by 1 
    723 !  npdl     Prandtl number   = 0 no vertical prandtl number (avt=avm) 
    724 !           !                = 1 prandtl number function of richarson number (avt=pdl*avm) 
    725 !           !                = 2 same as = 1 but a shapiro filter is applied on pdl 
    726 !  nave     =  horizontal averaged (=1) or not (=0) of avt  (default =1) 
    727 !  navb     = 0 cst background avt0, avm0 / =1 profile used on avtb 
    728 !  ln_lsfc  computation surface value of buoyancy lenght scale as function of wind stress 
    729 !  rn_lmin  buoyancy lenght scale minimum value if ln_lsfc = F 
    730 !  nn_etau  to add tke induced by wind  = 0 no add tke induced by wind 
    731 !           !                           = 1 add tke induced by wind 
    732 !           !                           = 2 add tke induced by wind only at the base of the mixed layer 
    733 !  nn_htau  TKE penetration profile type 
    734 !  rn_efr   fraction of TKE surface value which penetrates inside the thermocline 
    735 !  ln_lc    flag to take into account Langmuir Circulation (LC) 
    736 !  rn_lc    coef to compute vertical velocity of LC 
    737 !  nn_avtb_2d equatorial horizontal shape (=0 uniform) (=1 linear decrease and constant) 
     715!  ln_rstke    flag to restart with tke from a run without tke (default F) 
     716!  rn_ediff    coef. to compute vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) ) 
     717!  rn_ediss    coef. of the Kolmogoroff dissipation 
     718!  rn_ebb      coef. of the surface input of tke 
     719!  rn_efave    coef. to applied to the tke diffusion ( avtke=rn_efave*avm ) 
     720!  rn_emin     minimum value of tke (m^2/s^2) 
     721!  rn_emin0    surface minimum value of tke (m^2/s^2) 
     722!  nn_itke     number of restart iterative loops 
     723!  rn_cri      critic richardson number = 2/9 = 0.22222222 (hard coded) 
     724!  nn_mxl      length used = 0 bounded by the distance to surface and bottom 
     725!              !           = 1 bounded by the local vertical scale factor 
     726!              !           = 2 first vertical derivative of mixing length bounded by 1 
     727!              !           = 3 same criteria as case 2 but applied in a different way 
     728!  nn_pdl      Prandtl number = 0 no vertical prandtl number (avt=avm) 
     729!              !              = 1 prandtl number function of richarson number (avt=nn_pdl*avm) 
     730!  nn_ave      = horizontal averaged (=1) or not (=0) of avt  (default =1) 
     731!  nn_avb      = 0 cst background avt0, avm0 / =1 profile used on avtb 
     732!  ln_mxl0     mixing length scale surface value as function of wind stress or not 
     733!  rn_lmin     interior buoyancy lenght scale minimum value  
     734!  rn_lmin0    surface  buoyancy lenght scale minimum value  
     735!  nn_etau     to add tke induced by wind  = 0 no add tke induced by wind 
     736!              !                           = 1 add tke induced by wind 
     737!              !                           = 2 add tke induced by wind only at the mxl base 
     738!  nn_htau     TKE penetration profile type 
     739!  rn_efr      fraction of TKE surface value which penetrates inside the thermocline 
     740!  ln_lc       flag to take into account Langmuir Circulation (LC) 
     741!  rn_lc       coef to compute vertical velocity of LC 
     742!  nn_havtb    horizontal shape or not for avtb (=0/1) 
    738743&namtke 
    739744   ln_rstke = .false. 
    740    ediff =       0.1 
    741    ediss =       0.7 
    742    ebb   =      3.75 
    743    efave =        1. 
    744    emin  =     1.e-5 
    745    emin0 =     1.e-4 
    746    nitke =        50 
    747    nmxl  =         2 
    748    npdl  =         1 
    749    navb  =         0 
    750    nave  =         1 
    751    ln_lsfc =  .false. 
    752    rn_lmin =     0.4 
    753    nn_etau =       0 
    754    nn_htau =       2 
    755    rn_efr  =    0.05 
    756    ln_lc   =  .false. 
    757    rn_lc   =    0.15 
    758    nn_avtb_2d =    0 
     745   rn_ediff =     0.1 
     746   rn_ediss =     0.7 
     747   rn_ebb   =    3.75 
     748   rn_efave =      1. 
     749   rn_emin  =   1.e-5 
     750   rn_emin0 =   1.e-4 
     751   nn_itke  =      50 
     752   nn_mxl   =       2 
     753   nn_pdl   =       1 
     754   nn_avb   =       0 
     755   nn_ave   =       1 
     756   ln_mxl0  = .false. 
     757   rn_lmin  =     0.4 
     758   rn_lmin0 =     0.4 
     759   nn_etau  =       0 
     760   nn_htau  =       2 
     761   rn_efr   =    0.05 
     762   ln_lc    = .false. 
     763   rn_lc    =    0.15 
     764   nn_havtb =       0 
    759765/ 
    760766!----------------------------------------------------------------------- 
  • trunk/CONFIG/GYRE_LOBSTER/EXP00/namelist

    r1052 r1059  
    200200!                =1 observed ice-cover          , 
    201201!                =2 ice-model used                                      ("key_lim3" or "key_lim2) 
     202!  nn_ico_cpl    ice-ocean coupling : =0 LIM-3 old case  
     203!                                     =1 stresses computed using now ocean velocity 
     204!                                     =2 combination of 0 and 1 cases 
    202205!  ln_dm2dc     Daily mean to Diurnal Cycle short wave (qsr) 
    203206!  ln_ssr       Sea Surface Restoring on T and/or S (fill namsbc_ssr) 
     
    213216   ln_cpl      = .false. 
    214217   nn_ice      = 0 
     218   nn_ico_cpl  = 0 
    215219   ln_dm2dc    = .false. 
    216220   ln_rnf      = .false. 
     
    709713!       namtke   turbulent eddy kinetic dependent vertical diffusion    ("key_zdftke") 
    710714!----------------------------------------------------------------------- 
    711 !  ln_rstke flag to restart with tke from a run without tke (default F) 
    712 !  ediff    coef. to compute vertical eddy coef. (avt=ediff*mxl*sqrt(e) ) 
    713 !  ediss    coef. of the Kolmogoroff dissipation 
    714 !  ebb      coef. of the surface input of tke 
    715 !  efave    coef. to applied to the tke diffusion ( avtke=efave*avm ) 
    716 !  emin     minimum value of tke (m^2/s^2) 
    717 !  emin0    surface minimum value of tke (m^2/s^2) 
    718 !  nitke    number of restart iterative loops 
    719 !  ri_c     critic richardson number = 2/9 = 0.22222222 (hard coded) 
    720 !  nmxl     length used   = 0 bounded by the distance to surface and bottom 
    721 !           !             = 1 bounded by the local vertical scale factor 
    722 !           !             = 2 first vertical derivative of mixing length bounded by 1 
    723 !  npdl     Prandtl number   = 0 no vertical prandtl number (avt=avm) 
    724 !           !                = 1 prandtl number function of richarson number (avt=pdl*avm) 
    725 !           !                = 2 same as = 1 but a shapiro filter is applied on pdl 
    726 !  nave     =  horizontal averaged (=1) or not (=0) of avt  (default =1) 
    727 !  navb     = 0 cst background avt0, avm0 / =1 profile used on avtb 
    728 !  ln_lsfc  computation surface value of buoyancy lenght scale as function of wind stress 
    729 !  rn_lmin  buoyancy lenght scale minimum value if ln_lsfc = F 
    730 !  nn_etau  to add tke induced by wind  = 0 no add tke induced by wind 
    731 !           !                           = 1 add tke induced by wind 
    732 !           !                           = 2 add tke induced by wind only at the base of the mixed layer 
    733 !  nn_htau  TKE penetration profile type 
    734 !  rn_efr   fraction of TKE surface value which penetrates inside the thermocline 
    735 !  ln_lc    flag to take into account Langmuir Circulation (LC) 
    736 !  rn_lc    coef to compute vertical velocity of LC 
    737 !  nn_avtb_2d equatorial horizontal shape (=0 uniform) (=1 linear decrease and constant) 
     715!  ln_rstke    flag to restart with tke from a run without tke (default F) 
     716!  rn_ediff    coef. to compute vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) ) 
     717!  rn_ediss    coef. of the Kolmogoroff dissipation 
     718!  rn_ebb      coef. of the surface input of tke 
     719!  rn_efave    coef. to applied to the tke diffusion ( avtke=rn_efave*avm ) 
     720!  rn_emin     minimum value of tke (m^2/s^2) 
     721!  rn_emin0    surface minimum value of tke (m^2/s^2) 
     722!  nn_itke     number of restart iterative loops 
     723!  rn_cri      critic richardson number = 2/9 = 0.22222222 (hard coded) 
     724!  nn_mxl      length used = 0 bounded by the distance to surface and bottom 
     725!              !           = 1 bounded by the local vertical scale factor 
     726!              !           = 2 first vertical derivative of mixing length bounded by 1 
     727!              !           = 3 same criteria as case 2 but applied in a different way 
     728!  nn_pdl      Prandtl number = 0 no vertical prandtl number (avt=avm) 
     729!              !              = 1 prandtl number function of richarson number (avt=nn_pdl*avm) 
     730!  nn_ave      = horizontal averaged (=1) or not (=0) of avt  (default =1) 
     731!  nn_avb      = 0 cst background avt0, avm0 / =1 profile used on avtb 
     732!  ln_mxl0     mixing length scale surface value as function of wind stress or not 
     733!  rn_lmin     interior buoyancy lenght scale minimum value  
     734!  rn_lmin0    surface  buoyancy lenght scale minimum value  
     735!  nn_etau     to add tke induced by wind  = 0 no add tke induced by wind 
     736!              !                           = 1 add tke induced by wind 
     737!              !                           = 2 add tke induced by wind only at the mxl base 
     738!  nn_htau     TKE penetration profile type 
     739!  rn_efr      fraction of TKE surface value which penetrates inside the thermocline 
     740!  ln_lc       flag to take into account Langmuir Circulation (LC) 
     741!  rn_lc       coef to compute vertical velocity of LC 
     742!  nn_havtb    horizontal shape or not for avtb (=0/1) 
    738743&namtke 
    739744   ln_rstke = .false. 
    740    ediff =       0.1 
    741    ediss =       0.7 
    742    ebb   =      3.75 
    743    efave =        1. 
    744    emin  =     1.e-5 
    745    emin0 =     1.e-4 
    746    nitke =        50 
    747    nmxl  =         2 
    748    npdl  =         1 
    749    navb  =         0 
    750    nave  =         1 
    751    ln_lsfc =  .false. 
    752    rn_lmin =     0.4 
    753    nn_etau =       0 
    754    nn_htau =       2 
    755    rn_efr  =    0.05 
    756    ln_lc   =  .false. 
    757    rn_lc   =    0.15 
    758    nn_avtb_2d =    0 
     745   rn_ediff =     0.1 
     746   rn_ediss =     0.7 
     747   rn_ebb   =    3.75 
     748   rn_efave =      1. 
     749   rn_emin  =   1.e-5 
     750   rn_emin0 =   1.e-4 
     751   nn_itke  =      50 
     752   nn_mxl   =       2 
     753   nn_pdl   =       1 
     754   nn_avb   =       0 
     755   nn_ave   =       1 
     756   ln_mxl0  = .false. 
     757   rn_lmin  =     0.4 
     758   rn_lmin0 =     0.4 
     759   nn_etau  =       0 
     760   nn_htau  =       2 
     761   rn_efr   =    0.05 
     762   ln_lc    = .false. 
     763   rn_lc    =    0.15 
     764   nn_havtb =       0 
    759765/ 
    760766!----------------------------------------------------------------------- 
  • trunk/CONFIG/ORCA2_LIM/EXP00/namelist

    r1055 r1059  
    713713!       namtke   turbulent eddy kinetic dependent vertical diffusion    ("key_zdftke") 
    714714!----------------------------------------------------------------------- 
    715 !  ln_rstke flag to restart with tke from a run without tke (default F) 
    716 !  ediff    coef. to compute vertical eddy coef. (avt=ediff*mxl*sqrt(e) ) 
    717 !  ediss    coef. of the Kolmogoroff dissipation 
    718 !  ebb      coef. of the surface input of tke 
    719 !  efave    coef. to applied to the tke diffusion ( avtke=efave*avm ) 
    720 !  emin     minimum value of tke (m^2/s^2) 
    721 !  emin0    surface minimum value of tke (m^2/s^2) 
    722 !  nitke    number of restart iterative loops 
    723 !  ri_c     critic richardson number = 2/9 = 0.22222222 (hard coded) 
    724 !  nmxl     length used   = 0 bounded by the distance to surface and bottom 
    725 !           !             = 1 bounded by the local vertical scale factor 
    726 !           !             = 2 first vertical derivative of mixing length bounded by 1 
    727 !  npdl     Prandtl number   = 0 no vertical prandtl number (avt=avm) 
    728 !           !                = 1 prandtl number function of richarson number (avt=pdl*avm) 
    729 !           !                = 2 same as = 1 but a shapiro filter is applied on pdl 
    730 !  nave     =  horizontal averaged (=1) or not (=0) of avt  (default =1) 
    731 !  navb     = 0 cst background avt0, avm0 / =1 profile used on avtb 
    732 !  ln_lsfc  computation surface value of buoyancy lenght scale as function of wind stress 
    733 !  rn_lmin  buoyancy lenght scale minimum value if ln_lsfc = F 
    734 !  nn_etau  to add tke induced by wind  = 0 no add tke induced by wind 
    735 !           !                           = 1 add tke induced by wind 
    736 !           !                           = 2 add tke induced by wind only at the base of the mixed layer 
    737 !  nn_htau  TKE penetration profile type 
    738 !  rn_efr   fraction of TKE surface value which penetrates inside the thermocline 
    739 !  ln_lc    flag to take into account Langmuir Circulation (LC) 
    740 !  rn_lc    coef to compute vertical velocity of LC 
    741 !  nn_avtb_2d equatorial horizontal shape (=0 uniform) (=1 linear decrease and constant) 
     715!  ln_rstke    flag to restart with tke from a run without tke (default F) 
     716!  rn_ediff    coef. to compute vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) ) 
     717!  rn_ediss    coef. of the Kolmogoroff dissipation 
     718!  rn_ebb      coef. of the surface input of tke 
     719!  rn_efave    coef. to applied to the tke diffusion ( avtke=rn_efave*avm ) 
     720!  rn_emin     minimum value of tke (m^2/s^2) 
     721!  rn_emin0    surface minimum value of tke (m^2/s^2) 
     722!  nn_itke     number of restart iterative loops 
     723!  rn_cri      critic richardson number = 2/9 = 0.22222222 (hard coded) 
     724!  nn_mxl      length used = 0 bounded by the distance to surface and bottom 
     725!              !           = 1 bounded by the local vertical scale factor 
     726!              !           = 2 first vertical derivative of mixing length bounded by 1 
     727!              !           = 3 same criteria as case 2 but applied in a different way 
     728!  nn_pdl      Prandtl number = 0 no vertical prandtl number (avt=avm) 
     729!              !              = 1 prandtl number function of richarson number (avt=nn_pdl*avm) 
     730!  nn_ave      = horizontal averaged (=1) or not (=0) of avt  (default =1) 
     731!  nn_avb      = 0 cst background avt0, avm0 / =1 profile used on avtb 
     732!  ln_mxl0     mixing length scale surface value as function of wind stress or not 
     733!  rn_lmin     interior buoyancy lenght scale minimum value  
     734!  rn_lmin0    surface  buoyancy lenght scale minimum value  
     735!  nn_etau     to add tke induced by wind  = 0 no add tke induced by wind 
     736!              !                           = 1 add tke induced by wind 
     737!              !                           = 2 add tke induced by wind only at the mxl base 
     738!  nn_htau     TKE penetration profile type 
     739!  rn_efr      fraction of TKE surface value which penetrates inside the thermocline 
     740!  ln_lc       flag to take into account Langmuir Circulation (LC) 
     741!  rn_lc       coef to compute vertical velocity of LC 
     742!  nn_havtb    horizontal shape or not for avtb (=0/1) 
    742743&namtke 
    743744   ln_rstke = .false. 
    744    ediff =       0.1 
    745    ediss =       0.7 
    746    ebb   =      3.75 
    747    efave =        1. 
    748    emin  =     1.e-6 
    749    emin0 =     1.e-4 
    750    nitke =        50 
    751    nmxl  =         2 
    752    npdl  =         1 
    753    navb  =         0 
    754    nave  =         1 
    755    ln_lsfc =  .false. 
    756    rn_lmin =     0.4 
    757    nn_etau =       0 
    758    nn_htau =       2 
    759    rn_efr  =    0.05 
    760    ln_lc   =  .false. 
    761    rn_lc   =    0.15 
    762    nn_avtb_2d =    0 
     745   rn_ediff =     0.1 
     746   rn_ediss =     0.7 
     747   rn_ebb   =    3.75 
     748   rn_efave =      1. 
     749   rn_emin  =   1.e-6 
     750   rn_emin0 =   1.e-4 
     751   nn_itke  =      50 
     752   nn_mxl   =       2 
     753   nn_pdl   =       1 
     754   nn_avb   =       0 
     755   nn_ave   =       1 
     756   ln_mxl0  = .false. 
     757   rn_lmin  =     0.4 
     758   rn_lmin0 =     0.4 
     759   nn_etau  =       0 
     760   nn_htau  =       2 
     761   rn_efr   =    0.05 
     762   ln_lc    = .false. 
     763   rn_lc    =    0.15 
     764   nn_havtb =       0 
    763765/ 
    764766!----------------------------------------------------------------------- 
  • trunk/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist

    r1055 r1059  
    713713!       namtke   turbulent eddy kinetic dependent vertical diffusion    ("key_zdftke") 
    714714!----------------------------------------------------------------------- 
    715 !  ln_rstke flag to restart with tke from a run without tke (default F) 
    716 !  ediff    coef. to compute vertical eddy coef. (avt=ediff*mxl*sqrt(e) ) 
    717 !  ediss    coef. of the Kolmogoroff dissipation 
    718 !  ebb      coef. of the surface input of tke 
    719 !  efave    coef. to applied to the tke diffusion ( avtke=efave*avm ) 
    720 !  emin     minimum value of tke (m^2/s^2) 
    721 !  emin0    surface minimum value of tke (m^2/s^2) 
    722 !  nitke    number of restart iterative loops 
    723 !  ri_c     critic richardson number = 2/9 = 0.22222222 (hard coded) 
    724 !  nmxl     length used   = 0 bounded by the distance to surface and bottom 
    725 !           !             = 1 bounded by the local vertical scale factor 
    726 !           !             = 2 first vertical derivative of mixing length bounded by 1 
    727 !  npdl     Prandtl number   = 0 no vertical prandtl number (avt=avm) 
    728 !           !                = 1 prandtl number function of richarson number (avt=pdl*avm) 
    729 !           !                = 2 same as = 1 but a shapiro filter is applied on pdl 
    730 !  nave     =  horizontal averaged (=1) or not (=0) of avt  (default =1) 
    731 !  navb     = 0 cst background avt0, avm0 / =1 profile used on avtb 
    732 !  ln_lsfc  computation surface value of buoyancy lenght scale as function of wind stress 
    733 !  rn_lmin  buoyancy lenght scale minimum value if ln_lsfc = F 
    734 !  nn_etau  to add tke induced by wind  = 0 no add tke induced by wind 
    735 !           !                           = 1 add tke induced by wind 
    736 !           !                           = 2 add tke induced by wind only at the base of the mixed layer 
    737 !  nn_htau  TKE penetration profile type 
    738 !  rn_efr   fraction of TKE surface value which penetrates inside the thermocline 
    739 !  ln_lc    flag to take into account Langmuir Circulation (LC) 
    740 !  rn_lc    coef to compute vertical velocity of LC 
    741 !  nn_avtb_2d equatorial horizontal shape (=0 uniform) (=1 linear decrease and constant) 
     715!  ln_rstke    flag to restart with tke from a run without tke (default F) 
     716!  rn_ediff    coef. to compute vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) ) 
     717!  rn_ediss    coef. of the Kolmogoroff dissipation 
     718!  rn_ebb      coef. of the surface input of tke 
     719!  rn_efave    coef. to applied to the tke diffusion ( avtke=rn_efave*avm ) 
     720!  rn_emin     minimum value of tke (m^2/s^2) 
     721!  rn_emin0    surface minimum value of tke (m^2/s^2) 
     722!  nn_itke     number of restart iterative loops 
     723!  rn_cri      critic richardson number = 2/9 = 0.22222222 (hard coded) 
     724!  nn_mxl      length used = 0 bounded by the distance to surface and bottom 
     725!              !           = 1 bounded by the local vertical scale factor 
     726!              !           = 2 first vertical derivative of mixing length bounded by 1 
     727!              !           = 3 same criteria as case 2 but applied in a different way 
     728!  nn_pdl      Prandtl number = 0 no vertical prandtl number (avt=avm) 
     729!              !              = 1 prandtl number function of richarson number (avt=nn_pdl*avm) 
     730!  nn_ave      = horizontal averaged (=1) or not (=0) of avt  (default =1) 
     731!  nn_avb      = 0 cst background avt0, avm0 / =1 profile used on avtb 
     732!  ln_mxl0     mixing length scale surface value as function of wind stress or not 
     733!  rn_lmin     interior buoyancy lenght scale minimum value  
     734!  rn_lmin0    surface  buoyancy lenght scale minimum value  
     735!  nn_etau     to add tke induced by wind  = 0 no add tke induced by wind 
     736!              !                           = 1 add tke induced by wind 
     737!              !                           = 2 add tke induced by wind only at the mxl base 
     738!  nn_htau     TKE penetration profile type 
     739!  rn_efr      fraction of TKE surface value which penetrates inside the thermocline 
     740!  ln_lc       flag to take into account Langmuir Circulation (LC) 
     741!  rn_lc       coef to compute vertical velocity of LC 
     742!  nn_havtb    horizontal shape or not for avtb (=0/1) 
    742743&namtke 
    743744   ln_rstke = .false. 
    744    ediff =       0.1 
    745    ediss =       0.7 
    746    ebb   =      3.75 
    747    efave =        1. 
    748    emin  =     1.e-6 
    749    emin0 =     1.e-4 
    750    nitke =        50 
    751    nmxl  =         2 
    752    npdl  =         1 
    753    navb  =         0 
    754    nave  =         1 
    755    ln_lsfc =  .false. 
    756    rn_lmin =     0.4 
    757    nn_etau =       0 
    758    nn_htau =       2 
    759    rn_efr  =    0.05 
    760    ln_lc   =  .false. 
    761    rn_lc   =    0.15 
    762    nn_avtb_2d =    0 
     745   rn_ediff =     0.1 
     746   rn_ediss =     0.7 
     747   rn_ebb   =    3.75 
     748   rn_efave =      1. 
     749   rn_emin  =   1.e-6 
     750   rn_emin0 =   1.e-4 
     751   nn_itke  =      50 
     752   nn_mxl   =       2 
     753   nn_pdl   =       1 
     754   nn_avb   =       0 
     755   nn_ave   =       1 
     756   ln_mxl0  = .false. 
     757   rn_lmin  =     0.4 
     758   rn_lmin0 =     0.4 
     759   nn_etau  =       0 
     760   nn_htau  =       2 
     761   rn_efr   =    0.05 
     762   ln_lc    = .false. 
     763   rn_lc    =    0.15 
     764   nn_havtb =       0 
    763765/ 
    764766!----------------------------------------------------------------------- 
Note: See TracChangeset for help on using the changeset viewer.