/[MITgcm]/MITgcm/pkg/ctrl/ctrl_init.F
ViewVC logotype

Diff of /MITgcm/pkg/ctrl/ctrl_init.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.21 by heimbach, Thu Jan 5 17:05:42 2006 UTC revision 1.22 by gforget, Wed Mar 29 22:07:34 2006 UTC
# Line 187  c--   Heat flux. Line 187  c--   Heat flux.
187          call cal_TimePassed( xx_hfluxstartdate, modelstartdate,          call cal_TimePassed( xx_hfluxstartdate, modelstartdate,
188       &                       difftime, mythid )       &                       difftime, mythid )
189          call cal_ToSeconds ( difftime, diffsecs, mythid )          call cal_ToSeconds ( difftime, diffsecs, mythid )
190            if ( xx_hfluxperiod .EQ. 0 ) then
191            startrec=1
192            endrec=12
193            else
194          startrec = int((modelstart + startTime - diffsecs)/          startrec = int((modelstart + startTime - diffsecs)/
195       &                 xx_hfluxperiod) + 1       &                 xx_hfluxperiod) + 1
196          endrec   = int((modelend + startTime - diffsecs + modelstep/2)/          endrec   = int((modelend + startTime - diffsecs + modelstep/2)/
197       &                 xx_hfluxperiod) + 2       &                 xx_hfluxperiod) + 2
198            endif
199  # else  # else
200          startrec = 1          startrec = 1
201          endrec   = 1          endrec   = 1
# Line 209  c--   Atmos. temperature Line 214  c--   Atmos. temperature
214          call cal_TimePassed( xx_atempstartdate, modelstartdate,          call cal_TimePassed( xx_atempstartdate, modelstartdate,
215       &                       difftime, mythid )       &                       difftime, mythid )
216          call cal_ToSeconds ( difftime, diffsecs, mythid )          call cal_ToSeconds ( difftime, diffsecs, mythid )
217            if ( xx_atempperiod .EQ. 0 ) then
218            startrec=1
219            endrec=12
220            else
221          startrec = int((modelstart + startTime - diffsecs)/          startrec = int((modelstart + startTime - diffsecs)/
222       &                 xx_atempperiod) + 1       &                 xx_atempperiod) + 1
223          endrec   = int((modelend + startTime - diffsecs + modelstep/2)/          endrec   = int((modelend + startTime - diffsecs + modelstep/2)/
224       &                 xx_atempperiod) + 2       &                 xx_atempperiod) + 2
225            endif
226  # else  # else
227          startrec = 1          startrec = 1
228          endrec   = 1          endrec   = 1
# Line 241  c--   Salt flux. Line 251  c--   Salt flux.
251          call cal_TimePassed( xx_sfluxstartdate, modelstartdate,          call cal_TimePassed( xx_sfluxstartdate, modelstartdate,
252       &                       difftime, mythid )       &                       difftime, mythid )
253          call cal_ToSeconds ( difftime, diffsecs, mythid )          call cal_ToSeconds ( difftime, diffsecs, mythid )
254            if ( xx_sfluxperiod .EQ. 0 ) then
255            startrec=1
256            endrec=12
257            else
258          startrec = int((modelstart + startTime - diffsecs)/          startrec = int((modelstart + startTime - diffsecs)/
259       &                 xx_sfluxperiod) + 1       &                 xx_sfluxperiod) + 1
260          endrec   = int((modelend + startTime - diffsecs + modelstep/2)/          endrec   = int((modelend + startTime - diffsecs + modelstep/2)/
261       &                 xx_sfluxperiod) + 2       &                 xx_sfluxperiod) + 2
262            endif
263  # else  # else
264          startrec = 1          startrec = 1
265          endrec   = 1          endrec   = 1
# Line 263  c--   Atmos. humidity Line 278  c--   Atmos. humidity
278          call cal_TimePassed( xx_aqhstartdate, modelstartdate,          call cal_TimePassed( xx_aqhstartdate, modelstartdate,
279       &                       difftime, mythid )       &                       difftime, mythid )
280          call cal_ToSeconds ( difftime, diffsecs, mythid )          call cal_ToSeconds ( difftime, diffsecs, mythid )
281            if ( xx_aqhperiod .EQ. 0 ) then
282            startrec=1
283            endrec=12
284            else
285          startrec = int((modelstart + startTime - diffsecs)/          startrec = int((modelstart + startTime - diffsecs)/
286       &                 xx_aqhperiod) + 1       &                 xx_aqhperiod) + 1
287          endrec   = int((modelend + startTime - diffsecs + modelstep/2)/          endrec   = int((modelend + startTime - diffsecs + modelstep/2)/
288       &                 xx_aqhperiod) + 2       &                 xx_aqhperiod) + 2
289            endif
290  # else  # else
291          startrec = 1          startrec = 1
292          endrec   = 1          endrec   = 1
# Line 295  c--   Zonal wind stress. Line 315  c--   Zonal wind stress.
315          call cal_TimePassed( xx_tauustartdate, modelstartdate,          call cal_TimePassed( xx_tauustartdate, modelstartdate,
316       &                       difftime, mythid )       &                       difftime, mythid )
317          call cal_ToSeconds ( difftime, diffsecs, mythid )          call cal_ToSeconds ( difftime, diffsecs, mythid )
318            if ( xx_tauuperiod .EQ. 0 ) then
319            startrec=1
320            endrec=12
321            else
322          startrec = int((modelstart + startTime - diffsecs)/          startrec = int((modelstart + startTime - diffsecs)/
323       &                 xx_tauuperiod) + 1       &                 xx_tauuperiod) + 1
324          endrec   = int((modelend + startTime - diffsecs + modelstep/2)/          endrec   = int((modelend + startTime - diffsecs + modelstep/2)/
325       &                 xx_tauuperiod) + 2       &                 xx_tauuperiod) + 2
326            endif
327  # else  # else
328          startrec = 1          startrec = 1
329          endrec   = 1          endrec   = 1
# Line 317  c--   Zonal wind speed. Line 342  c--   Zonal wind speed.
342          call cal_TimePassed( xx_uwindstartdate, modelstartdate,          call cal_TimePassed( xx_uwindstartdate, modelstartdate,
343       &                       difftime, mythid )       &                       difftime, mythid )
344          call cal_ToSeconds ( difftime, diffsecs, mythid )          call cal_ToSeconds ( difftime, diffsecs, mythid )
345            if ( xx_uwindperiod .EQ. 0 ) then
346            startrec=1
347            endrec=12
348            else
349          startrec = int((modelstart + startTime - diffsecs)/          startrec = int((modelstart + startTime - diffsecs)/
350       &                 xx_uwindperiod) + 1       &                 xx_uwindperiod) + 1
351          endrec   = int((modelend + startTime - diffsecs + modelstep/2)/          endrec   = int((modelend + startTime - diffsecs + modelstep/2)/
352       &                 xx_uwindperiod) + 2       &                 xx_uwindperiod) + 2
353            endif
354  # else  # else
355          startrec = 1          startrec = 1
356          endrec   = 1          endrec   = 1
# Line 349  c--   Meridional wind stress. Line 379  c--   Meridional wind stress.
379          call cal_TimePassed( xx_tauvstartdate, modelstartdate,          call cal_TimePassed( xx_tauvstartdate, modelstartdate,
380       &                       difftime, mythid )       &                       difftime, mythid )
381          call cal_ToSeconds ( difftime, diffsecs, mythid )          call cal_ToSeconds ( difftime, diffsecs, mythid )
382            if ( xx_tauvperiod .EQ. 0 ) then
383            startrec=1
384            endrec=12
385            else
386          startrec = int((modelstart + startTime - diffsecs)/          startrec = int((modelstart + startTime - diffsecs)/
387       &                 xx_tauvperiod) + 1       &                 xx_tauvperiod) + 1
388          endrec   = int((modelend + startTime - diffsecs + modelstep/2)/          endrec   = int((modelend + startTime - diffsecs + modelstep/2)/
389       &                 xx_tauvperiod) + 2       &                 xx_tauvperiod) + 2
390            endif
391  # else  # else
392          startrec = 1          startrec = 1
393          endrec   = 1          endrec   = 1
# Line 371  c--   Meridional wind speed. Line 406  c--   Meridional wind speed.
406          call cal_TimePassed( xx_vwindstartdate, modelstartdate,          call cal_TimePassed( xx_vwindstartdate, modelstartdate,
407       &                       difftime, mythid )       &                       difftime, mythid )
408          call cal_ToSeconds ( difftime, diffsecs, mythid )          call cal_ToSeconds ( difftime, diffsecs, mythid )
409            if ( xx_vwindperiod .EQ. 0 ) then
410            startrec=1
411            endrec=12
412            else
413          startrec = int((modelstart + startTime - diffsecs)/          startrec = int((modelstart + startTime - diffsecs)/
414       &                 xx_vwindperiod) + 1       &                 xx_vwindperiod) + 1
415          endrec   = int((modelend + startTime - diffsecs + modelstep/2)/          endrec   = int((modelend + startTime - diffsecs + modelstep/2)/
416       &                 xx_vwindperiod) + 2       &                 xx_vwindperiod) + 2
417            endif
418  # else  # else
419          startrec = 1          startrec = 1
420          endrec   = 1          endrec   = 1
# Line 535  c-- Line 575  c--
575          call cal_TimePassed( xx_sststartdate, modelstartdate,          call cal_TimePassed( xx_sststartdate, modelstartdate,
576       &                       difftime, mythid )       &                       difftime, mythid )
577          call cal_ToSeconds ( difftime, diffsecs, mythid )          call cal_ToSeconds ( difftime, diffsecs, mythid )
578            if ( xx_sstperiod .EQ. 0 ) then
579            startrec=1
580            endrec=12
581            else
582          startrec = int((modelstart + startTime - diffsecs)/          startrec = int((modelstart + startTime - diffsecs)/
583       &                 xx_sstperiod) + 1       &                 xx_sstperiod) + 1
584          endrec   = int((modelend + startTime - diffsecs + modelstep/2)/          endrec   = int((modelend + startTime - diffsecs + modelstep/2)/
585       &                 xx_sstperiod) + 2       &                 xx_sstperiod) + 2
586            endif
587  # else  # else
588          startrec = 1          startrec = 1
589          endrec   = 1          endrec   = 1
# Line 566  c-- Line 611  c--
611          call cal_TimePassed( xx_sssstartdate, modelstartdate,          call cal_TimePassed( xx_sssstartdate, modelstartdate,
612       &                       difftime, mythid )       &                       difftime, mythid )
613          call cal_ToSeconds ( difftime, diffsecs, mythid )          call cal_ToSeconds ( difftime, diffsecs, mythid )
614            if ( xx_sssperiod .EQ. 0 ) then
615            startrec=1
616            endrec=12
617            else
618          startrec = int((modelstart + startTime - diffsecs)/          startrec = int((modelstart + startTime - diffsecs)/
619       &                 xx_sssperiod) + 1       &                 xx_sssperiod) + 1
620          endrec   = int((modelend + startTime - diffsecs + modelstep/2)/          endrec   = int((modelend + startTime - diffsecs + modelstep/2)/
621       &                 xx_sssperiod) + 2       &                 xx_sssperiod) + 2
622            endif
623  # else  # else
624          startrec = 1          startrec = 1
625          endrec   = 1          endrec   = 1
# Line 692  c--   Atmos. precipitation Line 742  c--   Atmos. precipitation
742          call cal_TimePassed( xx_precipstartdate, modelstartdate,          call cal_TimePassed( xx_precipstartdate, modelstartdate,
743       &                       difftime, mythid )       &                       difftime, mythid )
744          call cal_ToSeconds ( difftime, diffsecs, mythid )          call cal_ToSeconds ( difftime, diffsecs, mythid )
745            if ( xx_precipperiod .EQ. 0 ) then
746            startrec=1
747            endrec=12
748            else
749          startrec = int((modelstart + startTime - diffsecs)/          startrec = int((modelstart + startTime - diffsecs)/
750       &                 xx_precipperiod) + 1       &                 xx_precipperiod) + 1
751          endrec   = int((modelend + startTime - diffsecs + modelstep/2)/          endrec   = int((modelend + startTime - diffsecs + modelstep/2)/
752       &                 xx_precipperiod) + 2       &                 xx_precipperiod) + 2
753            endif
754  # else  # else
755          startrec = 1          startrec = 1
756          endrec   = 1          endrec   = 1
# Line 718  c--   Atmos. swflux Line 773  c--   Atmos. swflux
773          call cal_TimePassed( xx_swfluxstartdate, modelstartdate,          call cal_TimePassed( xx_swfluxstartdate, modelstartdate,
774       &                       difftime, mythid )       &                       difftime, mythid )
775          call cal_ToSeconds ( difftime, diffsecs, mythid )          call cal_ToSeconds ( difftime, diffsecs, mythid )
776            if ( xx_swfluxperiod .EQ. 0 ) then
777            startrec=1
778            endrec=12
779            else
780          startrec = int((modelstart + startTime - diffsecs)/          startrec = int((modelstart + startTime - diffsecs)/
781       &                 xx_swfluxperiod) + 1       &                 xx_swfluxperiod) + 1
782          endrec   = int((modelend + startTime - diffsecs + modelstep/2)/          endrec   = int((modelend + startTime - diffsecs + modelstep/2)/
783       &                 xx_swfluxperiod) + 2       &                 xx_swfluxperiod) + 2
784            endif
785  # else  # else
786          startrec = 1          startrec = 1
787          endrec   = 1          endrec   = 1
# Line 744  c--   Atmos. swdown Line 804  c--   Atmos. swdown
804          call cal_TimePassed( xx_swdownstartdate, modelstartdate,          call cal_TimePassed( xx_swdownstartdate, modelstartdate,
805       &                       difftime, mythid )       &                       difftime, mythid )
806          call cal_ToSeconds ( difftime, diffsecs, mythid )          call cal_ToSeconds ( difftime, diffsecs, mythid )
807            if ( xx_swdownperiod .EQ. 0 ) then
808            startrec=1
809            endrec=12
810            else
811          startrec = int((modelstart + startTime - diffsecs)/          startrec = int((modelstart + startTime - diffsecs)/
812       &                 xx_swdownperiod) + 1       &                 xx_swdownperiod) + 1
813          endrec   = int((modelend + startTime - diffsecs + modelstep/2)/          endrec   = int((modelend + startTime - diffsecs + modelstep/2)/
814       &                 xx_swdownperiod) + 2       &                 xx_swdownperiod) + 2
815            endif
816  # else  # else
817          startrec = 1          startrec = 1
818          endrec   = 1          endrec   = 1
Legend:
Removed from v.1.21  
changed lines
  Added in v.1.22
  ViewVC Help
Powered by ViewVC 1.1.22