Changeset 8127
- Timestamp:
- 2017-06-02T09:54:18+02:00 (8 years ago)
- Location:
- branches/2015/nemo_v3_6_STABLE/NEMOGCM
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/OPA_SRC/DIA/diadct.F90
r5505 r8127 249 249 !debug this section computing ? 250 250 lldebug=.FALSE. 251 IF( (jsec==nn_secdebug .OR. nn_secdebug==-1) .AND. kt==nit000+nn_dct-1 .AND. lwp) lldebug=.TRUE.251 IF( (jsec==nn_secdebug .OR. nn_secdebug==-1) .AND. kt==nit000+nn_dct-1 ) lldebug=.TRUE. 252 252 253 253 !Compute transport through section … … 258 258 IF( MOD(kt,nn_dctwri)==0 )THEN 259 259 260 IF( lwp .AND.kt==nit000+nn_dctwri-1 )WRITE(numout,*)" diadct: average transports and write at kt = ",kt260 IF( kt==nit000+nn_dctwri-1 )WRITE(numout,*)" diadct: average transports and write at kt = ",kt 261 261 262 262 !! divide arrays by nn_dctwri/nn_dct to obtain average … … 344 344 DO jsec=1,nb_sec_max !loop on the nb_sec sections 345 345 346 IF ( lwp .AND. ( jsec==nn_secdebug .OR. nn_secdebug==-1 )) &346 IF ( jsec==nn_secdebug .OR. nn_secdebug==-1 ) & 347 347 & WRITE(numout,*)'debuging for section number: ',jsec 348 348 … … 364 364 IF( jsec .NE. isec ) CALL ctl_stop( cltmp ) 365 365 366 IF( lwp .AND. ( jsec==nn_secdebug .OR. nn_secdebug==-1 ))WRITE(numout,*)"isec ",isec366 IF( jsec==nn_secdebug .OR. nn_secdebug==-1 )WRITE(numout,*)"isec ",isec 367 367 368 368 READ(numdct_in)secs(jsec)%name … … 383 383 !----- 384 384 385 IF( lwp .AND. ( jsec==nn_secdebug .OR. nn_secdebug==-1 ))THEN385 IF( jsec==nn_secdebug .OR. nn_secdebug==-1 )THEN 386 386 387 387 WRITE(clformat,'(a,i2,a)') '(A40,', nb_class_max,'(f8.3,1X))' … … 416 416 !debug 417 417 !----- 418 IF( lwp .AND. ( jsec==nn_secdebug .OR. nn_secdebug==-1 ))THEN418 IF( jsec==nn_secdebug .OR. nn_secdebug==-1 )THEN 419 419 WRITE(numout,*)" List of points in global domain:" 420 420 DO jpt=1,iptglo … … 450 450 !debug 451 451 !----- 452 IF( lwp .AND. ( jsec==nn_secdebug .OR. nn_secdebug==-1 ))THEN452 IF( jsec==nn_secdebug .OR. nn_secdebug==-1 )THEN 453 453 WRITE(numout,*)" List of points selected by the proc:" 454 454 DO jpt = 1,iptloc … … 468 468 !remove redundant points between processors 469 469 !------------------------------------------ 470 lldebug = .FALSE. ; IF ( (jsec==nn_secdebug .OR. nn_secdebug==-1) .AND. lwp) lldebug = .TRUE.470 lldebug = .FALSE. ; IF ( jsec==nn_secdebug .OR. nn_secdebug==-1 ) lldebug = .TRUE. 471 471 IF( iptloc .NE. 0 )THEN 472 472 CALL removepoints(secs(jsec),'I','top_list',lldebug) … … 484 484 !debug 485 485 !----- 486 IF( lwp .AND. ( jsec==nn_secdebug .OR. nn_secdebug==-1 ))THEN486 IF( jsec==nn_secdebug .OR. nn_secdebug==-1 )THEN 487 487 WRITE(numout,*)" List of points after removepoints:" 488 488 iptloc = secs(jsec)%nb_point … … 496 496 497 497 ELSE ! iptglo = 0 498 IF( lwp .AND. ( jsec==nn_secdebug .OR. nn_secdebug==-1 ))&498 IF( jsec==nn_secdebug .OR. nn_secdebug==-1 )& 499 499 WRITE(numout,*)' No points for this section.' 500 500 ENDIF -
branches/2015/nemo_v3_6_STABLE/NEMOGCM/TOOLS/SECTIONS_DIADCT/src/compute_sections.f90
r3632 r8127 681 681 IF( sec%nb_point .ne. 0 )THEN 682 682 IF ( sec%listPoint(sec%nb_point)%I .NE. sec%listPoint(1)%I ) THEN 683 sec%slopeSection = ( sec%listPoint(sec%nb_point)%J - sec%listPoint(1)%J ) / & 684 ( sec%listPoint(sec%nb_point)%I - sec%listPoint(1)%I ) 683 sec%slopeSection = REAL( ( sec%listPoint(sec%nb_point)%J - sec%listPoint(1)%J ) , wp )/ & 684 REAL( ( sec%listPoint(sec%nb_point)%I - sec%listPoint(1)%I ) , wp ) 685 685 686 ELSE 686 687 sec%slopeSection = 10000._wp
Note: See TracChangeset
for help on using the changeset viewer.