Changeset 489 for trunk/NEMO/TOP_SRC/trcdit.F90
- Timestamp:
- 2006-09-01T15:48:40+02:00 (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/TOP_SRC/trcdit.F90
r433 r489 146 146 147 147 ! define time axis 148 it = kt - nit 000 + 1148 it = kt - nittrc000 + 1 149 149 150 150 ! 1. Define NETCDF files and fields at beginning of first time step … … 152 152 153 153 IF(ll_print)WRITE(numout,*)'trcdit_wr kt=',kt,' kindic ',kindic 154 IF(kt == nit 000) THEN154 IF(kt == nittrc000) THEN 155 155 156 156 ! Compute julian date from starting date of the run … … 158 158 CALL ymds2ju(nyear,nmonth,nday,0.0,zjulian) 159 159 IF(lwp)WRITE(numout,*)' ' 160 IF(lwp)WRITE(numout,*)' Date 0 used :',nit 000 &160 IF(lwp)WRITE(numout,*)' Date 0 used :',nittrc000 & 161 161 & ,' YEAR ',nyear,' MONTH ',nmonth,' DAY ',nday & 162 162 & ,'Julian day : ',zjulian … … 177 177 ! Vertical grid for tracer : gdept 178 178 CALL histvert(nit5, 'deptht', 'Vertical T levels', & 179 & 'm', ipk, gdept , ndepit5)179 & 'm', ipk, gdept_0, ndepit5) 180 180 181 181 ! Index of ocean points in 3D and 2D (surface) … … 338 338 339 339 ! define time axis 340 it = kt - nit 000 + 1340 it = kt - nittrc000 + 1 341 341 342 342 ! Define the NETCDF files (one per tracer) 343 343 ! 344 344 IF(ll_print)WRITE(numout,*)'trcdid kt=',kt,' kindic ',kindic 345 IF(kt == nit 000) THEN345 IF(kt == nittrc000) THEN 346 346 347 347 DO jn=1,jptra … … 366 366 ! Vertical grid for tracer trend - one per each tracer IF needed 367 367 CALL histvert(nit6(jn), 'deptht', 'Vertical T levels', & 368 & 'm', ipk, gdept , ndepit6(jn))368 & 'm', ipk, gdept_0, ndepit6(jn)) 369 369 370 370 … … 647 647 648 648 ! define time axis 649 it = kt - nit 000 + 1649 it = kt - nittrc000 + 1 650 650 651 651 ! 1. Define NETCDF files and fields at beginning of first time step … … 653 653 654 654 IF(ll_print)WRITE(numout,*)'trcdii_wr kt=',kt,' kindic ',kindic 655 IF(kt == nit 000) THEN655 IF(kt == nittrc000) THEN 656 656 657 657 ! Define the NETCDF files for additional arrays : 2D or 3D … … 671 671 672 672 CALL histvert(nitd, 'deptht', 'Vertical T levels', & 673 & 'm', ipk, gdept , ndepitd)673 & 'm', ipk, gdept_0, ndepitd) 674 674 675 675 … … 729 729 DO jn=1,jpdia2d 730 730 cltra=ctrc2d(jn) ! short title for 2D diagnostic 731 CALL histwrite(nitd, cltra, kt, trc2d(:,:,jn), ndimt51 &731 CALL histwrite(nitd, cltra, it, trc2d(:,:,jn), ndimt51 & 732 732 & ,ndext51) 733 733 END DO … … 849 849 850 850 ! define time axis 851 it = kt - nit 000 + 1851 it = kt - nittrc000 + 1 852 852 853 853 ! 1. Define NETCDF files and fields at beginning of first time step … … 855 855 856 856 IF(ll_print)WRITE(numout,*)'trcdib_wr kt=',kt,' kindic ',kindic 857 IF(kt == nit 000) THEN857 IF(kt == nittrc000) THEN 858 858 859 859 ! Define the NETCDF files for biological trends … … 868 868 ! Vertical grid for biological trends 869 869 CALL histvert(nitb, 'deptht', 'Vertical T levels', & 870 & 'm', ipk, gdept , ndepitb)870 & 'm', ipk, gdept_0, ndepitb) 871 871 872 872 ! Declare all the output fields as NETCDF variables … … 905 905 DO jn=1,jpdiabio 906 906 cltra=ctrbio(jn) ! short title for biological diagnostic 907 CALL histwrite(nitb, cltra, kt, trbio(:,:,:,jn), ndimt50,ndext50)907 CALL histwrite(nitb, cltra, it, trbio(:,:,:,jn), ndimt50,ndext50) 908 908 END DO 909 909
Note: See TracChangeset
for help on using the changeset viewer.