- Timestamp:
- 2010-11-25T20:03:49+01:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/DTA/dtasal.F90
r2392 r2436 21 21 USE in_out_manager ! I/O manager 22 22 USE phycst ! physical constants 23 #if defined key_orca_lev1024 USE lbclnk ! ocean lateral boundary conditions (or mpp link)25 #endif26 23 27 24 IMPLICIT NONE … … 63 60 #endif 64 61 REAL(wp):: zl 65 #if defined key_orca_lev1066 INTEGER :: ikr, ikw, ikt, jjk67 REAL(wp):: zfac68 #endif69 62 REAL(wp), DIMENSION(jpk) :: zsaldta ! auxiliary array for interpolation 70 63 CHARACTER(len=100) :: cn_dir ! Root directory for location of ssr files … … 99 92 CALL ctl_stop( 'dta_sal: unable to allocate sf_sal structure' ) ; RETURN 100 93 ENDIF 101 102 #if defined key_orca_lev10103 ALLOCATE( sf_sal(1)%fnow(jpi,jpj,jpkdta) )104 IF( sn_sal%ln_tint ) ALLOCATE( sf_sal(1)%fdta(jpi,jpj,jpkdta,2) )105 #else106 94 ALLOCATE( sf_sal(1)%fnow(jpi,jpj,jpk) ) 107 95 IF( sn_sal%ln_tint ) ALLOCATE( sf_sal(1)%fdta(jpi,jpj,jpk,2) ) 108 #endif109 96 ! ! fill sf_sal with sn_sal and control print 110 97 CALL fld_fill( sf_sal, (/ sn_sal /), cn_dir, 'dta_sal', 'Salinity data', 'namdta_sal' ) … … 166 153 #endif 167 154 168 #if defined key_orca_lev10169 DO jjk = 1, 5170 s_dta(:,:,jjk) = sf_sal(1)%fnow(:,:,1)171 ENDDO172 DO jk = 1, jpk-20,10173 ikr = INT(jk/10) + 1174 ikw = (ikr-1) *10 + 1175 ikt = ikw + 5176 DO jjk=ikt,ikt+9177 zfac = ( gdept_0(jjk ) - gdepw_0(ikt) ) / ( gdepw_0(ikt+10) - gdepw_0(ikt) )178 s_dta(:,:,jjk) = sf_sal(1)%fnow(:,:,ikr) + ( sf_sal(1)%fnow(:,:,ikr+1) - sf_sal(1)%fnow(:,:,ikr) ) * zfac179 END DO180 END DO181 DO jjk = jpk-5, jpk182 s_dta(:,:,jjk) = sf_sal(1)%fnow(:,:,jpkdta-1)183 END DO184 ! fill the overlap areas185 CALL lbc_lnk (s_dta(:,:,:),'Z',-999.,'no0')186 #else187 155 s_dta(:,:,:)=sf_sal(1)%fnow(:,:,:) 188 #endif189 156 190 157 IF( ln_sco ) THEN
Note: See TracChangeset
for help on using the changeset viewer.