- Timestamp:
- 2021-06-04T18:00:42+02:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2021/ticket2680_C1D_PAPA/src/OCE/SBC/fldread.F90
r14437 r14954 38 38 USE ioipsl , ONLY : ymds2ju, ju2ymds ! for calendar 39 39 USE lib_mpp ! MPP library 40 USE lbclnk ! ocean lateral boundary conditions ( C1Dcase)40 USE lbclnk ! ocean lateral boundary conditions (online interpolation case) 41 41 42 42 IMPLICIT NONE … … 356 356 INTEGER :: idvar ! variable ID 357 357 INTEGER :: idmspc ! number of spatial dimensions 358 LOGICAL :: lmoor ! C1D case: point data359 358 REAL(wp), DIMENSION(:,:,:), POINTER :: dta_alias ! short cut 360 359 !!--------------------------------------------------------------------- … … 374 373 CALL lbc_lnk( 'fldread', dta_alias(:,:,:), sdjf%cltype, sdjf%zsgn, kfillmode = jpfillcopy ) 375 374 ELSE ! default case 376 ! C1D case: If product of spatial dimensions == ipk, then x,y are of377 ! size 1 (point/mooring data): this must be read onto the central grid point378 375 idvar = iom_varid( sdjf%num, sdjf%clvar ) 379 376 idmspc = iom_file ( sdjf%num )%ndims( idvar ) 380 377 IF( iom_file( sdjf%num )%luld( idvar ) ) idmspc = idmspc - 1 ! id of the last spatial dimension 381 lmoor = ( idmspc == 0 .OR. PRODUCT( iom_file( sdjf%num )%dimsz( 1:MAX(idmspc,1) ,idvar ) ) == ipk ) 382 ! 383 IF( lk_c1d .AND. lmoor ) THEN 384 CALL iom_get( sdjf%num, jpdom_unknown, sdjf%clvar, dta_alias(2,2,:), sdjf%nrec(1,iaa) ) ! jpdom_unknown -> no lbc_lnk 385 CALL lbc_lnk( 'fldread', dta_alias(:,:,:), 'T', 1., kfillmode = jpfillcopy ) 386 ELSE 387 CALL iom_get( sdjf%num, jpdom_global, sdjf%clvar, dta_alias(:,:,:), sdjf%nrec(1,iaa), & 388 & sdjf%cltype, sdjf%zsgn, kfill = jpfillcopy ) 389 ENDIF 378 CALL iom_get( sdjf%num, jpdom_global, sdjf%clvar, dta_alias(:,:,:), sdjf%nrec(1,iaa), & 379 & sdjf%cltype, sdjf%zsgn, kfill = jpfillcopy ) 390 380 ENDIF 391 381 !
Note: See TracChangeset
for help on using the changeset viewer.