Changeset 1466 for trunk/NEMO
- Timestamp:
- 2009-06-11T09:48:51+02:00 (15 years ago)
- Location:
- trunk/NEMO/OPA_SRC/SBC
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/SBC/sbc_ice.F90
r1465 r1466 24 24 LOGICAL , PUBLIC, PARAMETER :: lk_lim2 = .TRUE. !: LIM-2 ice model 25 25 LOGICAL , PUBLIC, PARAMETER :: lk_lim3 = .FALSE. !: no LIM-3 26 CHARACTER(len=1), PUBLIC :: ci ce_grid= 'B' !: 'B'-grid ice-velocity26 CHARACTER(len=1), PUBLIC :: cigr_type = 'B' !: 'B'-grid ice-velocity 27 27 # endif 28 28 # if defined key_lim3 29 29 LOGICAL , PUBLIC, PARAMETER :: lk_lim2 = .FALSE. !: no LIM-2 30 30 LOGICAL , PUBLIC, PARAMETER :: lk_lim3 = .TRUE. !: LIM-3 ice model 31 CHARACTER(len=1), PUBLIC :: ci ce_grid= 'C' !: 'C'-grid ice-velocity31 CHARACTER(len=1), PUBLIC :: cigr_type = 'C' !: 'C'-grid ice-velocity 32 32 # endif 33 33 … … 57 57 LOGICAL , PUBLIC, PARAMETER :: lk_lim2 = .FALSE. !: no LIM-2 ice model 58 58 LOGICAL , PUBLIC, PARAMETER :: lk_lim3 = .FALSE. !: no LIM-3 ice model 59 CHARACTER(len=1), PUBLIC :: ci ce_grid= '-' !: no grid ice-velocity59 CHARACTER(len=1), PUBLIC :: cigr_type = '-' !: no grid ice-velocity 60 60 61 61 #endif -
trunk/NEMO/OPA_SRC/SBC/sbccpl.F90
r1465 r1466 645 645 !! ** Method : transform the received stress from the atmosphere into 646 646 !! an atmosphere-ice stress in the (i,j) ocean referencial 647 !! and at the velocity point of the sea-ice model (ci ce_grid):647 !! and at the velocity point of the sea-ice model (cigr_type): 648 648 !! 'C'-grid : i- (j-) components given at U- (V-) point 649 649 !! 'B'-grid : both components given at I-point … … 660 660 !! first as 2 components on the sphere 661 661 !! second as 2 components oriented along the local grid 662 !! third as 2 components on the ci ce_gridpoint662 !! third as 2 components on the cigr_type point 663 663 !! 664 664 !! In 'oce and ice' case, only one vector stress field … … 668 668 !! transformation is done and only if the ice-grid is a 'B'-grid. 669 669 !! 670 !! ** Action : return ptau_i, ptau_j, the stress over the ice at ci ce_gridpoint670 !! ** Action : return ptau_i, ptau_j, the stress over the ice at cigr_type point 671 671 !!---------------------------------------------------------------------- 672 672 REAL(wp), INTENT(out), DIMENSION(jpi,jpj) :: p_taui ! i- & j-components of atmos-ice stress [N/m2] … … 730 730 ! 731 731 ! j+1 j -----V---F 732 ! ice stress on ice velocity point (ci ce_grid) ! |732 ! ice stress on ice velocity point (cigr_type) ! | 733 733 ! (C-grid ==>(U,V) or B-grid ==> I) j | T U 734 734 ! | | … … 737 737 ! i-1 i i 738 738 ! i i+1 (for I) 739 SELECT CASE ( ci ce_grid)739 SELECT CASE ( cigr_type ) 740 740 ! 741 741 CASE( 'B' ) ! B-grid ==> I … … 1046 1046 END DO 1047 1047 CASE( 'weighted oce and ice' ) 1048 IF( ci ce_grid== 'C' ) THEN ! 'C'-grid ice velocity1048 IF( cigr_type == 'C' ) THEN ! 'C'-grid ice velocity 1049 1049 DO jj = 2, jpjm1 1050 1050 DO ji = fs_2, fs_jpim1 ! vector opt. … … 1069 1069 CALL lbc_lnk( zitx1, 'T', -1. ) ; CALL lbc_lnk( zity1, 'T', -1. ) 1070 1070 CASE( 'mixed oce-ice' ) 1071 IF( ci ce_grid== 'C' ) THEN ! 'C'-grid ice velocity1071 IF( cigr_type == 'C' ) THEN ! 'C'-grid ice velocity 1072 1072 DO jj = 2, jpjm1 1073 1073 DO ji = fs_2, fs_jpim1 ! vector opt.
Note: See TracChangeset
for help on using the changeset viewer.