/[MITgcm]/MITgcm/pkg/seaice/seaice_mnc_init.F
ViewVC logotype

Annotation of /MITgcm/pkg/seaice/seaice_mnc_init.F

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


Revision 1.3 - (hide annotations) (download)
Thu Aug 11 02:50:12 2005 UTC (18 years, 9 months ago) by edhill
Branch: MAIN
CVS Tags: checkpoint57t_post, checkpoint57v_post, checkpoint57s_post, checkpoint57y_post, checkpoint57y_pre, checkpoint57r_post, checkpoint58, checkpoint57x_post, checkpoint57w_post, checkpint57u_post, checkpoint58a_post, checkpoint57q_post, checkpoint57z_post, checkpoint58b_post
Changes since 1.2: +21 -1 lines
 o fix an uninitialized variable problem introduced when seaice was
   MNC-ified a couple of months ago -- tested and works with both g77
   and ifort v8.1

1 edhill 1.3 C $Header: /u/gcmpack/MITgcm/pkg/seaice/seaice_mnc_init.F,v 1.2 2005/08/10 19:18:52 edhill Exp $
2 edhill 1.1 C $Name: $
3    
4     #include "SEAICE_OPTIONS.h"
5    
6     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
7     CBOP
8     C !ROUTINE: SEAICE_MNC_INIT
9    
10     C !INTERFACE:
11     SUBROUTINE SEAICE_MNC_INIT( myThid )
12    
13     C !DESCRIPTION:
14     C Initialize SEAICE variables within MNC
15    
16     C !USES:
17     IMPLICIT NONE
18     #include "SIZE.h"
19     #include "EEPARAMS.h"
20     #include "PARAMS.h"
21    
22     C !INPUT PARAMETERS:
23     INTEGER myThid
24     CEOP
25    
26     #ifdef ALLOW_SEAICE
27     #ifdef ALLOW_MNC
28     IF (useMNC) THEN
29    
30 edhill 1.2 CALL MNC_CW_ADD_VNAME('si_UWIND','U_xy_Hn__-__t',3,4,myThid)
31     CALL MNC_CW_ADD_VATTR_TEXT('si_UWIND','units','m/s', myThid)
32     CALL MNC_CW_ADD_VATTR_TEXT('si_UWIND','description',
33 edhill 1.1 & 'Surface (10m) zonal wind velocity', myThid)
34    
35 edhill 1.2 CALL MNC_CW_ADD_VNAME('si_VWIND','V_xy_Hn__-__t',3,4,myThid)
36     CALL MNC_CW_ADD_VATTR_TEXT('si_VWIND','units','m/s', myThid)
37     CALL MNC_CW_ADD_VATTR_TEXT('si_VWIND','description',
38 edhill 1.1 & 'Surface (10m) meridional wind velocity', myThid)
39    
40 edhill 1.2 CALL MNC_CW_ADD_VNAME('si_EVAP','Cen_xy_Hn__-__t',3,4,myThid)
41     CALL MNC_CW_ADD_VATTR_TEXT('si_EVAP','units','m/s', myThid)
42     CALL MNC_CW_ADD_VATTR_TEXT('si_EVAP','description',
43 edhill 1.1 & 'Evaporation', myThid)
44    
45 edhill 1.2 CALL MNC_CW_ADD_VNAME('si_PRECIP','Cen_xy_Hn__-__t',3,4,myThid)
46     CALL MNC_CW_ADD_VATTR_TEXT('si_PRECIP','units','m/s', myThid)
47     CALL MNC_CW_ADD_VATTR_TEXT('si_PRECIP','description',
48 edhill 1.1 & 'Precipitation', myThid)
49    
50 edhill 1.2 CALL MNC_CW_ADD_VNAME('si_RUNOFF','Cen_xy_Hn__-__t',3,4,myThid)
51     CALL MNC_CW_ADD_VATTR_TEXT('si_RUNOFF','units','m/s', myThid)
52     CALL MNC_CW_ADD_VATTR_TEXT('si_RUNOFF','description',
53 edhill 1.1 & 'River and glacier runoff', myThid)
54    
55 edhill 1.2 CALL MNC_CW_ADD_VNAME('si_ATEMP','Cen_xy_Hn__-__t',3,4,myThid)
56     CALL MNC_CW_ADD_VATTR_TEXT('si_ATEMP','units','deg K', myThid)
57     CALL MNC_CW_ADD_VATTR_TEXT('si_ATEMP','description',
58 edhill 1.1 & 'Surface (2m) air temperature', myThid)
59    
60 edhill 1.2 CALL MNC_CW_ADD_VNAME('si_AQH','Cen_xy_Hn__-__t',3,4,myThid)
61     CALL MNC_CW_ADD_VATTR_TEXT('si_AQH','units','kg/kg', myThid)
62     CALL MNC_CW_ADD_VATTR_TEXT('si_AQH','description',
63 edhill 1.1 & 'Surface (2m) specific humidity', myThid)
64    
65 edhill 1.2 CALL MNC_CW_ADD_VNAME('si_SWDOWN','Cen_xy_Hn__-__t',3,4,myThid)
66     CALL MNC_CW_ADD_VATTR_TEXT('si_SWDOWN','units','W/m^2', myThid)
67     CALL MNC_CW_ADD_VATTR_TEXT('si_SWDOWN','description',
68 edhill 1.1 & 'Downward shortwave radiation', myThid)
69    
70 edhill 1.2 CALL MNC_CW_ADD_VNAME('si_LWDOWN','Cen_xy_Hn__-__t',3,4,myThid)
71     CALL MNC_CW_ADD_VATTR_TEXT('si_LWDOWN','units','W/m^2', myThid)
72     CALL MNC_CW_ADD_VATTR_TEXT('si_LWDOWN','description',
73 edhill 1.1 & 'Downward longwave radiation', myThid)
74    
75 edhill 1.3 CALL MNC_CW_ADD_VNAME('si_UICE','U_xy_Hn__-__t',3,4,myThid)
76     CALL MNC_CW_ADD_VATTR_TEXT('si_UICE','units','m/s', myThid)
77     CALL MNC_CW_ADD_VATTR_TEXT('si_UICE','description',
78     & 'Ice U velocity', myThid)
79    
80     CALL MNC_CW_ADD_VNAME('si_VICE','V_xy_Hn__-__t',3,4,myThid)
81     CALL MNC_CW_ADD_VATTR_TEXT('si_VICE','units','m/s', myThid)
82     CALL MNC_CW_ADD_VATTR_TEXT('si_VICE','description',
83     & 'Ice V velocity', myThid)
84    
85     CALL MNC_CW_ADD_VNAME('si_HEFF','Cen_xy_Hn__-__t',3,4,myThid)
86     CALL MNC_CW_ADD_VATTR_TEXT('si_HEFF','units','m', myThid)
87     CALL MNC_CW_ADD_VATTR_TEXT('si_HEFF','description',
88     & 'Effective ice thickness',myThid)
89    
90     CALL MNC_CW_ADD_VNAME('si_AREA','Cen_xy_Hn__-__t',3,4,myThid)
91     CALL MNC_CW_ADD_VATTR_TEXT('si_AREA','units','m^2/m^2',myThid)
92     CALL MNC_CW_ADD_VATTR_TEXT('si_AREA','description',
93     & 'Fractional ice-covered area',myThid)
94    
95 edhill 1.1 C Time averages
96 edhill 1.2 CALL MNC_CW_ADD_VNAME('si_UICEtave','U_xy_Hn__-__t',3,4,myThid)
97     CALL MNC_CW_ADD_VATTR_TEXT('si_UICEtave','units','m/s', myThid)
98     CALL MNC_CW_ADD_VATTR_TEXT('si_UICEtave','description',
99 edhill 1.1 & 'Time-averaged ', myThid)
100    
101 edhill 1.2 CALL MNC_CW_ADD_VNAME('si_VICEtave','V_xy_Hn__-__t',3,4,myThid)
102     CALL MNC_CW_ADD_VATTR_TEXT('si_VICEtave','units','m/s', myThid)
103     CALL MNC_CW_ADD_VATTR_TEXT('si_VICEtave','description',
104 edhill 1.1 & 'Time-averaged ', myThid)
105    
106 edhill 1.2 CALL MNC_CW_ADD_VNAME('si_FUtave','U_xy_Hn__-__t',3,4,myThid)
107     CALL MNC_CW_ADD_VATTR_TEXT('si_FUtave','units','N/m^2', myThid)
108     CALL MNC_CW_ADD_VATTR_TEXT('si_FUtave','description',
109 edhill 1.1 & 'Time-averaged zonal surface wind stress', myThid)
110    
111 edhill 1.2 CALL MNC_CW_ADD_VNAME('si_FVtave','V_xy_Hn__-__t',3,4,myThid)
112     CALL MNC_CW_ADD_VATTR_TEXT('si_FVtave','units','N/m^2',myThid)
113     CALL MNC_CW_ADD_VATTR_TEXT('si_FVtave','description',
114 edhill 1.1 & 'Time-averaged Meridional surface wind stress', myThid)
115    
116 edhill 1.2 CALL MNC_CW_ADD_VNAME(
117     & 'si_EmPmRtave','Cen_xy_Hn__-__t',3,4,myThid)
118     CALL MNC_CW_ADD_VATTR_TEXT(
119     & 'si_EmPmRtave','units','m/s',myThid)
120     CALL MNC_CW_ADD_VATTR_TEXT('si_EmPmRtave','description',
121 edhill 1.1 & 'Time-averaged net upward freshwater flux', myThid)
122    
123 edhill 1.2 CALL MNC_CW_ADD_VNAME(
124     & 'si_QNETtave','Cen_xy_Hn__-__t',3,4,myThid)
125     CALL MNC_CW_ADD_VATTR_TEXT(
126     & 'si_QNETtave','units','W/m^2', myThid)
127     CALL MNC_CW_ADD_VATTR_TEXT('si_QNETtave','description',
128 edhill 1.1 & 'Time-averaged net upward surface heat flux', myThid)
129    
130 edhill 1.2 CALL MNC_CW_ADD_VNAME(
131     & 'si_QSWtave','Cen_xy_Hn__-__t',3,4,myThid)
132     CALL MNC_CW_ADD_VATTR_TEXT(
133     & 'si_QSWtave','units','W/m^2',myThid)
134     CALL MNC_CW_ADD_VATTR_TEXT('si_QSWtave','description',
135     & 'Time-averaged net upward shortwave radiation',myThid)
136    
137     CALL MNC_CW_ADD_VNAME(
138     & 'si_HEFFtave','Cen_xy_Hn__-__t',3,4,myThid)
139     CALL MNC_CW_ADD_VATTR_TEXT(
140     & 'si_HEFFtave','units','m', myThid)
141     CALL MNC_CW_ADD_VATTR_TEXT('si_HEFFtave','description',
142     & 'Time-averaged effective ice thickness',myThid)
143    
144     CALL MNC_CW_ADD_VNAME(
145     & 'si_AREAtave','Cen_xy_Hn__-__t',3,4,myThid)
146     CALL MNC_CW_ADD_VATTR_TEXT(
147     & 'si_AREAtave','units','m^2/m^2',myThid)
148     CALL MNC_CW_ADD_VATTR_TEXT('si_AREAtave','description',
149     & 'Time-averaged fractional ice-covered area',myThid)
150 edhill 1.1
151     ENDIF
152     #endif /* ALLOW_MNC */
153     #endif /* ALLOW_SEAICE */
154    
155     RETURN
156     END
157    

  ViewVC Help
Powered by ViewVC 1.1.22