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

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

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


Revision 1.3 - (show 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 C $Header: /u/gcmpack/MITgcm/pkg/thsice/thsice_write_pickup.F,v 1.2 2005/06/24 04:36:54 edhill Exp $
2 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 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
53 IF ( thSIce_pickup_write_mdsio ) THEN
54 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 IF ( stepFwd_oceMxL ) THEN
68 CALL MDSWRITEFIELD(fn,prec,lgf,'RL',1,tOceMxL,10,myIter,myThid)
69 CALL MDSWRITEFIELD(fn,prec,lgf,'RL',1,sOceMxL,11,myIter,myThid)
70 ENDIF
71 ENDIF
72
73 #ifdef ALLOW_MNC
74 IF ( thSIce_pickup_write_mnc ) THEN
75 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 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
107 #endif /* ALLOW_THSICE */
108 RETURN
109 END
110
111 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|

  ViewVC Help
Powered by ViewVC 1.1.22