Subject: Re: HAo Visit

From: [email-masked] (Steve Bougher x4900)
Date: Tue, 1 Oct 1996 14:03:58 +0700

----------
X-Sun-Data-Type: text
X-Sun-Data-Description: text
X-Sun-Data-Name: text
X-Sun-Charset: us-ascii
X-Sun-Content-Lines: 15

DEar Ben,

	Attached is a run script for my Mars TGCM.  I typicaly submit
jobs via migs from CEDAR. I have found that migs is better for me
(more versatile) than IRJE. It also allows me to visit NCAR/HAO without
disrupting my filesystems.

	I look forward to modifying the MTGCM and JTGCM codes to :
(1) avoid using nupdate, and (2) use in-core commands, increasing performance
on PAIUTE. Thanks for your help!

--Steve Bougher
  (520)-621-4900
----------
X-Sun-Data-Type: default
X-Sun-Data-Description: default
X-Sun-Data-Name: mtgcmsh5.mgs
X-Sun-Charset: us-ascii
X-Sun-Content-Lines: 261

# QSUB -lT 8000
# QSUB -lM 5Mw
# QSUB -lm 4Mw
# QSUB -q econ
# QSUB -nr
# QSUB -eo
# --------MARS GLOBAL SURVEYOR (1997-1998) VERSION : MTGCMSH5.PA/MGS---
# --------2-HEMISPHERE MESO-EXTENSION MODEL : 78-300 KM 
# --------4-MAJORS (O, CO, N2, CO2); 2-MINORS (AR, O2)
# --------GM  LBC FOR O AND CO FROM KRASNOPOLSKY [1993]
# --------(REVISED GM KZZ=1.5E+07 )
# -------------------------------------------------------------
# --------SMED1: F107=130 / d=1.401 AU / ds= -15.0/EUV=20% /CO2-O=1.5E-12
# --------SMED2: F107=130 / d=1.381 AU / ds= -22.6/EUV=20% /CO2-O=1.5E-12
# --------SLOW1: F107=70  / d=1.401 AU / ds= -15.0/EUV=20% /CO2-O=1.5E-12
# --------SMED3: F107=130 / d=1.659 AU / ds= +22.6/EUV=20% /CO2-O=1.5E-12
# ---*----SMED4: F107=130 / d=1.401 AU / ds= -24.7/EUV=20% /CO2-O=1.5E-12
# -------------------------------------------------------------
# --------NO TIDES BENCHMARK : 4-MAJORS/2-MINORS : CRATES/COMPN2/COMP2
# --------*****************************************************
# --------MGS97E DY:4                    I=SBIMSM  SMED1 : NO TIDES
# --------MGS97L DY:4                    I=SBIMSL  SMED2 : NO TIDES
# --------MGS97T DY:4                    I=SBIMST  SLOW1 : NO TIDES
# --------MGS98L DY:4                    I=SBIMAP  SMED3 : NO TIDES
# --------MGCM2A DY:4                    I=SBIM2A  SMED4 : NO TIDES
# -------------------------------------------------------------
# --------WITH TIDES CASES : Ls=250  SDEC=-22.6
# --------SMED2: DECEMBER 30, 1997 LATE IN MGS AEROBRAKING PHASE
# --------*****************************************************
# --------MGS97L DY:8    OLD TEST#1      (2,2) A=1.0E+05 PH=+1.0
# --------MS97LA DY:4    NEW TEST#2B     (2,2) A=1.5E+05 PH=+1.6
# --------                               (2,4) A=2.0E+05 PH=+2.8
# --------MS97LB DY:4    NEW TEST#2C     (2,2) A=1.5E+05 PH=+1.6
# --------MS97LC DY:0 HR:13 NEW TEST#2A  (2,2) A=3.0E+05 PH=+1.6
# --------              (Watch blowup)   (2,3) A=9.0E+04 PH=-5.3
# --------                               (2,4) A=4.0E+05 PH=+2.8
# -------------------------------------------------------------
# --------WITH TIDES TEST CASES : Ls=217 SDEC=-15
# --------SMED1 & SLOW1 : NOVEMBER 15, 1997 EARLY IN AEROBRAKING PHASE
# --------******************************************************
# --------MS97EA DY:4    NEW TEST#2B     (2,2) A=1.5E+05 PH=+1.6
# --------                               (2,4) A=2.0E+05 PH=+2.8
# --------MS97TA DY:4    NEW TEST#2B     (2,2) A=1.5E+05 PH=+1.6
# --------                               (2,4) A=2.0E+05 PH=+2.8
# --------MSGCM1 DY:4    GCM95.06 TEST#1
# --------                               (2,2) A=1.1E+05 PH=
# --------                               (2,3) A=0.2E+05 PH=
# --------                               (2,4) A=0.4E+05 PH=
# --------                               (2,5) A=0.1E+05 PH=
# --------                               (2,6) A=0.2E+05 PH=
# --------(COMPARE TO ORIGINAL SMED1 ABOVE)
# -------------------------------------------------------------
# --------WITH TIDES TEST CASES : Ls=278 SDEC=-24.7
# --------SMED4 FEB.20, 1998 EPOCH
# --------******************************************************
# --------MGCM2A DY:8     GCM96_01M (#2A) SEPARATE PHASES 
# --------MGCM2B DY:4     GCM96_01M (#2B) 
# --------MGCM2C DY:4     GCM96_01M (#2C) 
# --------MGCM2D DY:4     GCM96_01M (#2D) 
# --------******************************************************
# --------MGCM2A DY:12    GCM96_01M (#2A) ALL SAME PHASES AS (2,2)
# --------MGCM2B DY:8     GCM96_01M (#2B) 
# --------MGCM2D DY:8     GCM96_01M (#2D) 
# --------(COMPARE TO ORIGINAL SMED1 ABOVE)
# --------(COMPARE TO GCM95.06 ABOVE)
# --------******************************************************
ja
set echo timestamp
# --setup local tmp space on shavano ---------------
set dir = $TMPDIR
cd $dir
echo "$TMPDIR" > tempdir
fil.back tempdir "$QSUB_REQNAME".dir
# -------------------------------------------------------------
# --setup plib history needed ----------------------
set locl = "/usr/tmp/MTGCM2/SWBPST.MO"
set remot = "/BOUGHER/SWBPST/MTGCM5SH"
if (! -d /usr/tmp/MTGCM2) mkdir /usr/tmp/MTGCM2
if (! -f /usr/tmp/MTGCM2/SWBPST.MO) then
  msread -f TR $locl $remot
endif
ln $locl pl
# -------------------------------------------------------------
# --pull in small files needed from SHAVANO
cp ~/sflux/FLX130N iflx
#cp ~/sflux/F76200N iflx
# -------------------------------------------------------------
# -- mod decks via update utility ---------------
unset timestamp
cat > pdir << 'EOF'
*/
*/ MTGCM5SH MODS FOR (1) PARMS FOR MGS97-98 RUNS 
*/
*ID PARMS
*DELETE CON.204,212 
C ****  GM VALUE MAX EDDY COEFFICIENT : TURBOPAUSE (DIFF/VISCOSITY/COND)
      C(27)=1.5E+07 
*DELETE DTUV.198,200
C **** MESO LOWER LIMIT KZZ=1.0E+07
      IF(VERT.LE.0.667)VERT=0.667
*DELETE CON.232,235 
C ****  GM  LBC FOR O AND CO FROM KRASNOPOLSKY [1993; ICARUS, 101, p.33-41]. 
C ***   AT 80 KM (DOUBLED).
      C(46)=2.*2.500E-04*(16./44.)
      C(47)=2.*1.000E-03*(28./44.)
*DELETE QRJINT.16,19
C **** SOLAR DECLINATION IN DEGREES (SEASON) MGS SMED3 RUN
C     DATA SDEC/+22.6/  
C **** SOLAR DECLINATION IN DEGREES (SEASON) MGS SMED2 RUN
C     DATA SDEC/-22.6/  
C **** SOLAR DECLINATION IN DEGREES (SEASON) MGS SMED4 RUN
      DATA SDEC/-24.7/  
C **** SOLAR DECLINATION IN DEGREES (SEASON) MGS SMED1 RUN
C **** SOLAR DECLINATION IN DEGREES (SEASON) MGS SLOW1 RUN
C     DATA SDEC/-15.0/  
*DELETE IONPCE.29,30
C **** SOLAR DECLINATION IN DEGREES (SEASON) MGS SMED3 RUN
C     DATA SDEC/+22.6/  
C **** SOLAR DECLINATION IN DEGREES (SEASON) MGS SMED2 RUN
C     DATA SDEC/-22.6/  
C **** SOLAR DECLINATION IN DEGREES (SEASON) MGS SMED4 RUN
      DATA SDEC/-24.7/  
C **** SOLAR DECLINATION IN DEGREES (SEASON) MGS SMED1 RUN
C **** SOLAR DECLINATION IN DEGREES (SEASON) MGS SLOW1 RUN
C     DATA SDEC/-15.0/  
*DELETE VDIS.34,41
C ---------------- EUV(1-78)---------------------------------
C **** SOLAR MAX FLUXES (F107=200) READIN FROM SH: ~/sflux/F78348N 
C **** SOLAR MED FLUXES (F107=130) READIN FROM SH: ~/sflux/FLX130N
C **** SOLAR MIN FLUXES (F107=070) READIN FROM SH: ~/sflux/F76200N
C -----------------EUV(79-80)--------------------------------
C **** XRAY FLUXES ADDED AS NEEDED : F107=200/130/070
C     DATA EUV(79)/1.0022/,EUV(80)/0.3200/
      DATA EUV(79)/0.5601/,EUV(80)/0.1788/
C     DATA EUV(79)/0.1180/,EUV(80)/0.0376/
*DELETE VDIS.149,150
C **** HELIOCENTRIC SCALING FACTOR : MGS SMED3 RUN
C     SC=0.3633
C **** HELIOCENTRIC SCALING FACTOR : MGS SMED2 RUN
C     SC=0.52434
C **** HELIOCENTRIC SCALING FACTOR : MGS SMED1 RUN
C **** HELIOCENTRIC SCALING FACTOR : MGS SLOW1 RUN
C **** HELIOCENTRIC SCALING FACTOR : MGS SMED4 RUN
      SC=0.50948
*DELETE VDIS.174,179
C **** JCO2 IS STRAIGHT W/O ENHANCEMENT FACTOR FOR ICO2
C **** FULL IONPCE IMPACT ON O AND CO SRCS AND SINKS
*/
'EOF'
set timestamp
#
#--- edit/compile/load prog + libs/run --------------------
#
#--- edit update library
nupdate -p pl -i pdir -c mtgcm2 -f -o id,ns,um
rm pl pdir
#
# --assign formated  input/output ---------------------------
assign -a iflx fort.17
#---  compile on cft77 
cft77 -m 4 -e zi -a static -i 64 -l mtgcm2.cft mtgcm2.f
#fil.back mtgcm2.cft "$QSUB_REQNAME".cft
#--- *  debug (method 1) on cft77 * ------------------------ 
#cft77 -m 4 -o off -e sz -a static -i 64 -l mtgcm2.cft mtgcm2.f
#fil.back mtgcm2.cft "$QSUB_REQNAME".cft
#----------------------------------------------------------- 
#
#---  load using segldr
segldr -f zeros -o "$QSUB_REQNAME".out -i ~/segdir mtgcm2.o
# segldr -f indef -o "$QSUB_REQNAME".out -i ~/segdir mtgcm2.o
rm mtgcm2.f mtgcm2.o
unset timestamp
cp "$QSUB_REQNAME".out a.out
a.out  << 'EOF'
C*****************************************************************
C***  UNICOS TIMING: ON SH
C***  MTGCMSH5.MGS takes 134 sec to compile, and about 2.8 sec    ***
C***   per time step thereafter.  
C***   (A)     1800 sec per SOL     04/26/95    4-MAJORS/2-MINORS  *****
C***                                            (NEW COMP2 ADDED)
C*******************************************************************
C***  UNICOS TIMING: ON PA
C***  MTGCMSH5.AL takes 200 sec to start up, and about 5.1 sec   ***
C***   per time step thereafter.  
C***   (A)     3200 sec per SOL     10/XX/95    4-MAJORS/2-MINORS  *****
C***                                            (NEW COMP2 ADDED)
C*******************************************************************
C STEP=150 : ITER=576/SOL : 
C MAX SIZE FOR HISTORY ON MSS: 180 MB (OR 9-HISTORY SETS ONLY)
C*******************************************************************
 DATA='BOUGHER','SWBM96','SBIM2A'
 OUTPUT='BOUGHER','SWBM96','MGCM2D'
 START=4,0,0
 STOP=8,0,0
 HIST=1,0,0
 STEP=150
 SAVE=1,0,0
 DATE=1998,043
 LABEL='05/06/96 MTGCMSH5 : MGS-MESO'
 DISPOSE=1
 AURORA=0
 MAG=0.,0.,0.,0.
 POWER=0.
 CTPOTEN=0.
 BYIMF=0.
 DIFHOR=0
 IUIVI=0
 TIDE=8.2E+4,3.0E+4,2.0E+4,1.1E+4,1.5E+4,-1.1,-1.1,-1.1,-1.1,-1.1
 IPOWER=0
 ENDOFREAD
'EOF'
set timestamp
#---   * debug (method 1) on cft77 * ----------------------- 
#debug > mmgcm.dbg
#fil.back mmgcm.dbg "$QSUB_REQNAME".dbg
#----------------------------------------------------------- 
#  [GCM96_01M]
#  GCM2D  : DUST STORM 20-SOL INTEGRATION (TAU=0-2.5)
#  TIDE=8.2E+4,3.0E+4,2.0E+4,1.1E+4,1.5E+4,-1.1,-1.2,+0.8,+0.1,+1.5
#  TIDE=8.2E+4,3.0E+4,2.0E+4,1.1E+4,1.5E+4,-1.1,-1.1,-1.1,-1.1,-1.1
#  ------
#  GCM2C  : DUST STORM 20-SOL INTEGRATION (TAU=0-2.5)
#  TIDE=1.0E+5,4.4E+4,2.5E+4,1.8E+4,2.1E+4,+1.0,-1.3,-0.3,+0.4,+0.7
#  TIDE=1.0E+5,4.4E+4,2.5E+4,1.8E+4,2.1E+4,+1.0,+1.0,+1.0,+1.0,+1.0
#  ------
#  GCM2B  : DUST STORM 20-SOL INTEGRATION (TAU=0-2.5)
#  TIDE=1.1E+5,4.4E+4,3.8E+4,3.2E+4,2.6E+4,+0.3,-0.6,-0.9,+0.3,+0.1
#  TIDE=1.1E+5,4.4E+4,3.8E+4,3.2E+4,2.6E+4,+0.3,+0.3,+0.3,+0.3,+0.3
#  ------
#  GCM2A  : DUST STORM 20-SOL INTEGRATION (TAU=0-2.5)
#  TIDE=6.6E+4,3.0E+4,2.5E+4,1.8E+4,1.3E+4,+1.6,-5.1,+4.7,-5.6,+5.8
#  TIDE=6.6E+4,3.0E+4,2.5E+4,1.8E+4,1.3E+4,+1.6,+1.6,+1.6,+1.6,+1.6
#  ---------------------------------------------------------
#  [GCM95.06]
#  GCM1 TAU=1.0 GLOBALLY UNIFORM DUST TAU=1.0
#  TIDE=1.1E+5,0.2E+5,0.4E+5,0.1E+5,0.2E+5,+2.1,-1.2,+3.5,+5.6,+3.0
#  ---------------------------------------------------------
#  TEST#2C TAU=1.1 NON DUST-STORM CASE (NEW)
#  TIDE=1.5E+5,0.,0.,0.,0.,+1.6,0.,0.,0.,0.
#  ---------------------------------------------------------
#  TEST#2B TAU=1.1 NON DUST-STORM CASE (NEW)
#  TIDE=1.5E+5,0.,2.0E+5,0.,0.,+1.6,0.,+2.8,0.,0.
#  ---------------------------------------------------------
#  TEST#2A TAU=1.1 NON DUST-STORM CASE (NEW)
#  TIDE=3.0E+5,9.0E+4,4.0E+5,0.,0.,+1.6,-5.3,+2.8,0.,0.
#  ---------------------------------------------------------
#  TEST#1 TAU=1.5 MODERATE DUST-STORM CASE (OLD)
#  TIDE=1.0E+5,0.,0.,0.,0.,+1.0,0.,0.,0.,0.
#  ---------------------------------------------------------
#  TIDE=0.,0.,0.,0.,0.,0.,0.,0.,0.,0.
#
# --transfer extra data files back to cedar -----------------------------
#  QSUB_REQNAME automatically comes back ($OUT) to front-end --
#
#  TIDAL LBC CHECK :fil.back file.dat "$QSUB_REQNAME".dat
fil.back gmeta "$QSUB_REQNAME".gmeta
#
# -------------------------------------------------------------
#
#
ja -cslst