/[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.2 - (hide 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 edhill 1.2 C $Header: /u/gcmpack/MITgcm/pkg/thsice/thsice_read_pickup.F,v 1.1 2004/04/07 23:40:34 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     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 edhill 1.2
40     IF ( .NOT. thSIce_pickup_write_mnc ) THEN
41 jmc 1.1
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 edhill 1.2 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 jmc 1.1
80     #endif /* ALLOW_THSICE */
81     RETURN
82     END
  ViewVC Help
Powered by ViewVC 1.1.22