/[MITgcm]/MITgcm/pkg/icefront/icefront_init_fixed.F
ViewVC logotype

Diff of /MITgcm/pkg/icefront/icefront_init_fixed.F

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

revision 1.2 by dimitri, Fri Jan 22 00:51:54 2010 UTC revision 1.3 by dimitri, Fri Jan 29 01:05:27 2010 UTC
# Line 16  C     === Global variables === Line 16  C     === Global variables ===
16  #include "SIZE.h"  #include "SIZE.h"
17  #include "EEPARAMS.h"  #include "EEPARAMS.h"
18  #include "PARAMS.h"  #include "PARAMS.h"
19    #include "GRID.h"
20    #include "ICEFRONT.h"
21    
22  C     === Routine arguments ===  C     === Routine arguments ===
23  C     myThid -  Number of this instance of ICEFRONT_INIT_FIXED  C     myThid -  Number of this instance of ICEFRONT_INIT_FIXED
# Line 23  C     myThid -  Number of this instance Line 25  C     myThid -  Number of this instance
25    
26  #ifdef ALLOW_ICEFRONT  #ifdef ALLOW_ICEFRONT
27  C     === Local variables ===  C     === Local variables ===
28    C     I,J,K,bi,bj - Loop counters
29          INTEGER I, J, K, bi, bj
30          INTEGER ISinterface
31    
32          IF ( ICEFRONTlengthFile .NE. ' ' ) THEN
33           CALL READ_FLD_XY_RS( ICEFRONTlengthFile, ' ',
34         &      icefrontlength, 0, myThid )
35           _EXCH_XY_RS( icefrontlength, myThid )
36          ENDIF
37    
38          IF ( ICEFRONTdepthFile .NE. ' ' ) THEN
39           _BARRIER
40           CALL READ_FLD_XY_RS( ICEFRONTdepthFile, ' ',
41         &      R_icefront, 0, myThid )
42           _EXCH_XY_RS( R_icefront, myThid )
43          ENDIF
44    
45          DO bj = myByLo(myThid), myByHi(myThid)
46           DO bi = myBxLo(myThid), myBxHi(myThid)
47            DO J = 1-OLy, sNy+OLy
48             DO I = 1-OLx, sNx+OLx
49              K_icefront(i,j,bi,bj) = 0
50              DO K = 1 , Nr
51               IF ( ABS(R_icefront(I,J,bi,bj)) .GT. ABS(rF(K)))
52         &          K_icefront(I,J,bi,bj) = K
53              ENDDO
54             ENDDO
55            ENDDO
56           ENDDO
57          ENDDO
58    
59    #undef ALLOW_ICEFRONT_DEBUG
60    #ifdef ALLOW_ICEFRONT_DEBUG
61          DO bj = myByLo(myThid), myByHi(myThid)
62           DO bi = myBxLo(myThid), myBxHi(myThid)
63             DO J = 1, sNy
64              DO I = 1, sNx
65    C           IsInterface=0
66               IF (ICEFRONTlength(I,J,bi,bj) .GT. 0. _d 0) THEN
67    C              print*, 'IsInterface=', '2' , ',xuyun'
68                  IsInterface=Isinterface + K_icefront(I,J,bi,bj)
69               ENDIF
70              ENDDO
71             ENDDO
72            ENDDO
73           ENDDO
74           print*, 'Interface # =', IsInterface
75    #endif /* ALLOW_ICEFRONT_DEBUG */
76    
77  #ifdef ALLOW_DIAGNOSTICS  #ifdef ALLOW_DIAGNOSTICS
78        IF ( useDiagnostics ) THEN        IF ( useDiagnostics ) THEN
79          CALL ICEFRONT_DIAGNOSTICS_INIT( myThid )         CALL ICEFRONT_DIAGNOSTICS_INIT(myThid)
80        ENDIF        ENDIF
81  #endif /* ALLOW_DIAGNOSTICS */  #endif /* ALLOW_DIAGNOSTICS */
82  #endif /* ALLOW_ICEFRONT */  #endif /* ALLOW_ICEFRONT */
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22