pkg/seaice/seaice_cost_export.F

1e22f5fc71 Patr*0001 #include "SEAICE_OPTIONS.h"
                0002
                0003       subroutine seaice_cost_export( myThid )
03c669d1ab Jean*0004 C     *==========================================================*
1e22f5fc71 Patr*0005 C     | subroutine seaice_cost_export
88f72205aa Jean*0006 C     | o Poor man ice export through Fram Strait
1e22f5fc71 Patr*0007 C     |   Parameters currently hard-wired for arctic210x192x50
03c669d1ab Jean*0008 C     *==========================================================*
1e22f5fc71 Patr*0009        implicit none
                0010
                0011 C     == Global variables ===
                0012 #include "SIZE.h"
                0013 #include "EEPARAMS.h"
                0014 #include "PARAMS.h"
                0015 #include "GRID.h"
03c669d1ab Jean*0016 #include "SEAICE_SIZE.h"
1e22f5fc71 Patr*0017 #include "SEAICE.h"
                0018 #include "SEAICE_COST.h"
                0019
                0020 C     ======== Routine arguments ======================
                0021 C     myThid - Thread number for this instance of the routine.
                0022       integer myThid
                0023
                0024 #ifdef ALLOW_SEAICE_COST_EXPORT
                0025 C     ========= Local variables =========================
                0026       integer    isecbeg      , isecend      , jsec
                0027       integer    jsecbeg      , jsecend      , isec
                0028       integer    kmaxdepth
                0029       integer i, j, k
                0030       integer ig, jg
                0031       integer bi, bj
                0032
57867cd15e Patr*0033 #if (defined (SEAICE_COST_EXPORT_FRAM))
                0034       parameter ( isec = 74, jsecbeg = 60, jsecend = 89 )
                0035 #elif (defined (SEAICE_COST_EXPORT_CANARCH))
                0036       parameter ( isec = 202, jsecbeg = 236, jsecend = 242 )
                0037 #else
                0038       parameter ( isec = 1, jsecbeg = 1, jsecend = 1)
                0039 #endif
1e22f5fc71 Patr*0040       parameter ( kmaxdepth = 1 )
                0041
                0042       do bj=myByLo(myThid),myByHi(myThid)
                0043        do bi=myBxLo(myThid),myBxHi(myThid)
                0044         do i=1,sNx
                0045          ig = myXGlobalLo-1+(bi-1)*sNx+i
                0046          if (ig .eq. isec) then
                0047           do k = 1, kmaxdepth
                0048            do j=1,sNy
                0049             jg = myYGlobalLo-1+(bj-1)*sNy+j
                0050             if ((jg .ge. jsecbeg) .and. (jg .le. jsecend)) then
                0051              objf_ice_export(bi,bj) = objf_ice_export(bi,bj)
57867cd15e Patr*0052      &              - uHeffExportCell(i,j,bi,bj)
                0053 cph(
                0054 cph             print *, 'ph-canarch pos i, ig, j, jg, X, Y ',
                0055 cph     &        i, ig, j, jg, xC(i,j,bi,bj), yC(i,j,bi,bj)
                0056 cph)
1e22f5fc71 Patr*0057             end if
                0058            enddo
                0059           enddo
                0060          end if
                0061         end do
                0062 c--   end of bi,bj loop
                0063        end do
                0064       end do
                0065
                0066 #endif
7c7521a1da Jean*0067
03c669d1ab Jean*0068       return
1e22f5fc71 Patr*0069       end