27 |
|
|
28 |
C !INPUT/OUTPUT PARAMETERS: |
C !INPUT/OUTPUT PARAMETERS: |
29 |
C == Routine arguments == |
C == Routine arguments == |
30 |
C myTime - Current time of simulation ( s ) |
C myTime :: Current time of simulation ( s ) |
31 |
C myIter - Iteration number |
C myIter :: Iteration number |
32 |
C myThid - Number of this instance of INI_FORCING |
C myThid :: my Thread Id. number |
33 |
_RL myTime |
_RL myTime |
34 |
INTEGER myIter |
INTEGER myIter |
35 |
INTEGER myThid |
INTEGER myThid |
72 |
|
|
73 |
IF ( MASTER_CPU_IO(myThid) ) THEN |
IF ( MASTER_CPU_IO(myThid) ) THEN |
74 |
C-- only the master thread is allowed to switch On/Off mon_write_stdout |
C-- only the master thread is allowed to switch On/Off mon_write_stdout |
75 |
C & mon_write_mnc (since it's the only thread that uses those flags): |
C & mon_write_mnc (since it is the only thread that uses those flags): |
76 |
|
|
77 |
IF ( thSIce_mon_stdio ) THEN |
IF ( thSIce_mon_stdio ) THEN |
78 |
mon_write_stdout = .TRUE. |
mon_write_stdout = .TRUE. |
128 |
C-- Ice area and Ice thickness : |
C-- Ice area and Ice thickness : |
129 |
CALL MON_STATS_LATBND_RL( |
CALL MON_STATS_LATBND_RL( |
130 |
I 1, 1, 1, 2, yBand, |
I 1, 1, 1, 2, yBand, |
131 |
I iceHeight, locMask, maskH, rA, yC, locDr, |
I iceHeight, locMask, maskInC, rA, yC, locDr, |
132 |
O theMin, theMax, theMean, theVar, theVol, |
O theMin, theMax, theMean, theVar, theVol, |
133 |
I myThid ) |
I myThid ) |
134 |
theVolG= theVol(1)+theVol(2) |
theVolG= theVol(1)+theVol(2) |
150 |
C-- Snow thickness : |
C-- Snow thickness : |
151 |
CALL MON_STATS_LATBND_RL( |
CALL MON_STATS_LATBND_RL( |
152 |
I 1, 1, 1, 2, yBand, |
I 1, 1, 1, 2, yBand, |
153 |
I snowHeight, locMask, maskH, rA, yC, locDr, |
I snowHeight, locMask, maskInC, rA, yC, locDr, |
154 |
O theMin, theMax, theMean, theVar, theVol, |
O theMin, theMax, theMean, theVar, theVol, |
155 |
I myThid ) |
I myThid ) |
156 |
theVolG= theVol(1)+theVol(2) |
theVolG= theVol(1)+theVol(2) |
169 |
C-- Surface Temp. : |
C-- Surface Temp. : |
170 |
CALL MON_STATS_LATBND_RL( |
CALL MON_STATS_LATBND_RL( |
171 |
I 1, 1, 1, 2, yBand, |
I 1, 1, 1, 2, yBand, |
172 |
I Tsrf, locMask, maskH, rA, yC, locDr, |
I Tsrf, locMask, maskInC, rA, yC, locDr, |
173 |
O theMin, theMax, theMean, theVar, theVol, |
O theMin, theMax, theMean, theVar, theVol, |
174 |
I myThid ) |
I myThid ) |
175 |
theVolG= theVol(1)+theVol(2) |
theVolG= theVol(1)+theVol(2) |
201 |
C-- 1rst level (volume-mean) Temp. : |
C-- 1rst level (volume-mean) Temp. : |
202 |
CALL MON_STATS_LATBND_RL( |
CALL MON_STATS_LATBND_RL( |
203 |
I 1, 1, 1, 2, yBand, |
I 1, 1, 1, 2, yBand, |
204 |
I Tice1, locMask, maskH, rA, yC, locDr, |
I Tice1, locMask, maskInC, rA, yC, locDr, |
205 |
O theMin, theMax, theMean, theVar, theVol, |
O theMin, theMax, theMean, theVar, theVol, |
206 |
I myThid ) |
I myThid ) |
207 |
theVolG = theVol(1)+theVol(2) |
theVolG = theVol(1)+theVol(2) |
226 |
C-- 2nd level (volume-mean) Temp. : |
C-- 2nd level (volume-mean) Temp. : |
227 |
CALL MON_STATS_LATBND_RL( |
CALL MON_STATS_LATBND_RL( |
228 |
I 1, 1, 1, 2, yBand, |
I 1, 1, 1, 2, yBand, |
229 |
I Tice2, locMask, maskH, rA, yC, locDr, |
I Tice2, locMask, maskInC, rA, yC, locDr, |
230 |
O theMin, theMax, theMean, theVar, theVol, |
O theMin, theMax, theMean, theVar, theVol, |
231 |
I myThid ) |
I myThid ) |
232 |
theMeanG= theMean(1)*theVol(1)+theMean(2)*theVol(2) |
theMeanG= theMean(1)*theVol(1)+theMean(2)*theVol(2) |
244 |
CALL MON_OUT_RL(mon_var, theMax(2), mon_sufx(2), myThid) |
CALL MON_OUT_RL(mon_var, theMax(2), mon_sufx(2), myThid) |
245 |
|
|
246 |
C-- Total Energy : |
C-- Total Energy : |
247 |
CALL MON_STATS_RL( |
CALL MON_CALC_STATS_RL( |
248 |
I 1, Qice1, locMask, maskH, rA, locDr, |
I 1, Qice1, locMask, maskInC, rA, locDr, |
249 |
O theMin0,theMax0,theMean1,theSD,theDel2,theVolG, |
O theMin0,theMax0,theMean1,theSD,theDel2,theVolG, |
250 |
I myThid ) |
I myThid ) |
251 |
CALL MON_STATS_RL( |
CALL MON_CALC_STATS_RL( |
252 |
I 1, Qice2, locMask, maskH, rA, locDr, |
I 1, Qice2, locMask, maskInC, rA, locDr, |
253 |
O theMin0,theMax0,theMean2,theSD,theDel2,theVolG, |
O theMin0,theMax0,theMean2,theSD,theDel2,theVolG, |
254 |
I myThid ) |
I myThid ) |
255 |
theEnergy = theEnergy -rhoi*(theMean1+theMean2)*theVolG/2 |
theEnergy = theEnergy -rhoi*(theMean1+theMean2)*theVolG/2 |
258 |
|
|
259 |
IF ( MASTER_CPU_IO(myThid) ) THEN |
IF ( MASTER_CPU_IO(myThid) ) THEN |
260 |
C-- only the master thread is allowed to switch On/Off mon_write_stdout |
C-- only the master thread is allowed to switch On/Off mon_write_stdout |
261 |
C & mon_write_mnc (since it's the only thread that uses those flags): |
C & mon_write_mnc (since it is the only thread that uses those flags): |
262 |
|
|
263 |
IF (mon_write_stdout) THEN |
IF (mon_write_stdout) THEN |
264 |
WRITE(msgBuf,'(2A)') '// ==========================', |
WRITE(msgBuf,'(2A)') '// ==========================', |