- Timestamp:
- 02/28/20 10:07:19 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/trunk/src/sphere/geometry.f90
r963 r1014 227 227 INTEGER :: ind,i,j,n 228 228 REAL(rstd) :: schmidt_factor, schmidt_lon, schmidt_lat 229 LOGICAL :: read_metric=.FALSE. 230 231 CALL getin('read_metric', read_metric) 229 232 230 233 ! Schmidt transform parameters … … 241 244 schmidt_lat = schmidt_lat * pi/180. 242 245 243 DO ind=1,ndomain 244 IF (.NOT. assigned_domain(ind) .OR. .NOT. is_omp_level_master) CYCLE 245 CALL swap_dimensions(ind) 246 CALL swap_geometry(ind) 247 DO j=jj_begin,jj_end 248 DO i=ii_begin,ii_end 249 n=(j-1)*iim+i 250 CALL schmidt_transform(xyz_i(n,:), schmidt_factor, schmidt_lon, schmidt_lat) 246 IF (.NOT. read_metric) THEN 247 DO ind=1,ndomain 248 IF (.NOT. assigned_domain(ind) .OR. .NOT. is_omp_level_master) CYCLE 249 CALL swap_dimensions(ind) 250 CALL swap_geometry(ind) 251 DO j=jj_begin,jj_end 252 DO i=ii_begin,ii_end 253 n=(j-1)*iim+i 254 CALL schmidt_transform(xyz_i(n,:), schmidt_factor, schmidt_lon, schmidt_lat) 255 ENDDO 251 256 ENDDO 252 257 ENDDO 253 ENDDO 258 ENDIF 259 254 260 END SUBROUTINE remap_schmidt_loc 255 261 … … 263 269 USE getin_mod 264 270 USE omp_para 271 USE checksum_mod 265 272 IMPLICIT NONE 266 273 INTEGER :: nb_it=0 … … 271 278 REAL(rstd) :: sum 272 279 LOGICAL :: check 273 280 LOGICAL :: read_metric=.FALSE. 274 281 275 282 CALL getin('optim_it',nb_it) 283 CALL getin('read_metric', read_metric) 284 IF (read_metric) nb_it=0 276 285 277 286 DO ind=1,ndomain … … 425 434 USE getin_mod 426 435 USE omp_para 436 USE checksum_mod 427 437 IMPLICIT NONE 428 438 … … 446 456 ! so that XIOS gets the right values 447 457 CALL update_domain 448 458 449 459 DO ind=1,ndomain 450 460 IF (.NOT. assigned_domain(ind) .OR. .NOT. is_omp_level_master ) CYCLE
Note: See TracChangeset
for help on using the changeset viewer.