/[MITgcm]/MITgcm/verification/aim.5l_Equatorial_Channel/code/aim_surf_bc.F
ViewVC logotype

Diff of /MITgcm/verification/aim.5l_Equatorial_Channel/code/aim_surf_bc.F

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

revision 1.8 by jmc, Thu Jan 21 00:13:14 2010 UTC revision 1.9 by jmc, Mon Dec 12 19:05:41 2011 UTC
# Line 35  C-- size for MITgcm & Physics package : Line 35  C-- size for MITgcm & Physics package :
35  C-- MITgcm  C-- MITgcm
36  #include "EEPARAMS.h"  #include "EEPARAMS.h"
37  #include "PARAMS.h"  #include "PARAMS.h"
38  c #include "DYNVARS.h"  C_EqCh: start
39    #ifdef ALLOW_EXCH2
40    # include "W2_EXCH2_SIZE.h"
41    #endif /* ALLOW_EXCH2 */
42    #include "SET_GRID.h"
43    C_EqCh: end
44  #include "GRID.h"  #include "GRID.h"
45    c #include "DYNVARS.h"
46  c #include "SURFACE.h"  c #include "SURFACE.h"
47    
48  C-- Physics package  C-- Physics package
# Line 218  C-    Surface Albedo : (from F.M. FORDAT Line 224  C-    Surface Albedo : (from F.M. FORDAT
224  c_FM    DALB=ALBICE-ALBSEA  c_FM    DALB=ALBICE-ALBSEA
225          RSD=1. _d 0/SDALB          RSD=1. _d 0/SDALB
226          ALPHA= 2. _d 0*PI*(TYEAR+10. _d 0/365. _d 0)          ALPHA= 2. _d 0*PI*(TYEAR+10. _d 0/365. _d 0)
227    #ifdef ALLOW_INSOLATION
228            ZS = - SIN(OBLIQ * deg2rad) * COS(ALPHA)
229            ZC =   ASIN( ZS )
230            ZC =   COS(ZC)
231    #else /* ALLOW_INSOLATION */
232          RZEN = COS(ALPHA) * ( -23.45 _d 0 * deg2rad)          RZEN = COS(ALPHA) * ( -23.45 _d 0 * deg2rad)
233          ZC = COS(RZEN)          ZC = COS(RZEN)
234          ZS = SIN(RZEN)          ZS = SIN(RZEN)
235    #endif /* ALLOW_INSOLATION */
236          DO j=1,sNy          DO j=1,sNy
237           DO i=1,sNx           DO i=1,sNx
238  c_FM      SNOWC=MIN(1.,RSD*SNOW1(I,J))  c_FM      SNOWC=MIN(1.,RSD*SNOW1(I,J))
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22