#! /home/benno/multsfc/prg.ingrid/meqsfc2F
WINDSTR  (/home/keithr/public_html/band_i.cdf) readCDF
DENSITY  /fox/room/alexeyk/nwind/pacific.pro
ENDFILES                                
  0            =      NSTART  model start type(2=read from hist)    
 0.5           =      TSTART
  1	       =      MODELT      ! this flag chooses a force routine
Banded zonal pulse; One mode.
  1	       =      NUMMOD  number of modes 
  0	       =      HEQUIVF 0 means calculate modes from a profile
  1	       =      NSEG    number of horizontal segments
 .false.       =      CHKSUM  
  -1	       =      NPRINT  print every 3'rd timestep
  00	       =      NREWND  no rewind
  0	       =      NTAPE   write history after end
  -1000.       =      TPLMIN  put data out to SUMRY immediately
   3000.	       =      TDECAY1  time rate of decay for mode 1
  78	       =      NX      number of zonal grid points
  115	       =      NY      number of meridional grid points
  20	       =      NZ      number of points in depth
   0.5         =     TZERO                     ! start dec 1969
  172.5         =      TENDD                     ! stop  dec 1992
  .125          =      DTD     time step
  124.	       =      XWD     western limit
  280.	       =      XED     eastern limit
  -28.75       =      YSD     southern limit
   28.75       =      YND     northern limit
Depth          =      ZIS     vertical coord(could also be temperature)
  0.	       =      ZTOP    top limit
  1000.        =      ZBOT    bottom limit
  0	       =      MASKW   wind masking flag(for Greens's fn. studies)
  11	       =      MWNDGF 
  -60.	       =      WMXW   
   15.	       =      WMXE   
  -30.	       =      WMYS   
   30.	       =      WMYN   
 50.	       =      HSFC   
 2.	       =      RSINV   frictional decay time-suface layer(days)
  .333	       =      BCSCL   Eckman Horizontal Boundary Scale
 10000.	       =      ZCUT   Cut off depth for modes(1000 effectively disables)
  28.75	       =      YNORTH1  Sub basin information
  -28.75       =      YSOUTH1
  124.         =      XWEST1   
    0	       =      NATM    atmospheric model type(0= no atmosph. model)

184.3988210, .5275161, -.0186502 1	B0 , B1 , B2 Thermocline depth param
5 					Number of cubic spline coefficients
0.0,34.613866,50.0,123.871529,210.0,240.0				XK
17.076101,17.512403,19.958565,26.155706,28.846940			YI
 0.022973,-0.003364, 0.000089		 	COEFS in NCOEF lines of 3
 0.108295, 0.005829,-0.000165
 0.170686,-0.001774, 0.000008
 0.041315, 0.000022,-0.000002
 0.009193,-0.000395, 0.000138
26.  17.1 	HCUT , TCUT (minimum temperature for shallow thermocline)
 .5 .62 .0019  1.5 0.    4. 7.071 .72 0. GAMMA_ACLD AALPHA RH ASST CT0_UVMIN UNDRG RH HA
comments:
TDEEP parameterization to match Richard in the Pacific
/description (Cane - Patton Linear Equatorial Model \(JPO, 14, 1984\)) def

Richard's parameterization of 50m temperature 
      as a function of thermocline depth
Quadratic parameterization of thermocline depth
      as a function of model thermocline depth
 
Richard's heat flux parameters
:comments
SST plots -- equatorial pacific
.
\begin{ingrid}
GL
/XOVY 1.6 def
/plotsize 1200 def
SOURCES
MODEL 

H
%T 60.5 last 1. RANGESTEP
% Y 2. STEP % X Y 1 SM121
/description (Cane - Patton Linear Equatorial Model \(JPO, 14, 1984\)) def
dup

Z 0. VALUE Y -15 15 RANGE
/name (Sea_Level) cvn def
/fullname name cvntos def
%100. mul 
35. add
DATA 31 39 0.25 RANGESTEP
X Y COLOR
%X Y COLORCONTOUR
%(SL.cdf) writeCDF

%Z 250. VALUE
%/name (Thermocline_Depth) cvn def
%/fullname name cvntos def
%(/breadbox/wheat/alexeyk/CPmodel/THD2.cdf) writeCDF

\end{ingrid}

Local Variables:
compile-command: "sh 4senq2.in"
End: