1 |
C $Header: /u/gcmpack/MITgcm/pkg/seaice/seaice_cost_export.F,v 1.4 2010/03/16 00:23:18 jmc Exp $ |
2 |
C $Name: $ |
3 |
|
4 |
#include "SEAICE_OPTIONS.h" |
5 |
|
6 |
subroutine seaice_cost_export( myThid ) |
7 |
C *==========================================================* |
8 |
C | subroutine seaice_cost_export |
9 |
C | o Poor man ice export through Fram Strait |
10 |
C | Parameters currently hard-wired for arctic210x192x50 |
11 |
C *==========================================================* |
12 |
implicit none |
13 |
|
14 |
C == Global variables === |
15 |
#include "SIZE.h" |
16 |
#include "EEPARAMS.h" |
17 |
#include "PARAMS.h" |
18 |
#include "GRID.h" |
19 |
#include "SEAICE_SIZE.h" |
20 |
#include "SEAICE.h" |
21 |
#include "SEAICE_COST.h" |
22 |
|
23 |
C ======== Routine arguments ====================== |
24 |
C myThid - Thread number for this instance of the routine. |
25 |
integer myThid |
26 |
|
27 |
#ifdef ALLOW_SEAICE_COST_EXPORT |
28 |
C ========= Local variables ========================= |
29 |
integer isecbeg , isecend , jsec |
30 |
integer jsecbeg , jsecend , isec |
31 |
integer kmaxdepth |
32 |
integer i, j, k |
33 |
integer ig, jg |
34 |
integer bi, bj |
35 |
|
36 |
#if (defined (SEAICE_COST_EXPORT_FRAM)) |
37 |
parameter ( isec = 74, jsecbeg = 60, jsecend = 89 ) |
38 |
#elif (defined (SEAICE_COST_EXPORT_CANARCH)) |
39 |
parameter ( isec = 202, jsecbeg = 236, jsecend = 242 ) |
40 |
#else |
41 |
parameter ( isec = 1, jsecbeg = 1, jsecend = 1) |
42 |
#endif |
43 |
parameter ( kmaxdepth = 1 ) |
44 |
|
45 |
do bj=myByLo(myThid),myByHi(myThid) |
46 |
do bi=myBxLo(myThid),myBxHi(myThid) |
47 |
do i=1,sNx |
48 |
ig = myXGlobalLo-1+(bi-1)*sNx+i |
49 |
if (ig .eq. isec) then |
50 |
do k = 1, kmaxdepth |
51 |
do j=1,sNy |
52 |
jg = myYGlobalLo-1+(bj-1)*sNy+j |
53 |
if ((jg .ge. jsecbeg) .and. (jg .le. jsecend)) then |
54 |
objf_ice_export(bi,bj) = objf_ice_export(bi,bj) |
55 |
& - uHeffExportCell(i,j,bi,bj) |
56 |
cph( |
57 |
cph print *, 'ph-canarch pos i, ig, j, jg, X, Y ', |
58 |
cph & i, ig, j, jg, xC(i,j,bi,bj), yC(i,j,bi,bj) |
59 |
cph) |
60 |
end if |
61 |
enddo |
62 |
enddo |
63 |
end if |
64 |
end do |
65 |
c-- end of bi,bj loop |
66 |
end do |
67 |
end do |
68 |
|
69 |
#endif |
70 |
|
71 |
return |
72 |
end |