65 |
|
|
66 |
c == end of interface == |
c == end of interface == |
67 |
|
|
68 |
|
#ifdef ALLOW_CLIM_CYCLIC |
69 |
cph( |
c record numbers are assumed 1 to 12 corresponding to |
70 |
cph Routine exf_GetSSTClimRec has standard cyclic forcing |
c Jan. through Dec. |
71 |
cph assuming 12 monthly climatological records on file |
call cal_GetMonthsRec( |
72 |
cph with irec=1: Jan, ... irec=12: Dec |
O fac, first, changed, |
73 |
cph Will be included as optional choice of to the above. |
O count0, count1, |
74 |
cph |
I mycurrenttime, mycurrentiter, mythid |
75 |
cph call exf_GetSSTClimRec( |
& ) |
76 |
cph O fac, first, changed, |
#else |
|
cph O count0,count1, |
|
|
cph I mytime, myiter, mythid |
|
|
cph & ) |
|
|
cph) |
|
|
|
|
|
|
|
77 |
c get record numbers and interpolation factor for climsst |
c get record numbers and interpolation factor for climsst |
78 |
call exf_GetFFieldRec( |
call exf_GetFFieldRec( |
79 |
I climsststartdate, climsstperiod |
I climsststartdate, climsstperiod |
81 |
O , count0, count1 |
O , count0, count1 |
82 |
I , mycurrenttime, mycurrentiter, mythid |
I , mycurrenttime, mycurrentiter, mythid |
83 |
& ) |
& ) |
84 |
|
#endif |
85 |
|
|
86 |
|
#ifndef ALLOW_AUTODIFF_TAMC |
87 |
if ( first ) then |
if ( first ) then |
88 |
|
#endif |
89 |
call mdsreadfield( climsstfile, exf_clim_iprec |
call mdsreadfield( climsstfile, exf_clim_iprec |
90 |
& , exf_clim_yftype, 1 |
& , exf_clim_yftype, 1 |
91 |
& , climsst1, count0, mythid |
& , climsst1, count0, mythid |
95 |
else |
else |
96 |
call exf_filter_rs( climsst1, climsstmask, mythid ) |
call exf_filter_rs( climsst1, climsstmask, mythid ) |
97 |
end if |
end if |
98 |
|
#ifndef ALLOW_AUTODIFF_TAMC |
99 |
endif |
endif |
100 |
|
#endif |
101 |
|
|
102 |
|
#ifndef ALLOW_AUTODIFF_TAMC |
103 |
if (( first ) .or. ( changed )) then |
if (( first ) .or. ( changed )) then |
104 |
|
#endif |
105 |
call exf_SwapFFields( climsst0, climsst1, mythid ) |
call exf_SwapFFields( climsst0, climsst1, mythid ) |
106 |
|
|
107 |
call mdsreadfield( climsstfile, exf_clim_iprec |
call mdsreadfield( climsstfile, exf_clim_iprec |
113 |
else |
else |
114 |
call exf_filter_rs( climsst1, climsstmask, mythid ) |
call exf_filter_rs( climsst1, climsstmask, mythid ) |
115 |
end if |
end if |
116 |
|
#ifndef ALLOW_AUTODIFF_TAMC |
117 |
endif |
endif |
118 |
|
#endif |
119 |
|
|
120 |
c Loop over tiles. |
c Loop over tiles. |
121 |
do bj = mybylo(mythid),mybyhi(mythid) |
do bj = mybylo(mythid),mybyhi(mythid) |