- Timestamp:
- 2020-11-26T17:47:20+01:00 (4 years ago)
- Location:
- NEMO/branches/2020/dev_r12985_TOP-04_IMMERSE_BGC_interface
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12985_TOP-04_IMMERSE_BGC_interface
- Property svn:externals
-
old new 3 3 ^/utils/build/mk@HEAD mk 4 4 ^/utils/tools@HEAD tools 5 ^/vendors/AGRIF/dev @HEADext/AGRIF5 ^/vendors/AGRIF/dev_r12970_AGRIF_CMEMS ext/AGRIF 6 6 ^/vendors/FCM@HEAD ext/FCM 7 7 ^/vendors/IOIPSL@HEAD ext/IOIPSL 8 8 9 9 # SETTE 10 ^/utils/CI/sette@1 2931sette10 ^/utils/CI/sette@13559 sette
-
- Property svn:externals
-
NEMO/branches/2020/dev_r12985_TOP-04_IMMERSE_BGC_interface/src/SAO/sao_read.F90
r10069 r13892 10 10 USE netcdf 11 11 USE oce, ONLY: tsn, sshn 12 USE dom_oce, ONLY: n lci, nlcj, nimpp, njmpp, tmask12 USE dom_oce, ONLY: nimpp, njmpp, tmask 13 13 USE par_oce, ONLY: jpi, jpj, jpk 14 14 ! … … 94 94 IF (ifcst .LE. ntimes) THEN 95 95 ! Allocate temporary temperature array 96 ALLOCATE(temp_tn( nlci,nlcj,jpk))97 ALLOCATE(temp_sn( nlci,nlcj,jpk))98 ALLOCATE(temp_sshn( nlci,nlcj))96 ALLOCATE(temp_tn(jpi,jpj,jpk)) 97 ALLOCATE(temp_sn(jpi,jpj,jpk)) 98 ALLOCATE(temp_sshn(jpi,jpj)) 99 99 100 100 ! Set temp_tn, temp_sn to 0. … … 104 104 105 105 ! Create start and count arrays 106 start_n = (/ nimpp, njmpp, 1,ifcst /)107 count_n = (/ nlci, nlcj, jpk, 1/)108 start_s = (/ nimpp, njmpp , ifcst /)109 count_s = (/ nlci, nlcj, 1/)106 start_n = (/ nimpp, njmpp, 1, ifcst /) 107 count_n = (/ jpi, jpj, jpk, 1 /) 108 start_s = (/ nimpp, njmpp , ifcst /) 109 count_s = (/ jpi, jpj, 1 /) 110 110 111 111 ! Read information into temporary arrays … … 138 138 139 139 ! Mask out missing data index 140 tsn(1:nlci,1:nlcj,1:jpk,1) = temp_tn(:,:,:) * tmask(1:nlci,1:nlcj,1:jpk) 141 tsn(1:nlci,1:nlcj,1:jpk,2) = temp_sn(:,:,:) * tmask(1:nlci,1:nlcj,1:jpk) 142 sshn(1:nlci,1:nlcj) = temp_sshn(:,:) * tmask(1:nlci,1:nlcj,1) 143 144 ! Remove halo from tmask, tsn, sshn to prevent double obs counting 145 IF (jpi > nlci) THEN 146 tmask(nlci+1:,:,:) = 0 147 tsn(nlci+1:,:,:,1) = 0 148 tsn(nlci+1:,:,:,2) = 0 149 sshn(nlci+1:,:) = 0 150 END IF 151 IF (jpj > nlcj) THEN 152 tmask(:,nlcj+1:,:) = 0 153 tsn(:,nlcj+1:,:,1) = 0 154 tsn(:,nlcj+1:,:,2) = 0 155 sshn(:,nlcj+1:) = 0 156 END IF 157 140 tsn(1:jpi,1:jpj,1:jpk,1) = temp_tn(:,:,:) * tmask(1:jpi,1:jpj,1:jpk) 141 tsn(1:jpi,1:jpj,1:jpk,2) = temp_sn(:,:,:) * tmask(1:jpi,1:jpj,1:jpk) 142 sshn(1:jpi,1:jpj) = temp_sshn(:,:) * tmask(1:jpi,1:jpj,1) 143 158 144 ! Deallocate arrays 159 145 DEALLOCATE(temp_tn, temp_sn, temp_sshn)
Note: See TracChangeset
for help on using the changeset viewer.