/[MITgcm]/MITgcm_contrib/ESMF/global_ocean.128x60x15/code/ocn_export_sst.F
ViewVC logotype

Annotation of /MITgcm_contrib/ESMF/global_ocean.128x60x15/code/ocn_export_sst.F

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


Revision 1.1 - (hide annotations) (download)
Tue Mar 30 03:58:58 2004 UTC (21 years, 4 months ago) by cnh
Branch: MAIN
CVS Tags: adoption_1_0_pre_A, HEAD
New test with different size

1 cnh 1.1 C $Header: /u/gcmpack/MITgcm_contrib/ESMF/global_ocean.128x64x15/code/ocn_export_sst.F,v 1.1.1.1 2004/02/15 22:28:19 cnh Exp $
2     C $Name: $
3    
4     #include "CPP_OPTIONS.h"
5    
6     CStartOfInterface
7     SUBROUTINE OCN_EXPORT_SST(myThid)
8     C /==========================================================\
9     C | SUBROUTINE OCN_EXPORT_SST |
10     C | o Routine for exporting ocean sea-surface temperature |
11     C | to the coupling layer. |
12     C |==========================================================|
13     C | This version talks to the MIT Coupler. It uses the MIT |
14     C | Coupler "checkpoint1" library calls. |
15     C \==========================================================/
16     IMPLICIT NONE
17    
18     C == Global variables ==
19     #include "SIZE.h"
20     #include "EEPARAMS.h"
21     #include "GRID.h"
22     #include "DYNVARS.h"
23     #include "OCNIDS.h"
24     #include "OCNCPL.h"
25     #ifdef ALLOW_SEAICE
26     #include "ICE.h"
27     #endif
28    
29     C == Routine arguments ==
30     C myThid - Thread number for this instance of the routine
31     INTEGER myThid
32     CEndOfInterface
33    
34     C == Local variables ==
35     C I,J :: Loop variables
36     INTEGER bi,bj
37     INTEGER I,J
38    
39     DO bj = myByLo(myThid), myByHi(myThid)
40     DO bi = myBxLo(myThid), myBxHi(myThid)
41    
42     C Send my sea-surface temperatur to coupling layer
43     C o Information is sent on the ocean model grid.
44     DO J=1,sNy
45     DO I=1,sNx
46     IF ( hFacC(I,J,1,bi,bj) .EQ. 0. ) THEN
47     SSTocn2cpl(I,J,bi,bj) = theta(I,J,1,bi,bj)
48     ELSE
49     SSTocn2cpl(I,J,bi,bj) = theta(I,J,1,bi,bj)
50     ENDIF
51     ENDDO
52     ENDDO
53     #ifdef ALLOW_SEAICE
54     DO J=1,sNy
55     DO I=1,sNx
56     IF ( iceMask(I,J,bi,bj) .NE. 0. ) THEN
57     SSTocn2cpl(I,J,bi,bj) = 0.
58     ENDIF
59     ENDDO
60     ENDDO
61     #endif
62     ENDDO
63     ENDDO
64    
65     C CALL COMPSEND_R8TILES( ocnSSTName,
66     C I sNx, OLx, sNy, OLy, 1, nSx, nSy, SSTocn2cpl )
67    
68     RETURN
69     END

  ViewVC Help
Powered by ViewVC 1.1.22