New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 4599 for branches/2012 – NEMO

Changeset 4599 for branches/2012


Ignore:
Timestamp:
2014-03-27T10:00:19+01:00 (10 years ago)
Author:
pabouttier
Message:

Add W and F point treatments in grid_random.F90, see Ticket #1285

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/OPATAM_SRC/gridrandom.F90

    r3611 r4599  
    9797      !-------------------------------------------------------------------- 
    9898 
    99       IF ( cd_type == 'T' ) THEN 
     99      SELECT CASE ( cd_type ) 
     100      CASE( 'T' ) 
    100101  
    101102         pt2d(:,:) = pt2d(:,:) * tmask(:,:,1) 
    102103         CALL lbc_lnk( pt2d, 'T',  1.0 ) 
    103104 
    104       ELSEIF ( cd_type == 'U' ) THEN 
     105      CASE( 'U' )  
    105106 
    106107         pt2d(:,:) = pt2d(:,:) * umask(:,:,1) 
    107108         CALL lbc_lnk( pt2d, 'U', -1.0 ) 
    108109 
    109       ELSEIF ( cd_type == 'V' ) THEN 
     110      CASE ( 'V' ) 
    110111 
    111112         pt2d(:,:) = pt2d(:,:) * vmask(:,:,1) 
    112113         CALL lbc_lnk( pt2d, 'V', -1.0 ) 
    113114 
    114       ELSEIF ( cd_type == 'S' ) THEN 
     115      CASE ( 'S' ) !: AV: S ??? 
    115116         CALL lbc_lnk( pt2d, 'S', 1.0 ) 
    116                 
    117       ENDIF 
     117 
     118      CASE ( 'F' )          
     119       
     120         pt2d(:,:) = pt2d(:,:) * fmask(:,:,1) 
     121         CALL lbc_lnk( pt2d, 'F',  1.0 ) 
     122 
     123      CASE ( 'W' )          
     124       
     125         pt2d(:,:) = pt2d(:,:) * tmask(:,:,1) 
     126         CALL lbc_lnk( pt2d, 'W',  1.0 ) 
     127 
     128      END SELECT 
    118129             
    119130   END SUBROUTINE grid_2d_ran 
    … …  
    173184      !-------------------------------------------------------------------- 
    174185 
    175       IF ( cd_type == 'T' ) THEN 
     186      SELECT CASE ( cd_type ) 
     187      CASE( 'T' ) 
    176188 
    177189         pt3d(:,:,:) = pt3d(:,:,:) * tmask(:,:,:) 
    178190         CALL lbc_lnk( pt3d, 'T',  1.0 ) 
    179191          
    180       ELSEIF ( cd_type == 'U' ) THEN 
     192      CASE( 'U' ) 
    181193             
    182194         pt3d(:,:,:) = pt3d(:,:,:) * umask(:,:,:) 
    183195         CALL lbc_lnk( pt3d, 'U', -1.0 ) 
    184196             
    185       ELSEIF ( cd_type == 'V' ) THEN 
     197      CASE( 'V' ) 
    186198          
    187199         pt3d(:,:,:) = pt3d(:,:,:) * vmask(:,:,:) 
    188200         CALL lbc_lnk( pt3d, 'V', -1.0 ) 
    189201 
    190       ELSEIF ( cd_type == 'S' ) THEN 
     202      CASE( 'S' ) !: AV: S ??? 
     203 
    191204         CALL lbc_lnk( pt3d, 'S', 1.0 ) 
    192205 
    193       ENDIF 
     206      CASE( 'W' ) 
     207 
     208         pt3d(:,:,:) = pt3d(:,:,:) * tmask(:,:,:) 
     209         CALL lbc_lnk( pt3d, 'W',  1.0 ) 
     210 
     211      CASE( 'F' ) 
     212 
     213         pt3d(:,:,:) = pt3d(:,:,:) * fmask(:,:,:) 
     214         CALL lbc_lnk( pt3d, 'F',  1.0 ) 
     215 
     216      END SELECT 
    194217 
    195218   END SUBROUTINE grid_3d_ran 
Note: See TracChangeset for help on using the changeset viewer.