Changeset 3821


Ignore:
Timestamp:
2013-02-25T15:13:39+01:00 (8 years ago)
Author:
acc
Message:

Branch dev_MERGE_2012. #1060. Minor alterations to icbutl.F90 and icbclv.F90 (ICeBerg trajectory component)

Location:
branches/2012/dev_MERGE_2012/NEMOGCM/NEMO/OPA_SRC/ICB
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/2012/dev_MERGE_2012/NEMOGCM/NEMO/OPA_SRC/ICB/icbclv.F90

    r3614 r3821  
    5757      ! this assumes that input is given as equivalent water flux so that pure water density is appropriate 
    5858 
    59       zfact = ( (1000._wp)**3 / ( NINT(rday) * nyear_len(1) ) ) * 1000._wp 
     59      zfact = ( (1000._wp)**3 / ( NINT(rday) * nyear_len(1) ) ) * 850._wp 
    6060      berg_grid%calving(:,:) = src_calving(:,:) * tmask_i(:,:) * zfact 
    6161 
  • branches/2012/dev_MERGE_2012/NEMOGCM/NEMO/OPA_SRC/ICB/icbutl.F90

    r3785 r3821  
    196196            ! since we're looking for four T points containing quadrant we're in of  
    197197            ! current T cell 
    198             ii = INT( pi     ) 
    199             ij = INT( pj      )    ! T-point 
     198            ii = MAX(1, INT( pi     )) 
     199            ij = MAX(1, INT( pj     ))    ! T-point 
    200200            zi = pi - REAL(ii,wp) 
    201201            zj = pj - REAL(ij,wp) 
    202202         CASE ( 'U' ) 
    203             ii = INT( pi-0.5 ) 
    204             ij = INT( pj      )    ! U-point 
     203            ii = MAX(1, INT( pi-0.5 )) 
     204            ij = MAX(1, INT( pj     ))    ! U-point 
    205205            zi = pi - 0.5 - REAL(ii,wp) 
    206206            zj = pj - REAL(ij,wp) 
    207207         CASE ( 'V' ) 
    208             ii = INT( pi     ) 
    209             ij = INT( pj -0.5 )    ! V-point 
     208            ii = MAX(1, INT( pi     )) 
     209            ij = MAX(1, INT( pj-0.5 ))    ! V-point 
    210210            zi = pi - REAL(ii,wp) 
    211211            zj = pj - 0.5 - REAL(ij,wp) 
    212212         CASE ( 'F' ) 
    213             ii = INT( pi-0.5 ) 
    214             ij = INT( pj -0.5 )    ! F-point 
     213            ii = MAX(1, INT( pi-0.5 )) 
     214            ij = MAX(1, INT( pj-0.5 ))    ! F-point 
    215215            zi = pi - 0.5 - REAL(ii,wp) 
    216216            zj = pj - 0.5 - REAL(ij,wp) 
     
    250250            ! since we're looking for four T points containing quadrant we're in of  
    251251            ! current T cell 
    252             ii = INT( pi     ) 
    253             ij = INT( pj      )    ! T-point 
     252            ii = MAX(1, INT( pi     )) 
     253            ij = MAX(1, INT( pj     ))    ! T-point 
    254254            zi = pi - REAL(ii,wp) 
    255255            zj = pj - REAL(ij,wp) 
    256256         CASE ( 'U' ) 
    257             ii = INT( pi-0.5 ) 
    258             ij = INT( pj      )    ! U-point 
     257            ii = MAX(1, INT( pi-0.5 )) 
     258            ij = MAX(1, INT( pj     ))    ! U-point 
    259259            zi = pi - 0.5 - REAL(ii,wp) 
    260260            zj = pj - REAL(ij,wp) 
    261261         CASE ( 'V' ) 
    262             ii = INT( pi     ) 
    263             ij = INT( pj -0.5 )    ! V-point 
     262            ii = MAX(1, INT( pi     )) 
     263            ij = MAX(1, INT( pj-0.5 ))    ! V-point 
    264264            zi = pi - REAL(ii,wp) 
    265265            zj = pj - 0.5 - REAL(ij,wp) 
    266266         CASE ( 'F' ) 
    267             ii = INT( pi-0.5 ) 
    268             ij = INT( pj -0.5 )    ! F-point 
     267            ii = MAX(1, INT( pi-0.5 )) 
     268            ij = MAX(1, INT( pj-0.5 ))    ! F-point 
    269269            zi = pi - 0.5 - REAL(ii,wp) 
    270270            zj = pj - 0.5 - REAL(ij,wp) 
     
    304304      ! since we're looking for four T points containing quadrant we're in of  
    305305      ! current T cell 
    306       ii = INT( pi     ) 
    307       ij = INT( pj      )    ! T-point 
     306      ii = MAX(1, INT( pi     )) 
     307      ij = MAX(1, INT( pj     ))    ! T-point 
    308308      zi = pi - REAL(ii,wp) 
    309309      zj = pj - REAL(ij,wp) 
     
    348348      !!---------------------------------------------------------------------- 
    349349      ! 
    350       ii = INT( pi )   ;   ij = INT( pj )    ! left bottom T-point (i,j) indices 
     350      ii = MAX(1, INT( pi ))   ;   ij = MAX(1, INT( pj ))            ! left bottom T-point (i,j) indices 
    351351 
    352352      ! fractional box spacing 
Note: See TracChangeset for help on using the changeset viewer.