Changeset 6649 for branches/UKMO/dev_r5518_RH_MEDUSA_Stable
- Timestamp:
- 2016-06-01T10:31:19+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_RH_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/trcini.F90
r6648 r6649 155 155 ELSE 156 156 ! 157 #if !defined key_medusa158 157 IF( ln_trcdta .AND. nb_trcdta > 0 ) THEN ! Initialisation of tracer from a file that may also be used for damping 159 158 ! … … 165 164 CALL trc_dta( nit000, sf_trcdta(jl),rf_trfac(jl) ) ! read tracer data at nit000 166 165 ztrcdta(:,:,:) = sf_trcdta(jl)%fnow(:,:,:) 167 trn(:,:,:,jn) = ztrcdta(:,:,:) * tmask(:,:,:) 166 trn(:,:,:,jn) = ztrcdta(:,:,:) * tmask(:,:,:) 168 167 IF( .NOT.ln_trcdmp .AND. .NOT.ln_trcdmp_clo ) THEN !== deallocate data structure ==! 169 168 ! (data used only for initialisation) … … 177 176 CALL wrk_dealloc( jpi, jpj, jpk, ztrcdta ) 178 177 ENDIF 179 #else180 IF( ln_trcdta .AND. nb_trcdta > 0 ) THEN ! Initialisation of tracer from a file that may also be used for damping181 182 IF(lwp) WRITE(numout,*) 'AXY: calling wrk_alloc'183 IF(lwp) CALL flush(numout)184 185 CALL wrk_alloc( jpi, jpj, jpk, ztrcdta ) ! Memory allocation186 !187 IF(lwp) WRITE(numout,*) 'AXY: calling trc_dta'188 IF(lwp) CALL flush(numout)189 190 DO jn = 1, jptra191 IF( ln_trc_ini(jn) ) THEN ! update passive tracers arrays with input data read from file192 jl = n_trc_index(jn)193 CALL trc_dta( nit000, sf_trcdta(jl),rf_trfac(jl) ) ! read tracer data at nit000194 ztrcdta(:,:,:) = sf_trcdta(jl)%fnow(:,:,:)195 trn(:,:,:,jn) = ztrcdta(:,:,:) * tmask(:,:,:)196 IF( .NOT.ln_trcdmp .AND. .NOT. ln_trcdmp_clo ) THEN !== deallocate data structure ==!197 ! (data used only for initialisation)198 IF(lwp) WRITE(numout,*) 'trc_dta: deallocate data arrays as they are only used to initialize the run'199 DEALLOCATE( sf_trcdta(jl)%fnow ) ! arrays in the structure200 IF( sf_trcdta(jl)%ln_tint ) DEALLOCATE( sf_trcdta(jl)%fdta )201 !202 ENDIF203 ENDIF204 ENDDO205 206 IF(lwp) WRITE(numout,*) 'AXY: calling wrk_dealloc'207 IF(lwp) CALL flush(numout)208 CALL wrk_dealloc( jpi, jpj, jpk, ztrcdta )209 ENDIF210 #endif211 178 ! 212 179 # if defined key_debug_medusa
Note: See TracChangeset
for help on using the changeset viewer.