Changeset 11371 for NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/ICE/icedyn.F90
- Timestamp:
- 2019-07-30T16:16:42+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/ICE/icedyn.F90
r11362 r11371 75 75 !! 76 76 INTEGER :: ji, jj ! dummy loop indices 77 REAL(wp) :: zmiss_val ! missing value retrieved from xios78 77 REAL(wp) :: zcoefu, zcoefv 79 REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: zdivu_i , zmsk0078 REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: zdivu_i 80 79 !!-------------------------------------------------------------------- 81 80 ! … … 89 88 ENDIF 90 89 ! 91 ! get missing value from xml92 CALL iom_miss_val( "icethic", zmiss_val )93 !94 90 ! retrieve thickness from volume for landfast param. and UMx advection scheme 95 91 WHERE( a_i(:,:,:) >= epsi20 ) … … 159 155 CASE ( np_dynADV1D , np_dynADV2D ) 160 156 161 ALLOCATE( zdivu_i(jpi,jpj), zmsk00(jpi,jpj) ) 162 163 DO jj = 1, jpj 164 DO ji = 1, jpi 165 zmsk00(ji,jj) = MAX( 0._wp , SIGN( 1._wp , at_i(ji,jj) - epsi06 ) ) ! 1 if ice , 0 if no ice 166 END DO 167 END DO 157 ALLOCATE( zdivu_i(jpi,jpj) ) 168 158 DO jj = 2, jpjm1 169 159 DO ji = 2, jpim1 … … 174 164 CALL lbc_lnk( 'icedyn', zdivu_i, 'T', 1. ) 175 165 ! output 176 CALL iom_put( "icediv" , zdivu_i * zmsk00 + zmiss_val * ( 1._wp - zmsk00 ))177 178 DEALLOCATE( zdivu_i , zmsk00)166 CALL iom_put( 'icediv' , zdivu_i ) 167 168 DEALLOCATE( zdivu_i ) 179 169 180 170 END SELECT
Note: See TracChangeset
for help on using the changeset viewer.