- Timestamp:
- 2016-11-21T09:55:07+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_CNRS_2016/NEMOGCM/NEMO/LIM_SRC_2/limmsh_2.F90
r3625 r7277 70 70 ENDIF 71 71 72 IF( jphgr_msh == 2 .OR. jphgr_msh == 3 .OR. jphgr_msh == 5 ) &73 & CALL ctl_stop(' Coriolis parameter in LIM not set for f- or beta-plane' )74 75 72 !---------------------------------------------------------- 76 73 ! Initialization of local and some global (common) variables … … 79 76 njeq = INT( jpj / 2 ) !i bug mpp potentiel 80 77 njeqm1 = njeq - 1 81 82 fcor(:,:) = 2. * omega * SIN( gphit(:,:) * rad ) ! coriolis factor at T-point83 78 84 79 !i DO jj = 1, jpj … … 87 82 !i END DO 88 83 89 IF( f cor(1,1) * fcor(1,nlcj) < 0.e0 ) THEN ! local domain include both hemisphere84 IF( ff_t(1,1) * ff_t(1,nlcj) < 0.e0 ) THEN ! local domain include both hemisphere 90 85 l_jeq = .TRUE. 91 86 njeq = 1 92 DO WHILE ( njeq <= jpj .AND. f cor(1,njeq) < 0.e0 )87 DO WHILE ( njeq <= jpj .AND. ff_t(1,njeq) < 0.e0 ) 93 88 njeq = njeq + 1 94 89 END DO 95 90 IF(lwp ) WRITE(numout,*) ' the equator is inside the domain at about njeq = ', njeq 96 ELSEIF( f cor(1,1) < 0.e0 ) THEN91 ELSEIF( ff_t(1,1) < 0.e0 ) THEN 97 92 l_jeq = .FALSE. 98 93 njeq = jpj
Note: See TracChangeset
for help on using the changeset viewer.