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

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

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

Revision 1.2 - (show annotations) (download)
Fri Jun 24 04:36:54 2005 UTC (18 years, 10 months ago) by edhill
Branch: MAIN
CVS Tags: checkpoint57o_post, checkpoint57m_post, checkpoint57k_post, checkpoint57r_post, checkpoint57n_post, checkpoint57p_post, checkpoint57q_post, checkpoint57j_post, checkpoint57l_post
Changes since 1.1: +29 -5 lines 
 o mnc-ify the thsice package as requested by Daniel Enderton
   - the monitor--MNC output needs work
   - many attributes need to be added (most are currently blank)
   - does not break testreport (at least on IA32)
1 C $Header: /u/gcmpack/MITgcm/pkg/thsice/thsice_read_pickup.F,v 1.1 2004/04/07 23:40:34 jmc Exp $
2 C $Name: $
3
4 #include "THSICE_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: THSICE_READ_PICKUP
8 C !INTERFACE:
9 SUBROUTINE THSICE_READ_PICKUP( prec, myIter, myThid )
10
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 prec
31 INTEGER myIter
32 INTEGER myThid
33 CEOP
34
35 #ifdef ALLOW_THSICE
36
37 C == Local variables ==
38 CHARACTER*(MAX_LEN_FNAM) fn
39
40 IF ( .NOT. thSIce_pickup_write_mnc ) THEN
41
42 WRITE(fn,'(A,I10.10)') 'pickup_ic.',myIter
43 CALL MDSREADFIELD(fn,prec,'RL',1,iceMask, 1,myThid)
44 CALL MDSREADFIELD(fn,prec,'RL',1,iceHeight, 2,myThid)
45 CALL MDSREADFIELD(fn,prec,'RL',1,snowHeight,3,myThid)
46 CALL MDSREADFIELD(fn,prec,'RL',1,Tsrf, 4,myThid)
47 CALL MDSREADFIELD(fn,prec,'RL',1,Tice1, 5,myThid)
48 CALL MDSREADFIELD(fn,prec,'RL',1,Tice2, 6,myThid)
49 CALL MDSREADFIELD(fn,prec,'RL',1,Qice1, 7,myThid)
50 CALL MDSREADFIELD(fn,prec,'RL',1,Qice2, 8,myThid)
51 CALL MDSREADFIELD(fn,prec,'RL',1,snowAge, 9,myThid)
52
53 IF ( stepFwd_oceMxL .AND. startIceModel.EQ.0 ) THEN
54 CALL MDSREADFIELD(fn,prec,'RL',1,tOceMxL, 10,myThid)
55 CALL MDSREADFIELD(fn,prec,'RL',1,sOceMxL, 11,myThid)
56 ENDIF
57
58 ENDIF
59
60 #ifdef ALLOW_MNC
61 IF ( thSIce_pickup_write_mnc ) THEN
62 CALL MNC_FILE_CLOSE_ALL_MATCHING(fn, myThid)
63 CALL MNC_CW_SET_UDIM(fn, 1, myThid)
64 CALL MNC_CW_RL_R('D',fn,0,0,'iceMask', iceMask, myThid)
65 CALL MNC_CW_RL_R('D',fn,0,0,'iceHeight', iceHeight, myThid)
66 CALL MNC_CW_RL_R('D',fn,0,0,'snowHeight',snowHeight,myThid)
67 CALL MNC_CW_RL_R('D',fn,0,0,'Tsrf', Tsrf, myThid)
68 CALL MNC_CW_RL_R('D',fn,0,0,'Tice1', Tice1, myThid)
69 CALL MNC_CW_RL_R('D',fn,0,0,'Tice2', Tice1, myThid)
70 CALL MNC_CW_RL_R('D',fn,0,0,'Qice1', Qice1, myThid)
71 CALL MNC_CW_RL_R('D',fn,0,0,'Qice2', Qice2, myThid)
72 CALL MNC_CW_RL_R('D',fn,0,0,'snowAge', snowAge, myThid)
73 IF ( stepFwd_oceMxL ) THEN
74 CALL MNC_CW_RL_R('D',fn,0,0,'tOceMxL',tOceMxL,myThid)
75 CALL MNC_CW_RL_R('D',fn,0,0,'sOceMxL',sOceMxL,myThid)
76 ENDIF
77 ENDIF
78 #endif /* ALLOW_MNC */
79
80 #endif /* ALLOW_THSICE */
81 RETURN
82 END
  ViewVC Help
Powered by ViewVC 1.1.22