Changeset 12101 for utils/tools_UKMO_MERGE_2019/DOMAINcfg/src/domwri.F90
- Timestamp:
- 2019-12-06T18:45:39+01:00 (5 years ago)
- File:
-
- 1 copied
Legend:
- Unmodified
- Added
- Removed
-
utils/tools_UKMO_MERGE_2019/DOMAINcfg/src/domwri.F90
r12100 r12101 154 154 CALL iom_rstput( 0, 0, inum, 'ff_t', ff_t, ktype = jp_r8 ) 155 155 156 ! note that mbkt is set to 1 over land ==> use surface tmask156 ! note that mbkt and mikt is set to 1 over land ==> use surface tmask 157 157 zprt(:,:) = ssmask(:,:) * REAL( mbkt(:,:) , wp ) 158 CALL iom_rstput( 0, 0, inum, 'mbathy', zprt, ktype = jp_i4 ) ! ! nb of ocean T-points 158 CALL iom_rstput( 0, 0, inum, 'mbathy', zprt, ktype = jp_i4 ) ! ! nb of ocean T-points 159 CALL iom_rstput( 0, 0, inum, 'bathy_metry', bathy(:,:) * ssmask(:,:), ktype = jp_r8 ) ! ! bathymetry 159 160 zprt(:,:) = ssmask(:,:) * REAL( mikt(:,:) , wp ) 160 CALL iom_rstput( 0, 0, inum, 'misf', zprt, ktype = jp_i4 ) ! ! nb of ocean T-points 161 zprt(:,:) = ssmask(:,:) * REAL( risfdep(:,:) , wp ) 162 CALL iom_rstput( 0, 0, inum, 'isfdraft', zprt, ktype = jp_r8 ) ! ! nb of ocean T-points 163 ! ! vertical mesh 161 CALL iom_rstput( 0, 0, inum, 'misf', zprt, ktype = jp_i4 ) ! ! first wet level 162 CALL iom_rstput( 0, 0, inum, 'isfdraft' , risfdep(:,:) * ssmask(:,:), ktype = jp_r8 ) ! ! ice shelf draft 163 zprt(:,:) = ssmask(:,:) * REAL( mbkt(:,:) - mikt(:,:) + 1, wp ) 164 CALL iom_rstput( 0, 0, inum, 'mhw',zprt, ktype = jp_i4 ) 165 CALL iom_rstput( 0, 0, inum, 'hw' ,(bathy-risfdep)*ssmask, ktype = jp_r8 ) 166 167 ! ! vertical mesh 164 168 CALL iom_rstput( 0, 0, inum, 'e3t_0', e3t_0, ktype = jp_r8 ) ! ! scale factors 165 169 CALL iom_rstput( 0, 0, inum, 'e3u_0', e3u_0, ktype = jp_r8 ) … … 177 181 ENDIF 178 182 ! 179 ! IF( ll_wd ) CALL iom_rstput( 0, 0, inum, 'ht_0' , ht_0 , ktype = jp_r8 )180 181 183 ! ! ============================ 182 184 CALL iom_close( inum ) ! close the files
Note: See TracChangeset
for help on using the changeset viewer.