/[MITgcm]/MITgcm/pkg/thsice/thsice_read_pickup.F
ViewVC logotype

Annotation of /MITgcm/pkg/thsice/thsice_read_pickup.F

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


Revision 1.6 - (hide annotations) (download)
Fri Mar 24 23:51:14 2017 UTC (7 years, 1 month ago) by jmc
Branch: MAIN
CVS Tags: checkpoint66g, checkpoint66f, checkpoint66o, checkpoint66n, checkpoint66m, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, HEAD
Changes since 1.5: +14 -7 lines
use new S/R RW_GET_SUFFIX to get file suffix (according to "rwSuffixType")

1 jmc 1.6 C $Header: /u/gcmpack/MITgcm/pkg/thsice/thsice_read_pickup.F,v 1.5 2007/10/19 03:22:25 jmc Exp $
2 jmc 1.1 C $Name: $
3    
4     #include "THSICE_OPTIONS.h"
5    
6     CBOP
7     C !ROUTINE: THSICE_READ_PICKUP
8     C !INTERFACE:
9 jmc 1.6 SUBROUTINE THSICE_READ_PICKUP( myIter, myThid )
10 jmc 1.1
11     C !DESCRIPTION: \bv
12     C *==========================================================*
13     C | S/R THSICE_READ_PICKUP
14     C | o Read thsice pickup file
15     C *==========================================================*
16     C \ev
17    
18     C !USES:
19     IMPLICIT NONE
20    
21     C === Global variables ===
22     #include "SIZE.h"
23     #include "EEPARAMS.h"
24     #include "PARAMS.h"
25     #include "THSICE_PARAMS.h"
26     #include "THSICE_VARS.h"
27    
28     C !INPUT/OUTPUT PARAMETERS:
29     C == Routine arguments ==
30     INTEGER myIter
31     INTEGER myThid
32     CEOP
33    
34     #ifdef ALLOW_THSICE
35    
36     C == Local variables ==
37 jmc 1.6 INTEGER prec
38     CHARACTER*(10) suff
39 jmc 1.1 CHARACTER*(MAX_LEN_FNAM) fn
40 jmc 1.6
41 edhill 1.2 IF ( .NOT. thSIce_pickup_write_mnc ) THEN
42 jmc 1.1
43 jmc 1.6 prec = precFloat64
44     IF ( pickupSuff.EQ.' ' ) THEN
45     IF ( rwSuffixType.EQ.0 ) THEN
46     WRITE(fn,'(A,I10.10)') 'pickup_ic.', myIter
47     ELSE
48     CALL RW_GET_SUFFIX( suff, startTime, myIter, myThid )
49     WRITE(fn,'(A,A)') 'pickup_ic.', suff
50     ENDIF
51 jmc 1.5 ELSE
52 jmc 1.6 WRITE(fn,'(A,A10)') 'pickup_ic.', pickupSuff
53 jmc 1.5 ENDIF
54     CALL READ_REC_3D_RL( fn, prec, 1, iceMask, 1, myIter,myThid )
55     CALL READ_REC_3D_RL( fn, prec, 1, iceHeight,2, myIter,myThid )
56     CALL READ_REC_3D_RL( fn, prec, 1,snowHeight,3, myIter,myThid )
57     CALL READ_REC_3D_RL( fn, prec, 1, Tsrf, 4, myIter,myThid )
58     CALL READ_REC_3D_RL( fn, prec, 1, Tice1, 5, myIter,myThid )
59     CALL READ_REC_3D_RL( fn, prec, 1, Tice2, 6, myIter,myThid )
60     CALL READ_REC_3D_RL( fn, prec, 1, Qice1, 7, myIter,myThid )
61     CALL READ_REC_3D_RL( fn, prec, 1, Qice2, 8, myIter,myThid )
62     CALL READ_REC_3D_RL( fn, prec, 1, snowAge, 9, myIter,myThid )
63 jmc 1.1
64 edhill 1.2 IF ( stepFwd_oceMxL .AND. startIceModel.EQ.0 ) THEN
65 jmc 1.5 CALL READ_REC_3D_RL( fn,prec,1, tOceMxL, 10, myIter,myThid )
66     CALL READ_REC_3D_RL( fn,prec,1, sOceMxL, 11, myIter,myThid )
67 edhill 1.2 ENDIF
68    
69     ENDIF
70    
71     #ifdef ALLOW_MNC
72     IF ( thSIce_pickup_write_mnc ) THEN
73 edhill 1.3 WRITE(fn,'(A)') 'pickup_ic'
74 edhill 1.2 CALL MNC_FILE_CLOSE_ALL_MATCHING(fn, myThid)
75     CALL MNC_CW_SET_UDIM(fn, 1, myThid)
76 mlosch 1.4 CALL MNC_CW_SET_CITER(fn, 3, 3, myIter, 0, myThid)
77 edhill 1.2 CALL MNC_CW_RL_R('D',fn,0,0,'iceMask', iceMask, myThid)
78     CALL MNC_CW_RL_R('D',fn,0,0,'iceHeight', iceHeight, myThid)
79     CALL MNC_CW_RL_R('D',fn,0,0,'snowHeight',snowHeight,myThid)
80     CALL MNC_CW_RL_R('D',fn,0,0,'Tsrf', Tsrf, myThid)
81     CALL MNC_CW_RL_R('D',fn,0,0,'Tice1', Tice1, myThid)
82     CALL MNC_CW_RL_R('D',fn,0,0,'Tice2', Tice1, myThid)
83     CALL MNC_CW_RL_R('D',fn,0,0,'Qice1', Qice1, myThid)
84     CALL MNC_CW_RL_R('D',fn,0,0,'Qice2', Qice2, myThid)
85     CALL MNC_CW_RL_R('D',fn,0,0,'snowAge', snowAge, myThid)
86     IF ( stepFwd_oceMxL ) THEN
87     CALL MNC_CW_RL_R('D',fn,0,0,'tOceMxL',tOceMxL,myThid)
88     CALL MNC_CW_RL_R('D',fn,0,0,'sOceMxL',sOceMxL,myThid)
89     ENDIF
90     ENDIF
91     #endif /* ALLOW_MNC */
92 jmc 1.1
93     #endif /* ALLOW_THSICE */
94     RETURN
95     END

  ViewVC Help
Powered by ViewVC 1.1.22