- Timestamp:
- 2011-06-27T13:18:25+02:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/dev_r2787_LOCEAN3_TRA_TRP/NEMOGCM/NEMO/OPA_SRC/IOM/restart.F90
r2528 r2789 24 24 USE trdmld_oce ! ocean active mixed layer tracers trends variables 25 25 USE domvvl ! variable volume 26 USE traswp ! swap from 4D T-S to 3D T & S and vice versa27 26 28 27 IMPLICIT NONE … … 117 116 CALL iom_rstput( kt, nitrst, numrow, 'ub' , ub ) ! before fields 118 117 CALL iom_rstput( kt, nitrst, numrow, 'vb' , vb ) 119 CALL iom_rstput( kt, nitrst, numrow, 'tb' , t b)120 CALL iom_rstput( kt, nitrst, numrow, 'sb' , sb)118 CALL iom_rstput( kt, nitrst, numrow, 'tb' , tsb(:,:,:,jp_tem) ) 119 CALL iom_rstput( kt, nitrst, numrow, 'sb' , tsb(:,:,:,jp_sal) ) 121 120 CALL iom_rstput( kt, nitrst, numrow, 'rotb' , rotb ) 122 121 CALL iom_rstput( kt, nitrst, numrow, 'hdivb' , hdivb ) … … 126 125 CALL iom_rstput( kt, nitrst, numrow, 'un' , un ) ! now fields 127 126 CALL iom_rstput( kt, nitrst, numrow, 'vn' , vn ) 128 CALL iom_rstput( kt, nitrst, numrow, 'tn' , t n)129 CALL iom_rstput( kt, nitrst, numrow, 'sn' , sn)127 CALL iom_rstput( kt, nitrst, numrow, 'tn' , tsn(:,:,:,jp_tem) ) 128 CALL iom_rstput( kt, nitrst, numrow, 'sn' , tsn(:,:,:,jp_sal) ) 130 129 CALL iom_rstput( kt, nitrst, numrow, 'rotn' , rotn ) 131 130 CALL iom_rstput( kt, nitrst, numrow, 'hdivn' , hdivn ) … … 186 185 CALL iom_get( numror, jpdom_autoglo, 'ub' , ub ) ! before fields 187 186 CALL iom_get( numror, jpdom_autoglo, 'vb' , vb ) 188 CALL iom_get( numror, jpdom_autoglo, 'tb' , t b)189 CALL iom_get( numror, jpdom_autoglo, 'sb' , sb)187 CALL iom_get( numror, jpdom_autoglo, 'tb' , tsb(:,:,:,jp_tem) ) 188 CALL iom_get( numror, jpdom_autoglo, 'sb' , tsb(:,:,:,jp_sal) ) 190 189 CALL iom_get( numror, jpdom_autoglo, 'rotb' , rotb ) 191 190 CALL iom_get( numror, jpdom_autoglo, 'hdivb' , hdivb ) … … 195 194 CALL iom_get( numror, jpdom_autoglo, 'un' , un ) ! now fields 196 195 CALL iom_get( numror, jpdom_autoglo, 'vn' , vn ) 197 CALL iom_get( numror, jpdom_autoglo, 'tn' , t n)198 CALL iom_get( numror, jpdom_autoglo, 'sn' , sn)196 CALL iom_get( numror, jpdom_autoglo, 'tn' , tsn(:,:,:,jp_tem) ) 197 CALL iom_get( numror, jpdom_autoglo, 'sn' , tsn(:,:,:,jp_sal) ) 199 198 CALL iom_get( numror, jpdom_autoglo, 'rotn' , rotn ) 200 199 CALL iom_get( numror, jpdom_autoglo, 'hdivn' , hdivn ) … … 205 204 CALL iom_get( numror, jpdom_autoglo, 'rhd' , rhd ) ! now in situ density anomaly 206 205 ELSE 207 CALL tra_swap208 206 CALL eos( tsn, rhd ) ! compute rhd 209 207 ENDIF … … 211 209 ! 212 210 IF( neuler == 0 ) THEN ! Euler restart (neuler=0) 213 tb (:,:,:) = tn (:,:,:) ! all before fields set to now values 214 sb (:,:,:) = sn (:,:,:) 215 ub (:,:,:) = un (:,:,:) 216 vb (:,:,:) = vn (:,:,:) 217 rotb (:,:,:) = rotn (:,:,:) 218 hdivb(:,:,:) = hdivn(:,:,:) 219 sshb (:,:) = sshn (:,:) 211 tsb (:,:,:,:) = tsn (:,:,:,:) ! all before fields set to now values 212 ub (:,:,:) = un (:,:,:) 213 vb (:,:,:) = vn (:,:,:) 214 rotb (:,:,:) = rotn (:,:,:) 215 hdivb(:,:,:) = hdivn(:,:,:) 216 sshb (:,:) = sshn (:,:) 220 217 IF( lk_vvl ) THEN 221 218 DO jk = 1, jpk
Note: See TracChangeset
for help on using the changeset viewer.