C C Test low_proton routine, called from orora: C CLABEL='TIEGCM EQUINOX.P7 PLANET COMPARISON 3/24/97 ' LABEL='TIEGCM TEST LOW_PROTON 5/98 ' DATE=1997,80 COUTPUT='ROBLE','RGR97','TEN2X10','ROBLE','RGR97','TEN2X11', C 'ROBLE','RGR97','TEN2X12','ROBLE','RGR97','TEN2X13', C 'ROBLE','RGR97','TEN2X14','ROBLE','RGR97','TEN2X15', C 'ROBLE','RGR97','TEN2X16','ROBLE','RGR97','TEN2X17', C 'ROBLE','RGR97','TEN2X18','ROBLE','RGR97','TEN2X19', C 'ROBLE','RGR97','TEN2X20','ROBLE','RGR97','TEN2X21', C 'ROBLE','RGR97','TEN2X22','ROBLE','RGR97','TEN2X23' SOURCE = 'ROBLE','RGR97','TEN2X19',10,0,0 OUTPUT = 'FOSTER','PROTON1','TEN2X19' MAGVOL='ECRIDLEY','ECR90','ECRMG6' START= 10,0,0 STOP= 10,0,3 HIST= 0,0,3 SAVE=0,0,3 STEP= 180 INCORE = 1 C C Secondary histories (comment all SECxxxx out if sec hist not needed): C (these histories contain a max of 29 fields at the current time step) C C SECOUT = mss secondary history output volumes C SECHIST = disk write frequency of secondary histories (day:hr:min) C SECSAVE = mss write frequency of secondary history files (day:hr:min) C SECSTART,STOP = start,stop times for writing secondary histories C (day:hr:min) C COUTPUT = 'FOSTER','PROTON','TEN2X19' CSECOUT = 'FOSTER','SPROTON','TEN2X19' CSECSTART = 10,0,0 CSECSTOP = 10,0,3 CSECHIST = 0,0,3 CSECSAVE = 0,0,3 C C SECFLDS = list of character names of secondary history fields. C (may be primary fields or user-specified by addfsech calls) C Below is list of primary history fields on tgcm13: C C /'TN ','UN ','VN ','O2 ','O1 ', C 'N4S ','NO ','O+ ','N2D ','TI ', C 'TE ','NE ','O2+ ','W ','Z ', C 'POTEN ' C C QIA[1-5]_xx is output ion production for N2+, O2+, O+, NO+, N+ C CSECFLDS = 'XNO2','XNO','XNN2','XRHO' CSECFLDS = 'QIA1_N2+','QIA2_O2+','QIA3_O+','QIA5_N+', C 'XNO2','XNO','XNN2','XRHO','ZPHT' CSECFLDS = 'EFLUX_LP' CSECFLDS = 'QIA1_N2+','QIA2_O2+','QIA3_O+','QIA4_NO+','QIA5_N+', C 'NQO2P','NQOP','NQN2P','NQNP' CSECFLDS = 'S7:T6*T8','S8:T7*T9','S9','S10' C C*****************************P CDISPOSE=1 DISPOSE=0 C*****************************P MAG=-74.5,127.,79.,-70. DIFHOR=0 IUIVI=1 CTIDE=0.,0.,0.,0.,0.0,0.0,0.0,0.0,0.0,0.0 C JANUARY TIDES FORBES(NEW) CTIDE=2.8E+4,0.7E+4,1.0E+4,1.5E+4,0.6E+4,10.4,1.2,7.8,9.0,8.3 C FEBRUARY TIDES FORBES(NEW) CTIDE=2.2E+4,0.7E+4,1.1E+4,1.3E+4,0.3E+4,10.2,0.8,8.1,8.8,7.8 C MARCH TIDES FORBES(NEW) TIDE=1.5E+4,0.1E+4,2.8E+4,0.1E+4,0.2E+4,10.3,5.8,7.6,8.3,8.4 C APRIL TIDES FORBES(NEW) CTIDE=1.8E+4,0.3E+4,2.2E+4,0.6E+4,0.3E+4,10.2,5.6,7.7,1.4,6.2 C MAY TIDES FORBES(NEW) CTIDE=1.9E+4,0.3E+4,0.5E+4,0.9E+4,0.5E+4,10.3,8.5,7.7,2.2,7.4 C JUNE TIDES FORBES(NEW) CTIDE=2.3E+4,0.5E+4,0.4E+4,0.8E+4,0.4E+4,10.5,8.6,6.6,3.0,8.3 C JULY TIDES FORBES(NEW) CTIDE=2.3E+4,0.8E+4,0.7E+4,1.1E+4,0.6E+4,10.2,10.5,3.8,3.2,8.1 C AUGUST TIDES FORBES(NEW) CTIDE=2.2E+4,1.1E+4,0.5E+4,1.3E+4,0.6E+4,10.2,10.7,2.1,2.7,7.4 C SEPTEMBER TIDES FORBES(NEW) CTIDE=1.5E+4,0.4E+4,1.4E+4,0.9E+4,0.6E+4,10.4,1.6,8.2,1.3,0.8 C OCTOBER TIDES FORBES(NEW) CTIDE=1.7E+4,0.4E+4,1.8E+4,0.2E+4,0.4E+4,10.3,3.4,7.7,4.0,2.9 C NOVEMBER TIDES FORBES(NEW) CTIDE=2.4E+4,0.4E+4,1.2E+4,1.0E+4,0.7E+4,10.5,1.8,8.0,8.5,6.8 C DECEMBER TIDES FORBES(NEW) CTIDE=2.7E+4,0.7E+4,0.9E+4,1.2E+4,0.5E+4,10.5,0.9,7.4,8.8,8.1 CTIDE2=4.0E+4,4.0 CTIDE2=1.0E+4,4.0 TIDE2=0.0,0.0 TIDEANN=0 AURORA=1 COLFAC=1.5 ENDOFREAD F107=200. F107A=200. CF107=130. CF107A=130. CF107=70. CF107A=70. CAMIVOL='EMERY','PVIEN','S86F3B' OLDALF=1 C H POWER = 2 3 5 7 11 16 23 33 48 82 : 115 155 255 C CTPOTEN = 20 30 40 50 60 70 80 90 100 115 : 130 150 200 POWER=5. CTPOTEN=45. CPOWER=23 CCTPOTEN=80. BYIMF=0. ENDOFREAD