Changeset 2864
- Timestamp:
- 2011-09-26T17:02:11+02:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/dev_r2802_MERCATOR10_diadct/NEMOGCM/NEMO/OPA_SRC/DIA/diadct.F90
r2857 r2864 335 335 secs(jsec)%listPoint(iptloc) = POINT_SECTION(mi0(iiglo),mj0(ijglo)) ! store local coordinates 336 336 secs(jsec)%direction(iptloc) = directemp(jpt) ! store local direction 337 !WRITE(narea+200,*)' # I J : ',iiloc,ijloc 337 338 ENDIF 338 339 … … 343 344 !debug 344 345 !----- 345 IF( lwp .AND. ( jsec==nn_secdebug .OR. nn_secdebug==-1 ) )THEN 346 !IF( lwp .AND. ( jsec==nn_secdebug .OR. nn_secdebug==-1 ) )THEN 347 IF( ( jsec==nn_secdebug .OR. nn_secdebug==-1 ) )THEN 346 348 WRITE(numout,*)" List of points selected by the proc:" 347 349 DO jpt = 1,iptloc … … 349 351 ijglo = secs(jsec)%listPoint(jpt)%J + jpjzoom - 1 + njmpp - 1 350 352 WRITE(numout,*)' # I J : ',iiglo,ijglo 353 !WRITE(narea+200,*)' # I J : ',iiglo,ijglo 351 354 ENDDO 352 355 ENDIF … … 582 585 CASE(3) ; k = POINT_SECTION(sec%listPoint(jseg)%I,sec%listPoint(jseg)%J+1) 583 586 END SELECT 584 587 585 588 !------------------------------- 586 589 ! LOOP ON THE DENSITY CLASSES | … … 620 623 !TEST ON THE DENSITY/SALINITY/TEMPERATURE/LEVEL! 621 624 !----------------------------------------------! 625 IF(ld_debug)write(narea+200,*)k ,jk 622 626 623 627 IF ( ( ((( zrhop .GE. (sec%zsigp(jclass)+1000. )) .AND. & 624 628 ( zrhop .LE. (sec%zsigp(jclass+1)+1000. ))) .OR. & 625 629 ( sec%zsigp(jclass) .EQ. 99.)) .AND. & 626 ((( zrhoi .GE. (sec%zsigi(jclass) )) .AND. &627 ( zrhoi .LE. (sec%zsigi(jclass+1) ))) .OR. &630 ((( zrhoi .GE. (sec%zsigi(jclass) + 1000. )) .AND. & 631 ( zrhoi .LE. (sec%zsigi(jclass+1)+1000. ))) .OR. & 628 632 ( sec%zsigi(jclass) .EQ. 99.)) .AND. & 629 633 ((( zsn .GT. sec%zsal(jclass)) .AND. & … … 659 663 ENDIF 660 664 #endif 665 !IF(ld_debug)write(narea+200,*)k ,jk,zTnorm 661 666 !COMPUTE TRANSPORT 662 667 !zTnorm=transport through one cell for one class … … 1054 1059 !!---------------------------------------------------------------------- 1055 1060 LOGICAL, PUBLIC, PARAMETER :: lk_diadct = .FALSE. !: diamht flag 1061 PUBLIC 1056 1062 CONTAINS 1063 1064 SUBROUTINE dia_dct_init ! Dummy routine 1065 WRITE(*,*) 'dia_dct_init: You should not have seen this print! error?', kt 1066 END SUBROUTINE dia_dct_init 1067 1057 1068 SUBROUTINE dia_dct( kt ) ! Dummy routine 1058 1069 INTEGER, INTENT( in ) :: kt ! ocean time-step index
Note: See TracChangeset
for help on using the changeset viewer.