1 |
#include "EXF_CPPOPTIONS.h" |
#include "EXF_OPTIONS.h" |
2 |
|
|
3 |
subroutine exf_set_gen( |
subroutine exf_set_gen( |
4 |
& genfile, genstartdate, genperiod, exf_inscal_gen, |
& genfile, genstartdate, genperiod, exf_inscal_gen, |
40 |
|
|
41 |
c == routine arguments == |
c == routine arguments == |
42 |
|
|
43 |
integer genstartdate(4) |
_RL genstartdate, genperiod |
|
_RL genperiod |
|
44 |
_RL exf_inscal_gen |
_RL exf_inscal_gen |
45 |
_RL genfld(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy) |
_RL genfld(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy) |
46 |
_RL gen0 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy) |
_RL gen0 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy) |
71 |
_RL fac |
_RL fac |
72 |
|
|
73 |
integer bi, bj |
integer bi, bj |
74 |
integer i, j |
integer i, j, interp_method |
75 |
|
|
76 |
c == end of interface == |
c == end of interface == |
77 |
|
|
87 |
|
|
88 |
if ( first ) then |
if ( first ) then |
89 |
#ifdef USE_EXF_INTERPOLATION |
#ifdef USE_EXF_INTERPOLATION |
90 |
|
interp_method = 2 |
91 |
call exf_interp( genfile, exf_iprec |
call exf_interp( genfile, exf_iprec |
92 |
& , gen1, count0, gen_xout, gen_yout |
& , gen1, count0, gen_xout, gen_yout |
93 |
& , gen_lon0,gen_lon_inc |
& , gen_lon0,gen_lon_inc |
94 |
& , gen_lat0,gen_lat_inc |
& , gen_lat0,gen_lat_inc |
95 |
& , gen_nlon,gen_nlat,mythid |
& , gen_nlon,gen_nlat,interp_method,mythid |
96 |
& ) |
& ) |
97 |
#else |
#else |
98 |
call mdsreadfield( genfile, exf_iprec, exf_yftype, 1 |
call mdsreadfield( genfile, exf_iprec, exf_yftype, 1 |
111 |
call exf_SwapFFields( gen0, gen1, mythid ) |
call exf_SwapFFields( gen0, gen1, mythid ) |
112 |
|
|
113 |
#ifdef USE_EXF_INTERPOLATION |
#ifdef USE_EXF_INTERPOLATION |
114 |
|
interp_method = 2 |
115 |
call exf_interp( genfile, exf_iprec |
call exf_interp( genfile, exf_iprec |
116 |
& , gen1, count1, gen_xout, gen_yout |
& , gen1, count1, gen_xout, gen_yout |
117 |
& , gen_lon0,gen_lon_inc |
& , gen_lon0,gen_lon_inc |
118 |
& , gen_lat0,gen_lat_inc |
& , gen_lat0,gen_lat_inc |
119 |
& , gen_nlon,gen_nlat,mythid |
& , gen_nlon,gen_nlat,interp_method,mythid |
120 |
& ) |
& ) |
121 |
#else |
#else |
122 |
call mdsreadfield( genfile, exf_iprec, exf_yftype, 1 |
call mdsreadfield( genfile, exf_iprec, exf_yftype, 1 |