Changeset 12377 for NEMO/trunk/src/OCE/DOM/domwri.F90
- Timestamp:
- 2020-02-12T15:39:06+01:00 (4 years ago)
- Location:
- NEMO/trunk
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk
- Property svn:externals
-
old new 3 3 ^/utils/build/mk@HEAD mk 4 4 ^/utils/tools@HEAD tools 5 ^/vendors/AGRIF/dev @HEAD ext/AGRIF5 ^/vendors/AGRIF/dev_r11615_ENHANCE-04_namelists_as_internalfiles_agrif@HEAD ext/AGRIF 6 6 ^/vendors/FCM@HEAD ext/FCM 7 7 ^/vendors/IOIPSL@HEAD ext/IOIPSL
-
- Property svn:externals
-
NEMO/trunk/src/OCE/DOM/domwri.F90
r11532 r12377 16 16 !! dom_stiff : diagnose maximum grid stiffness/hydrostatic consistency (s-coordinate) 17 17 !!---------------------------------------------------------------------- 18 ! 18 19 USE dom_oce ! ocean space and time domain 19 20 USE phycst , ONLY : rsmall … … 32 33 33 34 !! * Substitutions 34 # include " vectopt_loop_substitute.h90"35 # include "do_loop_substitute.h90" 35 36 !!---------------------------------------------------------------------- 36 37 !! NEMO/OCE 4.0 , NEMO Consortium (2018) … … 99 100 100 101 CALL dom_uniq( zprw, 'T' ) 101 DO jj = 1, jpj 102 DO ji = 1, jpi 103 zprt(ji,jj) = ssmask(ji,jj) * zprw(ji,jj) ! ! unique point mask 104 END DO 105 END DO ! ! unique point mask 102 DO_2D_11_11 103 zprt(ji,jj) = ssmask(ji,jj) * zprw(ji,jj) ! ! unique point mask 104 END_2D 106 105 CALL iom_rstput( 0, 0, inum, 'tmaskutil', zprt, ktype = jp_i1 ) 107 106 CALL dom_uniq( zprw, 'U' ) 108 DO jj = 1, jpj 109 DO ji = 1, jpi 110 zprt(ji,jj) = ssumask(ji,jj) * zprw(ji,jj) ! ! unique point mask 111 END DO 112 END DO 107 DO_2D_11_11 108 zprt(ji,jj) = ssumask(ji,jj) * zprw(ji,jj) ! ! unique point mask 109 END_2D 113 110 CALL iom_rstput( 0, 0, inum, 'umaskutil', zprt, ktype = jp_i1 ) 114 111 CALL dom_uniq( zprw, 'V' ) 115 DO jj = 1, jpj 116 DO ji = 1, jpi 117 zprt(ji,jj) = ssvmask(ji,jj) * zprw(ji,jj) ! ! unique point mask 118 END DO 119 END DO 112 DO_2D_11_11 113 zprt(ji,jj) = ssvmask(ji,jj) * zprw(ji,jj) ! ! unique point mask 114 END_2D 120 115 CALL iom_rstput( 0, 0, inum, 'vmaskutil', zprt, ktype = jp_i1 ) 121 116 !!gm ssfmask has been removed ==>> find another solution to defined fmaskutil … … 155 150 156 151 ! note that mbkt is set to 1 over land ==> use surface tmask 157 zprt(:,:) = ssmask(:,:) *REAL( mbkt(:,:) , wp )152 zprt(:,:) = REAL( mbkt(:,:) , wp ) 158 153 CALL iom_rstput( 0, 0, inum, 'mbathy', zprt, ktype = jp_i4 ) ! ! nb of ocean T-points 159 zprt(:,:) = ssmask(:,:) *REAL( mikt(:,:) , wp )154 zprt(:,:) = REAL( mikt(:,:) , wp ) 160 155 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-points163 156 ! ! vertical mesh 164 157 CALL iom_rstput( 0, 0, inum, 'e3t_1d', e3t_1d, ktype = jp_r8 ) ! ! scale factors
Note: See TracChangeset
for help on using the changeset viewer.