/[MITgcm]/MITgcm/pkg/atm2d/sum_ocn_fluxes.F
ViewVC logotype

Contents of /MITgcm/pkg/atm2d/sum_ocn_fluxes.F

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

Revision 1.3 - (show annotations) (download)
Mon Oct 8 23:48:28 2007 UTC (16 years, 7 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint64y, checkpoint64x, checkpoint64z, checkpoint64q, checkpoint64p, checkpoint64s, checkpoint64r, checkpoint64u, checkpoint64t, checkpoint64w, checkpoint64v, checkpoint64i, checkpoint64h, checkpoint64k, checkpoint64j, checkpoint64m, checkpoint64l, checkpoint64o, checkpoint64n, checkpoint64a, checkpoint64c, checkpoint64b, checkpoint64e, checkpoint64d, checkpoint64g, checkpoint64f, checkpoint63p, checkpoint63q, checkpoint63r, checkpoint63s, checkpoint63l, checkpoint63m, checkpoint63n, checkpoint63o, checkpoint63h, checkpoint63i, checkpoint63j, checkpoint63k, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint63g, checkpoint63a, checkpoint63b, checkpoint63c, checkpoint64, checkpoint65, checkpoint60, checkpoint61, checkpoint62, checkpoint63, checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint66o, checkpoint66n, checkpoint66m, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, checkpoint65z, checkpoint65x, checkpoint65y, checkpoint65r, checkpoint65s, checkpoint65p, checkpoint65q, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, checkpoint65j, checkpoint65k, checkpoint65h, checkpoint65i, checkpoint65n, checkpoint65o, checkpoint65l, checkpoint65m, checkpoint65b, checkpoint65c, checkpoint65a, checkpoint65f, checkpoint65g, checkpoint65d, checkpoint65e, checkpoint59q, checkpoint59p, checkpoint59r, checkpoint59m, checkpoint59l, checkpoint59o, checkpoint59n, checkpoint59i, checkpoint59k, checkpoint59j, checkpoint62c, checkpoint62b, checkpoint62a, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62s, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint62w, checkpoint62v, checkpoint62u, checkpoint62t, checkpoint62z, checkpoint62y, checkpoint62x, checkpoint61f, checkpoint61g, checkpoint61d, checkpoint61e, checkpoint61b, checkpoint61c, checkpoint61a, checkpoint61n, checkpoint61o, checkpoint61l, checkpoint61m, checkpoint61j, checkpoint61k, checkpoint61h, checkpoint61i, checkpoint61v, checkpoint61w, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61p, checkpoint61q, checkpoint61z, checkpoint61x, checkpoint61y, HEAD
Changes since 1.2: +14 -11 lines 
add missing cvs $Header:$ or $Name:$
1 C $Header: $
2 C $Name: $
3
4 #include "ctrparam.h"
5 #include "ATM2D_OPTIONS.h"
6
7 C !INTERFACE:
8 SUBROUTINE SUM_OCN_FLUXES( myThid )
9 C *==========================================================*
10 C | Sums the atmos-> ocn fluxes. Note the stress reduction |
11 c | (wind forcing) which occurs as icefract > windice_thres |
12 C *==========================================================*
13 IMPLICIT NONE
14
15 #include "ATMSIZE.h"
16 #include "SIZE.h"
17 #include "GRID.h"
18 #include "EEPARAMS.h"
19
20 C === Global SeaIce Variables ===
21 #include "THSICE_VARS.h"
22
23 C === Atmos/Ocean/Seaice Interface Variables ===
24 #include "ATM2D_VARS.h"
25
26 c start phasing out wind stress to ocean at this ice fraction
27 _RS windice_thres
28 PARAMETER ( windice_thres= 0.5 )
29
30 C !INPUT/OUTPUT PARAMETERS:
31 C === Routine arguments ===
32 C myThid - Thread no. that called this routine.
33 INTEGER myThid
34
35 C LOCAL VARIABLES:
36 INTEGER i,j
37
38 DO j=1, sNy
39 DO i=1,sNx
40
41 IF (maskC(i,j,1,1,1) .EQ. 1.) THEN
42
43 C Ad hoc phase out wind stress if sufficient ice coverage
44 C similar idea as stressReduction used in thsice_main
45 IF (iceMask(i,j,1,1) .GT. windice_thres) THEN
46
47 fu_2D(i,j)= fu_2D(i,j) * (1. _d 0 -
48 & (iceMask(i,j,1,1) + iceMask(i-1,j,1,1)) * 0.5 _d 0 )
49 & / (1. _d 0 - windice_thres)
50 fv_2D(i,j)= fv_2D(i,j) * (1. _d 0 -
51 & (iceMask(i,j,1,1) + iceMask(i,j-1,1,1)) * 0.5 _d 0 )
52 & / (1. _d 0 - windice_thres)
53 wspeed_2D(i,j)= wspeed_2D(i,j)
54 & * (1. _d 0 - iceMask(i,j,1,1))
55 & / (1. _d 0 - windice_thres)
56 ENDIF
57
58 sum_runoff(i,j)= sum_runoff(i,j) + runoff_2D(i,j)
59 sum_precip(i,j)= sum_precip(i,j) +
60 & precipo_2D(i,j)*(1. _d 0-iceMask(i,j,1,1))
61 sum_evap(i,j)= sum_evap(i,j) +
62 & evapo_2D(i,j)*(1. _d 0-iceMask(i,j,1,1))
63 sum_qnet(i,j)= sum_qnet(i,j) +
64 & qneto_2D(i,j)*(1. _d 0-iceMask(i,j,1,1))
65 sum_fu(i,j)= sum_fu(i,j) + fu_2D(i,j)
66 sum_fv(i,j)= sum_fv(i,j) + fv_2D(i,j)
67 sum_wspeed(i,j)= sum_wspeed(i,j) + wspeed_2D(i,j)
68 sum_solarnet(i,j)= sum_solarnet(i,j) +
69 & solarnet_ocn_2D(i,j)*(1. _d 0-iceMask(i,j,1,1))
70 sum_slp(i,j)= sum_slp(i,j) + slp_2D(i,j)
71 sum_pCO2(i,j)= sum_pCO2(i,j) + pCO2_2D(i,j)
72 ENDIF
73 ENDDO
74 ENDDO
75
76 C PRINT *,'Sum_ocn_fluxes:',JBUGI,JBUGJ,fu_2D(JBUGI,JBUGJ),
77 C & fv_2D(JBUGI,JBUGJ), runoff_2D(JBUGI,JBUGJ),
78 C & precipo_2D(JBUGI,JBUGJ), evapo_2D(JBUGI,JBUGJ),
79 C & qneto_2D(JBUGI,JBUGJ)
80
81 RETURN
82 END
  ViewVC Help
Powered by ViewVC 1.1.22