/[MITgcm]/MITgcm/pkg/dic/dic_ini_forcing.F
ViewVC logotype

Diff of /MITgcm/pkg/dic/dic_ini_forcing.F

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

revision 1.6 by stephd, Tue May 1 22:51:12 2007 UTC revision 1.9 by dfer, Fri Oct 26 21:08:13 2007 UTC
# Line 75  C First call requires that we initialize Line 75  C First call requires that we initialize
75  # endif  # endif
76    
77  #ifdef READ_PAR  #ifdef READ_PAR
78  #ifdef PAR_FROMSW  #ifdef USE_QSW
79          WRITE(msgBuf,'(2A)')          WRITE(msgBuf,'(2A)')
80       &       ' DIC_INI_FORCING: You can not use READ_PAR ',       &       ' DIC_INI_FORCING: You can not use READ_PAR ',
81       &       ' and PAR_FROMSW'       &       ' and USE_QSW'
82          CALL PRINT_ERROR( msgBuf, myThid )          CALL PRINT_ERROR( msgBuf, myThid )
83          STOP 'ABNORMAL END: S/R DIC_INI_FORCING'          STOP 'ABNORMAL END: S/R DIC_INI_FORCING'
84  #endif  #endif
# Line 89  c set reasonable values to those that ne Line 89  c set reasonable values to those that ne
89          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
90           DO j=1-Oly,sNy+Oly           DO j=1-Oly,sNy+Oly
91            DO i=1-Olx,sNx+Olx            DO i=1-Olx,sNx+Olx
92               WIND(i,j,bi,bj) = 5.d0*maskC(i,j,1,bi,bj)               WIND(i,j,bi,bj)   = 5. _d 0*maskC(i,j,1,bi,bj)
93               ATMOSP(i,j,bi,bj)   =1.d0*maskC(i,j,1,bi,bj)               AtmosP(i,j,bi,bj) = 1. _d 0*maskC(i,j,1,bi,bj)
94               SILICA(i,j,bi,bj)   =7.6838e-3*maskC(i,j,1,bi,bj)               SILICA(i,j,bi,bj) = 7.6838 _d -3*maskC(i,j,1,bi,bj)
95               FIce(i,j,bi,bj) =0.d0               FIce(i,j,bi,bj)   = 0. _d 0
96  #ifdef READ_PAR  #ifdef READ_PAR
97               PAR(i,j,bi,bj) =100.d0*maskC(i,j,1,bi,bj)               PAR(i,j,bi,bj)    = 100. _d 0*maskC(i,j,1,bi,bj)
98  #endif  #endif
99  #ifdef ALLOW_FE  #ifdef ALLOW_FE
100               InputFe(i,j,bi,bj)  = 1e-11*maskC(i,j,1,bi,bj)               InputFe(i,j,bi,bj)= 1. _d -11*maskC(i,j,1,bi,bj)
101  #endif  #endif
102            ENDDO            ENDDO
103           ENDDO           ENDDO
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22