/[MITgcm]/MITgcm_contrib/darwin2/pkg/darwin/darwin_fields_load.F
ViewVC logotype

Diff of /MITgcm_contrib/darwin2/pkg/darwin/darwin_fields_load.F

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

revision 1.7 by jahn, Wed Dec 4 21:21:49 2013 UTC revision 1.8 by jahn, Wed Dec 4 21:27:54 2013 UTC
# Line 4  C $Name$ Line 4  C $Name$
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
5  #include "PTRACERS_OPTIONS.h"  #include "PTRACERS_OPTIONS.h"
6  #include "DARWIN_OPTIONS.h"  #include "DARWIN_OPTIONS.h"
7    #ifdef ALLOW_EXF
8    #include "EXF_OPTIONS.h"
9    #endif
10    
11  #ifdef ALLOW_PTRACERS  #ifdef ALLOW_PTRACERS
12  #ifdef ALLOW_DARWIN  #ifdef ALLOW_DARWIN
# Line 37  c#include "GCHEM.h" Line 40  c#include "GCHEM.h"
40  #ifdef ALLOW_THSICE  #ifdef ALLOW_THSICE
41  #include "THSICE_VARS.h"  #include "THSICE_VARS.h"
42  #endif  #endif
43  #ifdef ALLOW_OFFLINE  #ifdef ALLOW_EXF
44  #include "OFFLINE.h"  #include "EXF_FIELDS.h"
45  #endif  #endif
46  #ifdef OASIM  #ifdef OASIM
47  #include "SPECTRAL.h"  #include "SPECTRAL.h"
# Line 458  C Line 461  C
461          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
462           DO j=1-Oly,sNy+Oly           DO j=1-Oly,sNy+Oly
463            DO i=1-Olx,sNx+Olx            DO i=1-Olx,sNx+Olx
464  cQQ need to include ice model here, if used             IF ( darwin_useiceFile ) THEN
465  #ifdef ALLOW_THSICE              FIce(i,j,bi,bj) = bWght*fice0(i,j,bi,bj)
466               FIce(i,j,bi,bj) = iceMask(i,j,bi,bj)       &                       +aWght*fice1(i,j,bi,bj)
467  #else             ENDIF
468  #ifdef ALLOW_SEAICE  #ifdef ALLOW_EXF
469               FIce(i,j,bi,bj) = AREA(i,j,bi,bj)  #ifdef EXF_SEAICE_FRACTION
470  #else             IF ( darwin_useareamask ) THEN
471             IF ( darwin_iceFile .NE. ' '  ) THEN              FIce(i,j,bi,bj) = areamask(i,j,bi,bj)
             fice(i,j,bi,bj)   = bWght*fice0(i,j,bi,bj)  
      &                         +aWght*fice1(i,j,bi,bj)  
            ELSE  
             fice(i,j,bi,bj)  = 0. _d 0  
472             ENDIF             ENDIF
473  #endif  #endif
474  #endif  #endif
475  c or use offline fields if provided  #ifdef ALLOW_THSICE
476  #ifdef ALLOW_OFFLINE             IF ( darwin_useiceMask ) THEN
477             IF (IceFile      .NE. ' ') THEN              FIce(i,j,bi,bj) = iceMask(i,j,bi,bj)
478               fice(i,j,bi,bj)  = ICEM(i,j,bi,bj)             ENDIF
479    #endif
480    #ifdef ALLOW_SEAICE
481               IF ( darwin_useAREA ) THEN
482                FIce(i,j,bi,bj) = AREA(i,j,bi,bj)
483             ENDIF             ENDIF
484  #endif  #endif
485    
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22