Ignore:
Timestamp:
Aug 2, 2011, 11:13:07 AM (14 years ago)
Author:
emillour
Message:

Generic GCM

  • Massive update to version 0.7

EM+RW

Location:
trunk/LMDZ.GENERIC/deftank
Files:
4 added
4 deleted
10 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.GENERIC/deftank/autospinup

    r135 r253  
    1 #!/bin/bash
     1# @ job_name=gcm_3D
     2# @ output = $(job_name).$(jobid)
     3# @ error = $(output)
     4# @ job_type = serial
     5# @ class = D
     6# @ queue
    27
    3 #DEF_DIR=/san/home/rdword/gcm/spinups
    4 DEF_DIR=/u/rwlmd/spinups/
     8#grep "^# @" $0
     9#set -vx
     10#
    511
    6 #makegcm -d 32x24x20 -b 32x36 -t 2 -p std gcm
    7 #makegcm -d 32x24x20 -b 32x36 -t 2 -p std newstart
    8 
    9 newstart.e < start_earth
     12cp -f run_spinup.def run.def
     13newstart.e < start_earlymars
    1014startrestart
    1115
    12 echo "Starting short GCM spinup run #1..."
    13 cp -f $DEF_DIR/run1.def run.def
    14 gcm.e > lruns1
    15 startrestart
    16 
    17 echo "Starting short GCM spinup run #2..."
    18 cp -f $DEF_DIR/run2.def run.def
    19 gcm.e > lruns2
    20 startrestart
    21 
    22 echo "Finishing GCM spinup year..."
    23 cp -f $DEF_DIR/run_1styr_earth.def run.def
    24 gcm.e > lrun0
     16echo "Starting GCM spinup year..."
     17gcm.e
    2518startrestart0
    2619
    27 exit 0
     20echo "Running for multiple years to reach equilibrium."
     21cp -f run_normal.def run.def
     22run0
    2823
    29 # run for many years to reach equilibrium
    30 #cp -f run_1yr_earth.def run.def
    31 #rm -f num_run
    32 #run0
    33 
    34 # start the final 1-yr run with intensive OLR save
    35 #cp -f run_1yr_earth_OLR.def run.def
    36 #run_OLR

  • trunk/LMDZ.GENERIC/deftank/callphys.earlymars

    r135 r253  
    44tracer    = .true.
    55# Diurnal cycle ?  if diurnal=false, diurnally averaged solar heating
    6 diurnal   = .false.
     6diurnal   = .true.
    77# Seasonal cycle ? if season=false, Ls stays constant, to value set in "start"
    88season    = .true.
    … …  
    1414lwrite    = .false.
    1515# Save statistics in file "stats.nc" ?
    16 callstats = .false.
     16callstats = .true.
    1717# Test energy conservation of model physics ?
    1818enertest  = .false.
    … …  
    3535# This is used only in the aerosol radiative transfer setup. (see aerave.F)
    3636tplanet    = 215.
    37 # Output spectral OLR in 3D?
     37# Output spectral OLR in 1D/3D?
    3838specOLR    = .false.
    39 # Output global radiative balance in file 'olr.txt' - slow for 1D!!
    40 meanOLR    = .false.
     39# Output global radiative balance in file 'rad_bal.out' - slow for 1D!!
     40meanOLR    = .true.
    4141# Variable gas species: Radiatively active ?
    42 varactive  = .false.
     42varactive  = .true.
    4343# Variable gas species: Fixed vertical distribution ?
    44 varfixed   = .true.
     44varfixed   = .false.
    4545# Variable gas species: Saturation percentage value at ground ?
    46 satval     = 0.8
     46satval     = 0.0
    4747
    4848## Star type
    … …  
    5858# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    5959# Stellar flux at 1 AU. Examples:
    60 # 1356.0 W m-2          Sol today
    61 # 1017.0 W m-2          Sol today x 0.75 = weak early Sun
    62 # 18.462 W m-2          The feeble Gl581
    63 Fat1AU = 1017.0
     60# 1366.0 W m-2          Sol today
     61# 1024.5 W m-2          Sol today x 0.75 = weak Sun 3.8 GYa
     62# 18.462 W m-2          The feeble GJ581
     63# 19.960 W m-2          GJ581 with e=0.38 orbital average
     64Fat1AU = 1024.5
    6465
    6566## Tracer and aerosol options
    … …  
    7879
    7980#########################################################################
    80 ## extra non-standard definitions for Earth / Early Mars
     81## extra non-standard definitions for Gliese 581d
    8182#########################################################################
    82  
    83 ## Radiative transfer options
    84 ## ~~~~~~~~~~~~~~~~~~~~~~~~~~
    85 # Parameterized Earth-like ozone absorption ?
    86 ozone         = .false.
    8783
    8884## Tracer and aerosol options
    8985## ~~~~~~~~~~~~~~~~~~~~~~~~~~
    9086# Fixed aerosol distributions?
    91 aerofixed     = .true.
     87aerofixed     = .false.
     88# Varying H2O cloud fraction?
     89CLFvarying    = .false.
     90# H2O cloud fraction?
     91CLFfixval     = 0.5
    9292# number mixing ratio of CO2 ice particles
    9393Nmix_co2      = 100000.
    … …  
    105105# WATER: Precipitation threshold (simple scheme only) ?
    106106rainthreshold = 0.011
     107# WATER: hydrology ?
     108hydrology     = .true.
    107109
    108110## CO2 options
    … …  
    114116# Set initial temperature profile to 1 K above CO2 condensation everywhere?
    115117nearco2cond   = .false.
     118
     119## N2 options
     120## ~~~~~~~~~~
     121# include N2 in cp and Rayleigh scattering calculation ?
     122addn2 = .false.
     123# N2 mixing ratio
     124n2mixratio = 0.0
     125
     126## H2 options
     127## ~~~~~~~~~~
     128# include H2 in cp and Rayleigh scattering calculation ?
     129# for now we assume H2 dominant if included
     130addh2 = .false.

  • trunk/LMDZ.GENERIC/deftank/callphys.earth

    r135 r253  
    2323callrad    = .true.
    2424# the rad. transfer is computed every "iradia" physical timestep
    25 iradia     = 4
     25iradia     = 6
    2626# call multilayer correlated-k radiative transfer ?
    2727corrk      = .true.
    2828# folder in which correlated-k data is stored ?
    29 corrkdir   = earth
     29corrkdir   = N2CO2poor_H2Ovar
    3030# call visible gaseous absorption in radiative transfer ?
    3131callgasvis = .true.
    … …  
    3737# Output spectral OLR in 3D?
    3838specOLR    = .false.
    39 # Output global radiative balance in file 'olr.txt' - slow for 1D!!
    40 meanOLR    = .false.
     39# Output global radiative balance in file 'rad_bal.out' - slow for 1D!!
     40meanOLR    = .true.
    4141# Variable gas species: Radiatively active ?
    42 varactive  = .false.
     42varactive  = .true.
    4343# Variable gas species: Fixed vertical distribution ?
    44 varfixed   = .true.
     44varfixed   = .false.
    4545# Variable gas species: Saturation percentage value at ground ?
    46 satval     = 0.8
     46satval     = 0.0
    4747
    4848## Star type
    … …  
    5858# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    5959# Stellar flux at 1 AU. Examples:
    60 # 1356.0 W m-2          Sol today
    61 # 1017.0 W m-2          Sol today x 0.75 = weak early Sun
     60# 1366.0 W m-2          Sol today
     61# 1024.5 W m-2          Sol today x 0.75 = weak early Sun
    6262# 18.462 W m-2          The feeble Gl581
    63 Fat1AU = 1356.0
     63Fat1AU = 1366.0
    6464
    6565## Tracer and aerosol options
    … …  
    7878
    7979#########################################################################
    80 ## extra non-standard definitions for Earth / Early Mars
     80## extra non-standard definitions for Earth
    8181#########################################################################
    8282 
    83 ## Radiative transfer options
    84 ## ~~~~~~~~~~~~~~~~~~~~~~~~~~
    85 # Parameterized Earth-like ozone absorption ?
    86 ozone         = .false.
    87 
    8883## Tracer and aerosol options
    8984## ~~~~~~~~~~~~~~~~~~~~~~~~~~
    9085# Fixed aerosol distributions?
    91 aerofixed     = .true.
     86aerofixed     = .false.
     87# Varying H2O cloud fraction?
     88CLFvarying    = .true.
     89# H2O cloud fraction if fixed?
     90CLFfixval     = 0.5
    9291# number mixing ratio of CO2 ice particles
    9392Nmix_co2      = 100000.
    … …  
    105104# WATER: Precipitation threshold (simple scheme only) ?
    106105rainthreshold = 0.011
     106# WATER: hydrology ?
     107hydrology     = .true.
    107108
    108109## CO2 options
    … …  
    114115# Set initial temperature profile to 1 K above CO2 condensation everywhere?
    115116nearco2cond   = .false.
     117
     118## N2 options
     119## ~~~~~~~~~~
     120# include N2 in cp and Rayleigh scattering calculation ?
     121addn2 = .false.
     122# N2 mixing ratio
     123n2mixratio = 0.0
     124
     125## H2 options
     126## ~~~~~~~~~~
     127# include H2 in cp and Rayleigh scattering calculation ?
     128# for now we assume H2 dominant if included
     129addh2 = .false.

  • trunk/LMDZ.GENERIC/deftank/run0

    r135 r253  
    44# SET HERE the maximum total number of simulations
    55
    6 set nummax=2
     6set nummax=5
    77
    88###########################################################################
    … …  
    4949
    5050# Run GCM
    51 nice gcm.e >! lrun${numnew}
     51gcm.e >! lrun${numnew}
    5252
    5353
    … …  
    6666  endif
    6767 echo 99999 >! num_run
    68 ############## To receive an Email message if the run crashes ########
    69 mail -s "crash run GCM" $address <<ENDMAIL
    70 The run on $machine in $dir has just crashed.
    71 ENDMAIL
     68############## To receive an email message if the run crashes ########
     69#mail -s "crash run GCM" $address <<ENDMAIL
     70#The run on $machine in $dir has just crashed.
     71#ENDMAIL
    7272############################################""
    7373 exit
    … …  
    7777if ( -f diagfi.nc ) then
    7878  \mv -f diagfi.nc diagfi${numnew}.nc
     79endif
     80if ( -f diagspec.nc ) then
     81  \mv -f diagspec.nc diagspec${numnew}.nc
    7982endif
    8083if ( -f diagsoil.nc ) then
    … …  
    9598
    9699if ( -f rad_bal.out ) then
    97     mv rad_bal.out radbal${numnew}
     100    mv rad_bal.out rad_bal${numnew}.out
    98101endif
    99102
     103if ( -f h2o_bal.out ) then
     104    mv h2o_bal.out h2o_bal${numnew}.out
     105endif
    100106
    101107# If we are over nummax : stop

  • trunk/LMDZ.GENERIC/deftank/run1d.earlymars

    r135 r253  
    11#
    22#-----------------------------------------------------------------------
    3 # Run parameters for the 1D 'testphys1d.e' model                                         
     3# Run parameters for the 1D 'rcm1d.e' model                                         
    44#-----------------------------------------------------------------------                                         
    55

  • trunk/LMDZ.GENERIC/deftank/run1d.earth

    r135 r253  
    11#
    22#-----------------------------------------------------------------------
    3 # Run parameters for the 1D 'testphys1d.e' model                                         
     3# Run parameters for the 1D 'rcm1d.e' model                                         
    44#-----------------------------------------------------------------------                                         
    55
    … …  
    1111time=0
    1212# Number of time steps per sol
    13 day_step=48
     13day_step=24
    1414# Number of sols to run
    15 ndt=2000
     15ndt =2000
    1616
    1717#### Physical parameters
    … …  
    2323# Molar mass of atmosphere (g)
    2424mugaz=28.97
    25 # Specific heat capacity of atmosphere (J kg-1 K-1)
     25# Specific heat capacity of atmosphere?
    2626cpp=1003.0
    2727# latitude (in degrees)
    2828latitude=60.0
    2929# orbital distance at perihelion (AU)
    30 periheli=1.0
     30periastr=1.0
    3131# orbital distance at aphelion (AU)
    32 aphelie=1.0
     32apoastr=1.0
    3333# obliquity (degrees)
    3434obliquit=0.0
    3535
    3636# Albedo of bare ground
    37 albedo=0.1
     37albedo=0.16
    3838# Emissivity of bare ground
    3939emis=1.0
    … …  
    4848# hybrid vertical coordinate ? (.true. for hybrid and .false. for sigma levels)
    4949hybrid=.false.
    50 
     50# autocompute vertical discretisation? (useful for exoplanet runs)
     51autozlevs=.true.
    5152
    5253###### Initial atmospheric temperature profile
    … …  
    7374hauteur=30.
    7475
    75 
    7676# some definitions for the physics, in file 'callphys.def'
    7777INCLUDEDEF=callphys.def

  • trunk/LMDZ.GENERIC/deftank/start_earlymars

    r135 r253  
    77peri_day
    880.0
    9 periheli
     9periastr
    1010227.9
    11 aphelie
     11apoastr
    1212227.9
    1313year_day
    … …  
    2525daysec
    262688775.00
     27albedice
     280.5
     290.5
    2730
    2831surface_mars.nc
    … …  
    3235250
    3336ptot
    34 200000.0
     37500000.0
    3538y
    3639n
    … …  
    38411.0
    3942isotherm
    40 260
     43260.0
    4144co2ice=0
    4245q=0
    43 watercapn
    44 watercaps
     46oborealis
    4547
    4648exit 0

  • trunk/LMDZ.GENERIC/deftank/start_earth

    r135 r253  
    77peri_day
    880.0
    9 periheli
     9periastr
    1010149.0
    11 aphelie
     11apoastr
    1212149.0
    1313year_day
    1414365
     15rad
     166378137.0
    1517omeg
    16187.2722E-05
    … …  
    2628surface_earth.nc
    27290
    28 bilball
    29 0.1
    30 250
    3130ptot
    3231101325.0
    … …  
    3433y
    3534101325.0
    36 emis
    37 1.0
    38 isotherm
    39 300
     35radequi
    4036q=0
    41 watercapn
    42 watercaps
    4337
    4438exit 0

  • trunk/LMDZ.GENERIC/deftank/z2sig.earth

    r135 r253  
    2020   3.6962495e+01
    2121   4.0453872e+01
     22   45
     23   50
     24   55
     25   60
     26   65
     27   70
     28   80
     29   90

  • trunk/LMDZ.GENERIC/deftank/z2sig.mars

    r135 r253  
    1 10.00000     H: atmospheric scale height (km) (used as a reference only)
    2 0.0040       Typical pseudo-altitude (m) for 1st layer (z=H*log(sigma))
     18.00000     H: atmospheric scale height (km) (used as a reference only)
    320.018       ,, ,,  ,, ,, ,, ,,  ,, ,, ,,  2nd layer, etc...
    430.0400
Note: See TracChangeset for help on using the changeset viewer.