- Timestamp:
- 2016-04-21T18:15:17+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/IOM/iom.F90
r6487 r6491 710 710 CHARACTER(LEN=256) :: clname ! file name 711 711 CHARACTER(LEN=1) :: clrankpv, cldmspc ! 712 LOGICAL :: ll_depth_spec ! T => if kstart, kcount present then *only* use values for 3rd spatial dimension. 712 713 !--------------------------------------------------------------------- 713 714 ! … … 722 723 IF( PRESENT(kcount) .AND. (.NOT. PRESENT(kstart)) ) CALL ctl_stop(trim(clinfo), 'kcount present needs kstart present') 723 724 IF( PRESENT(kstart) .AND. (.NOT. PRESENT(kcount)) ) CALL ctl_stop(trim(clinfo), 'kstart present needs kcount present') 724 IF( PRESENT(kstart) .AND. idom /= jpdom_unknown ) CALL ctl_stop(trim(clinfo), 'kstart present needs kdom = jpdom_unknown') 725 IF( PRESENT(kstart) .AND. idom /= jpdom_unknown .AND. idom /= jpdom_autoglo_xy ) & 726 & CALL ctl_stop(trim(clinfo), 'kstart present needs kdom = jpdom_unknown or kdom = jpdom_autoglo_xy') 725 727 726 728 luse_jattr = .false. … … 755 757 ! update idom definition... 756 758 ! Identify the domain in case of jpdom_auto(glo/dta) definition 759 IF( idom == jpdom_autoglo_xy ) THEN 760 ll_depth_spec = .TRUE. 761 idom = jpdom_autoglo 762 ELSE 763 ll_depth_spec = .FALSE. 764 ENDIF 757 765 IF( idom == jpdom_autoglo .OR. idom == jpdom_autodta ) THEN 758 766 IF( idom == jpdom_autoglo ) THEN ; idom = jpdom_global … … 808 816 istart(idmspc+1) = itime 809 817 810 IF( PRESENT(kstart)) THEN ; istart(1:idmspc) = kstart(1:idmspc) ; icnt(1:idmspc) = kcount(1:idmspc)818 IF( PRESENT(kstart) .AND. .NOT. ll_depth_spec ) THEN ; istart(1:idmspc) = kstart(1:idmspc) ; icnt(1:idmspc) = kcount(1:idmspc) 811 819 ELSE 812 IF( idom == jpdom_unknown ) THEN ; icnt(1:idmspc) = idimsz(1:idmspc)820 IF( idom == jpdom_unknown ) THEN ; icnt(1:idmspc) = idimsz(1:idmspc) 813 821 ELSE 814 822 IF( .NOT. PRESENT(pv_r1d) ) THEN ! not a 1D array … … 833 841 ENDIF 834 842 IF( PRESENT(pv_r3d) ) THEN 835 IF( idom == jpdom_data ) THEN ; icnt(3) = jpkdta 836 ELSE ; icnt(3) = jpk 843 IF( idom == jpdom_data ) THEN ; icnt(3) = jpkdta 844 ELSE IF( ll_depth_spec .AND. PRESENT(kstart) ) THEN ; istart(3) = kstart(3); icnt(3) = kcount(3) 845 ELSE ; icnt(3) = jpk 837 846 ENDIF 838 847 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.