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

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

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


Revision 1.3 - (hide annotations) (download)
Sat Sep 17 03:17:06 2005 UTC (18 years, 8 months ago) by edhill
Branch: MAIN
CVS Tags: checkpoint57t_post, checkpoint57s_post, checkpoint57v_post, checkpoint57w_post, checkpint57u_post
Changes since 1.2: +15 -3 lines
 o fix mnc checkpoint writing problem reported by Baylor -- now works
   correctly with all the MLAdjust inputs

1 edhill 1.3 C $Header: /u/gcmpack/MITgcm/pkg/thsice/thsice_write_pickup.F,v 1.2 2005/06/24 04:36:54 edhill Exp $
2 jmc 1.1 C $Name: $
3    
4     #include "THSICE_OPTIONS.h"
5     CBOP
6     C !ROUTINE: THSICE_WRITE_PICKUP
7     C !INTERFACE:
8     SUBROUTINE THSICE_WRITE_PICKUP(
9     I prec, lgf, permCheckPoint,
10     I myIter, myThid )
11    
12     C !DESCRIPTION: \bv
13     C *==========================================================*
14     C | S/R THSICE_WRITE_PICKUP
15     C | o Write thsice pickup file
16     C *==========================================================*
17     C \ev
18    
19     C !USES:
20     IMPLICIT NONE
21    
22     C === Global variables ===
23     #include "SIZE.h"
24     #include "EEPARAMS.h"
25     #include "PARAMS.h"
26     #include "THSICE_PARAMS.h"
27     #include "THSICE_VARS.h"
28    
29     C == Common blocks ==
30     COMMON /PCKP_GBLFLS/ globalFile
31     LOGICAL globalFile
32    
33     C !INPUT/OUTPUT PARAMETERS:
34     C == Routine arguments ==
35     INTEGER prec
36     LOGICAL lgf
37     LOGICAL permCheckPoint
38     INTEGER myIter
39     INTEGER myThid
40     CEOP
41    
42     #ifdef ALLOW_THSICE
43    
44     C == Local variables ==
45     CHARACTER*(MAX_LEN_FNAM) fn
46    
47 edhill 1.2 IF ( permCheckPoint ) THEN
48     WRITE(fn,'(A,I10.10)') 'pickup_ic.',myIter
49     ELSE
50     WRITE(fn,'(A,A)') 'pickup_ic.',checkPtSuff(nCheckLev)
51     ENDIF
52 jmc 1.1
53 edhill 1.2 IF ( thSIce_pickup_write_mdsio ) THEN
54 jmc 1.1 CALL MDSWRITEFIELD(fn,prec,lgf,'RL',1,iceMask, 1,
55     & myIter,myThid)
56     CALL MDSWRITEFIELD(fn,prec,lgf,'RL',1,iceHeight, 2,
57     & myIter,myThid)
58     CALL MDSWRITEFIELD(fn,prec,lgf,'RL',1,snowHeight,3,
59     & myIter,myThid)
60     CALL MDSWRITEFIELD(fn,prec,lgf,'RL',1,Tsrf, 4, myIter,myThid)
61     CALL MDSWRITEFIELD(fn,prec,lgf,'RL',1,Tice1, 5, myIter,myThid)
62     CALL MDSWRITEFIELD(fn,prec,lgf,'RL',1,Tice2, 6, myIter,myThid)
63     CALL MDSWRITEFIELD(fn,prec,lgf,'RL',1,Qice1, 7, myIter,myThid)
64     CALL MDSWRITEFIELD(fn,prec,lgf,'RL',1,Qice2, 8, myIter,myThid)
65     CALL MDSWRITEFIELD(fn,prec,lgf,'RL',1,snowAge,9, myIter,myThid)
66    
67 edhill 1.2 IF ( stepFwd_oceMxL ) THEN
68 jmc 1.1 CALL MDSWRITEFIELD(fn,prec,lgf,'RL',1,tOceMxL,10,myIter,myThid)
69     CALL MDSWRITEFIELD(fn,prec,lgf,'RL',1,sOceMxL,11,myIter,myThid)
70 edhill 1.2 ENDIF
71     ENDIF
72    
73     #ifdef ALLOW_MNC
74     IF ( thSIce_pickup_write_mnc ) THEN
75 edhill 1.3 IF ( permCheckPoint ) THEN
76     WRITE(fn,'(A)') 'pickup_ic'
77     ELSE
78     WRITE(fn,'(A,A)') 'pickup_ic.',checkPtSuff(nCheckLev)
79     ENDIF
80     C First ***define*** the file group name
81     CALL MNC_CW_SET_UDIM(fn, 0, myThid)
82     IF ( permCheckPoint ) THEN
83     CALL MNC_CW_SET_CITER(fn, 3, 3, myIter, 0, myThid)
84     ELSE
85     CALL MNC_CW_SET_CITER(fn, 2, -1, -1, -1, myThid)
86     ENDIF
87     C Then set the actual unlimited dimension
88     CALL MNC_CW_SET_UDIM(fn, 1, myThid)
89 edhill 1.2 CALL MNC_CW_I_W_S('I',fn,0,0,'iter', myIter, myThid)
90     C CALL MNC_CW_RL_W_S('D',fn,0,0,'T', myTime, myThid)
91     CALL MNC_CW_RL_W('D',fn,0,0,'iceMask', iceMask, myThid)
92     CALL MNC_CW_RL_W('D',fn,0,0,'iceHeight', iceHeight, myThid)
93     CALL MNC_CW_RL_W('D',fn,0,0,'snowHeight',snowHeight,myThid)
94     CALL MNC_CW_RL_W('D',fn,0,0,'Tsrf', Tsrf, myThid)
95     CALL MNC_CW_RL_W('D',fn,0,0,'Tice1', Tice1, myThid)
96     CALL MNC_CW_RL_W('D',fn,0,0,'Tice2', Tice1, myThid)
97     CALL MNC_CW_RL_W('D',fn,0,0,'Qice1', Qice1, myThid)
98     CALL MNC_CW_RL_W('D',fn,0,0,'Qice2', Qice2, myThid)
99     CALL MNC_CW_RL_W('D',fn,0,0,'snowAge', snowAge, myThid)
100     IF ( stepFwd_oceMxL ) THEN
101     CALL MNC_CW_RL_W('D',fn,0,0,'tOceMxL',tOceMxL,myThid)
102     CALL MNC_CW_RL_W('D',fn,0,0,'sOceMxL',sOceMxL,myThid)
103     ENDIF
104     ENDIF
105     #endif /* ALLOW_MNC */
106 jmc 1.1
107     #endif /* ALLOW_THSICE */
108     RETURN
109     END
110 edhill 1.2
111     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|

  ViewVC Help
Powered by ViewVC 1.1.22