Changeset 1330 for trunk/LMDZ.GENERIC


Ignore:
Timestamp:
Aug 26, 2014, 8:15:53 AM (10 years ago)
Author:
aslmd
Message:

LMDZ.GENERIC. a few corrections in edu version (within deftank).

Location:
trunk/LMDZ.GENERIC
Files:
1 added
6 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.GENERIC/deftank/edu/INIT/planet_start

    r1322 r1330  
    77g                    !! * 
    889.81                 !! <-- Acceleration gravite (m/s2)
     9mugaz                !! *
     1028.0                 !! <-- Masse molaire atmosphérique (g/mol)
     11rcp                  !! *
     120.286                !! <-- Rapport r/cp (massiques)
    913daysec               !! *
    101486400                !! <-- Duree du jour (s)
    … …  
    20240.                   !! <-- Date periastre (jours depuis printemps)
    2125obliquit             !! *
    22 25.0                 !! <-- Obliquite (degres)
     2623.9                 !! <-- Obliquite (degres)
    2327                 
    2428none             

  • trunk/LMDZ.GENERIC/deftank/edu/RUN/def_benj_earth/callphys.def

    r1322 r1330  
    3838corrk      = .true.
    3939# Include continuum absorption in radiative transfer (note CO2 is treated separately) ?
    40 continuum  = .true.
     40continuum  = .false.
    4141# folder in which correlated-k data is stored ?
    42  corrkdir    = earth
     42# corrkdir    = earth
    4343# corrkdir   = N2CO2poor_H2Ovar
    4444# corrkdir   = megaCO2
    45 # corrkdir   = null
     45 corrkdir   = null
    4646# call visible gaseous absorption in radiative transfer ?
    47 callgasvis = .true.
     47callgasvis = .false.
    4848# Include Rayleigh scattering in the visible ?
    49 rayleigh   = .true.
     49rayleigh   = .false.
    5050# Uniform absorption coefficient in radiative transfer?
    51 graybody   = .false.
     51graybody   = .true.
    5252# Constant absorption coefficient in visible
    5353#      (in m^2/kg; only if graybody=true):
    5454#      tau_surf= kappa*P/g
    55 kappa_VI = 5.e-7
     55#kappa_VI = 5.e-7
     56kappa_VI = 5.e-20
    5657# Constant absorption coefficient in IR
    5758#      (in m^2/kg; only if graybody=true):
    … …  
    7980# Force specific heat capacity and molecular mass values
    8081force_cpp      = .true.
     82force_cpp      = .false.
    8183# Specific heat capacity in J K-1 kg-1 [only used if force_cpp = T]
    8284cpp            = 1039.918
    … …  
    158160
    159161# To not call abort when temperature is outside boundaries:
    160 strictboundcorrk = .false.
     162strictboundcorrk = .true.
    161163
    162164

  • trunk/LMDZ.GENERIC/deftank/edu/RUN/etu.def

    r1323 r1330  
    11# Nombre de jours d'intégration               
    2 nday = 730
     2nday = 1825
     3#nday = 3650
     4#nday = 7300
     5#nday = 14600
    36
    47# Nombre de pas d'intégration dynamique par jour
    … …  
    1114iphysiq = 1
    1215
    13 ##########
    14 # mode simple "gris" (coefficient d'absorption constant selon lambda)
    15 graybody  = .true.
    16 continuum = .false.
    17 corrkdir  = null
    18 # coefficient d'absorption IR et VIS en m2/kg [--> kappa ~ tau_surf*g/P]
     16# Coefficient d'absorption IR en m2/kg [--> kappa ~ tau_surf*g/P]
    1917kappa_IR = 5.e-5
    20 kappa_VI = 5.e-6
    21 ##########
    2218
    2319

  • trunk/LMDZ.GENERIC/deftank/edu/atlas.exe

    r1323 r1330  
    99#########################################################################
    1010
     11# arguments
     12from optparse import OptionParser ### TBR by argparse
     13parser = OptionParser()
     14(opt,args) = parser.parse_args()
     15if len(args) == 0: args = "resultat.nc"
     16
    1117# OK. then the actual plotting commands.
    1218# -- people should be encouraged to use the command line pp.py too
    13 
    1419from ppclass import pp
    1520
    16 fi="resultat.nc"
     21fi=args
    1722
    18 t = pp(var="tsurf",file =fi,x="-180,180",y=[-80,-60,-20,0],superpose=True,quiet=True,ylabel=r'surface temperature ($^{\circ}$C)').get()-273.15
     23t = pp(var="tsurf",file=fi,x="-180,180",y=[-80,-60,-20,0],superpose=True,quiet=True,ylabel=r'surface temperature ($^{\circ}$C)').get()-273.15
    1924t.plot()
    2025

  • trunk/LMDZ.GENERIC/deftank/edu/run.exe

    r1322 r1330  
    1414b1=3 # recompile with full if you changed this
    1515b2=2 # recompile with full if you changed this
     16b1=1 # recompile with full if you changed this
     17b2=1 # recompile with full if you changed this
    1618
    17 echo "*** COMPILATION GCM ***"
     19echo "*** COMPILATION GCM *** ne pas interrompre ***"
    1820cd $wheresource/LMDZ.COMMON ; \rm gcm.e
    1921makelmdz -cpp NODYN -d $nx"x"$ny"x"$nz -b $b1"x"$b2 -t $tr -s 1 -p generic -arch gfortran gcm > logcompilegcm 2> logcompilegcm
    … …  
    2224fi
    2325
    24 echo "*** COMPILATION NEWSTART ***"
     26echo "*** COMPILATION NEWSTART *** ne pas interrompre ***"
    2527cd $wheresource/LMDZ.GENERIC ; \rm newstart.e
    2628makegcm_gfortran_local -d $nx"x"$ny"x"$nz -debug newstart > logcompilerestart 2> logcompilerestart
    … …  
    3941cd $thisfolder/RUN ; \rm resultat.nc > /dev/null 2> /dev/null
    4042time gcm.e | tee log_gcm | grep "Ls ="
     43#time gcm.e > /dev/null 2> /dev/null
     44#time gcm.e #| grep "WRITE"
    4145mv diagfi.nc ../resultat.nc
    4246

  • trunk/LMDZ.GENERIC/libf/phystd/newstart.F

    r1321 r1330  
    430430      mugaz = p_mugaz
    431431      daysec = p_daysec
     432
     433      if (p_omeg .eq. -9999.) then
     434        p_omeg = 8.*atan(1.)/p_daysec
     435        print*,"new value of omega",p_omeg
     436      endif
    432437
    433438      kappa = 8.314*1000./(p_mugaz*p_cpp) ! added by RDW
Note: See TracChangeset for help on using the changeset viewer.