- Timestamp:
- 2017-12-01T14:51:50+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r8126_LIM3_couple/NEMOGCM/NEMO/OPA_SRC/SBC/geo2ocean.F90
r8877 r8878 27 27 PRIVATE 28 28 29 PUBLIC repcmo ! called in sbccpl 29 30 PUBLIC rot_rep ! called in sbccpl, fldread, and cyclone 30 31 PUBLIC geo2oce ! called in sbccpl … … 51 52 !!---------------------------------------------------------------------- 52 53 CONTAINS 54 55 SUBROUTINE repcmo ( pxu1, pyu1, pxv1, pyv1, & 56 px2 , py2 , kchoix ) 57 !!---------------------------------------------------------------------- 58 !! *** ROUTINE repcmo *** 59 !! 60 !! ** Purpose : Change vector componantes from a geographic grid to a 61 !! stretched coordinates grid. 62 !! 63 !! ** Method : Initialization of arrays at the first call. 64 !! 65 !! ** Action : - px2 : first componante (defined at u point) 66 !! - py2 : second componante (defined at v point) 67 !!---------------------------------------------------------------------- 68 REAL(wp), INTENT(in ), DIMENSION(jpi,jpj) :: pxu1, pyu1 ! geographic vector componantes at u-point 69 REAL(wp), INTENT(in ), DIMENSION(jpi,jpj) :: pxv1, pyv1 ! geographic vector componantes at v-point 70 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: px2 ! i-componante (defined at u-point) 71 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: py2 ! j-componante (defined at v-point) 72 !!---------------------------------------------------------------------- 73 INTEGER, INTENT( IN ) :: & 74 kchoix ! type of transformation 75 ! = 1 change from geographic to model grid. 76 ! =-1 change from model to geographic grid 77 !!---------------------------------------------------------------------- 78 79 SELECT CASE (kchoix) 80 CASE ( 1) 81 ! Change from geographic to stretched coordinate 82 ! ---------------------------------------------- 83 84 CALL rot_rep( pxu1, pyu1, 'U', 'en->i',px2 ) 85 CALL rot_rep( pxv1, pyv1, 'V', 'en->j',py2 ) 86 CASE (-1) 87 ! Change from stretched to geographic coordinate 88 ! ---------------------------------------------- 89 90 CALL rot_rep( pxu1, pyu1, 'U', 'ij->e',px2 ) 91 CALL rot_rep( pxv1, pyv1, 'V', 'ij->n',py2 ) 92 END SELECT 93 94 END SUBROUTINE repcmo 53 95 54 96 SUBROUTINE rot_rep ( pxin, pyin, cd_type, cdtodo, prot )
Note: See TracChangeset
for help on using the changeset viewer.