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 11371 for NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/ICE/icedyn.F90 – NEMO

Ignore:
Timestamp:
2019-07-30T16:16:42+02:00 (5 years ago)
Author:
clem
Message:

continue cleaning the ice outputs

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/ICE/icedyn.F90

    r11362 r11371  
    7575      !! 
    7676      INTEGER  ::   ji, jj        ! dummy loop indices 
    77       REAL(wp) ::   zmiss_val     ! missing value retrieved from xios  
    7877      REAL(wp) ::   zcoefu, zcoefv 
    79       REAL(wp), ALLOCATABLE, DIMENSION(:,:) ::   zdivu_i, zmsk00 
     78      REAL(wp), ALLOCATABLE, DIMENSION(:,:) ::   zdivu_i 
    8079      !!-------------------------------------------------------------------- 
    8180      ! 
     
    8988      ENDIF 
    9089      !                       
    91       ! get missing value from xml 
    92       CALL iom_miss_val( "icethic", zmiss_val ) 
    93       ! 
    9490      ! retrieve thickness from volume for landfast param. and UMx advection scheme 
    9591      WHERE( a_i(:,:,:) >= epsi20 ) 
     
    159155         CASE ( np_dynADV1D , np_dynADV2D ) 
    160156 
    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) ) 
    168158            DO jj = 2, jpjm1 
    169159               DO ji = 2, jpim1 
     
    174164            CALL lbc_lnk( 'icedyn', zdivu_i, 'T', 1. ) 
    175165            ! 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 ) 
    179169 
    180170         END SELECT 
Note: See TracChangeset for help on using the changeset viewer.