|
|
|
GISS AGCM Ocean Boundary Conditions
The pconvert.f FORTRAN code can be used to convert boundary
conditions from one grid resolution to another taking into account the area
weighting of the grid cells on a spherical surface.
This project uses pconvert.f to convert the PRISM2 SST data set from
the supplied 2 X 2 degree (Lat. X Lon.) resolution to a 4 X 5 degree
grid used by the GISS family of GCMs. However, the code is fairly
generic and can be used to convert from any grid resolution to any
other resulotion. The grids must be uniform, except the polar cell is
handled with a special case. Optionally, SICE can be calculated rather
than SST.
The steps used for this conversion process are as follows:
- Download the PRISM
SST data
- Open in Excel and for each tab (jan, feb, ..., dec) save as a
tab-delimeted file (jan.txt, feb.txt, ..., dec.txt)
- If necessary, convert files to the line endings for your
operating system. On OS X, the
mac2unix should be used.
mac2unix *.txt
- Modify the
pconvert.f source to suit your needs. Set
your input and output grid sizes, set the flags to operate using SICE
rather than SST, etc.
- Compile the code. We use Absoft
F77 although the free and cross platform GNU G77
should work.
f77 pconvert.f
- Run the code
./a.out
Converting PRISM to E
Converting 2nd array to EOM data
- The end result should be a
p2e.dat file created in
the current directory. This can be used in Model E
|
