Changeset 2912
- Timestamp:
- 2011-10-13T15:40:45+02:00 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/dev_r2802_MERCATOR10_diadct/NEMOGCM/NEMO/OPA_SRC/DIA/diadct.F90
r2909 r2912 277 277 READ(numdct_in,iostat=iost)isec 278 278 IF (iost .NE. 0 )EXIT !end of file 279 WRITE(cltmp,'(a,i4.4,a,i4.4)')'diadct: read sections : Problem of section number: isec= ',isec,' etjsec= ',jsec279 WRITE(cltmp,'(a,i4.4,a,i4.4)')'diadct: read sections : Problem of section number: isec= ',isec,' and jsec= ',jsec 280 280 IF( jsec .NE. isec ) CALL ctl_stop( cltmp ) 281 281 … … 306 306 WRITE(numout,*) " Number of points in the section: ",iptglo 307 307 WRITE(numout,*) " Number of classes ",secs(jsec)%nb_class 308 WRITE(numout,'(A40,10(f8.3,1X))')" Insitu density classes : ",secs(jsec)%zsigi 308 ! WRITE(numout,'(A40,10(f8.3,1X))')" Insitu density classes : ",secs(jsec)%zsigi 309 WRITE(numout,'(A40,nb_class_max(f8.3,1X))')" Insitu density classes : ",secs(jsec)%zsigi 309 310 WRITE(numout,'(A40,10(f8.3,1X))')" Potential density classes : ",secs(jsec)%zsigp 310 311 WRITE(numout,'(A40,10(f8.3,1X))')" Salinity classes : ",secs(jsec)%zsal … … 353 354 iptloc = iptloc + 1 ! count local points 354 355 secs(jsec)%listPoint(iptloc) = POINT_SECTION(mi0(iiglo),mj0(ijglo)) ! store local coordinates 355 secs(jsec)%direction(iptloc) = directemp(jpt) ! store local direction356 secs(jsec)%direction(iptloc) = directemp(jpt) ! store local direction 356 357 ENDIF 357 358 … … 362 363 !debug 363 364 !----- 364 !IF( lwp .AND. ( jsec==nn_secdebug .OR. nn_secdebug==-1 ) )THEN 365 IF( ( jsec==nn_secdebug .OR. nn_secdebug==-1 ) )THEN 365 IF( lwp .AND. ( jsec==nn_secdebug .OR. nn_secdebug==-1 ) )THEN 366 366 WRITE(numout,*)" List of points selected by the proc:" 367 367 DO jpt = 1,iptloc … … 369 369 ijglo = secs(jsec)%listPoint(jpt)%J + jpjzoom - 1 + njmpp - 1 370 370 WRITE(numout,*)' # I J : ',iiglo,ijglo 371 !WRITE(narea+200,*)' # I J : ',iiglo,ijglo372 371 ENDDO 373 372 ENDIF … … 396 395 ELSE ! iptglo = 0 397 396 IF( lwp .AND. ( jsec==nn_secdebug .OR. nn_secdebug==-1 ) )& 398 WRITE(numout,*)' No tpoints for this section.'397 WRITE(numout,*)' No points for this section.' 399 398 ENDIF 400 399 … … 425 424 itest ,& !indice value of the side of the domain 426 425 !where points could be redundant 427 isgn ,& !way of course in listpoint 428 ipoint ,& !way of course in listpoint 426 isgn ,& ! isgn= 1 : scan listpoint from start to end 427 ! isgn=-1 : scan listpoint from end to start 428 ipoint ,& !way of course in listpoint 429 429 istart,iend !first and last points selected in listpoint 430 430 INTEGER :: jpoint =0 !loop on list points 431 INTEGER,DIMENSION(nb_point_max) :: idirec !contains temporar esec%direction432 INTEGER,DIMENSION(2,nb_point_max) :: icoord !contains temporar esec%listpoint431 INTEGER,DIMENSION(nb_point_max) :: idirec !contains temporary sec%direction 432 INTEGER,DIMENSION(2,nb_point_max) :: icoord !contains temporary sec%listpoint 433 433 !---------------------------------------------------------------------------- 434 434 IF( ld_debug )WRITE(numout,*)' -------------------------' … … 482 482 483 483 SUBROUTINE transport(sec,ld_debug) 484 !!--------------------------------------------------------------------- 484 !!------------------------------------------------------------------------------------------- 485 485 !! *** ROUTINE transport *** 486 486 !! … … 497 497 !! test on the density/temperature/salinity/level 498 498 !! 499 !! ** Output: sec%transport: transport in the 2 direction and temperature, 500 !! salinity, density meaned by the transport 501 !! 502 !! 503 !!------------------------------------------------------------------- 499 !! ** Output: sec%transport: volume/mass/ice/heat/salt transport in the 2 directions 500 !! 501 !! 502 !!------------------------------------------------------------------------------------------- 504 503 !! * Arguments 505 504 TYPE(SECTION),INTENT(INOUT) :: sec … … 513 512 zumid_ice , zvmid_ice ,&!U/V ice velocity 514 513 zTnorm ,&!transport of velocity through one cell's sides 515 ztransp1 , ztransp2 ,&!total transport in sens 1 and 2514 ztransp1 , ztransp2 ,&!total transport in directions 1 and 2 516 515 ztemp1 , ztemp2 ,&!temperature transport " 517 516 zrhoi1 , zrhoi2 ,&!mass transport " 518 517 zrhop1 , zrhop2 ,&!mass transport " 519 518 zsal1 , zsal2 ,&!salinity transport " 520 zice_vol_pos , zice_vol_neg ,&!volum ic icetransport "521 zice_surf_pos, zice_surf_neg !surfac ic icetransport "519 zice_vol_pos , zice_vol_neg ,&!volume ice transport " 520 zice_surf_pos, zice_surf_neg !surface ice transport " 522 521 REAL(wp):: ztn, zsn, zrhoi, zrhop, zsshn, zfsdep ! temperature/salinity/ssh/potential density /depth at u/v point 523 522 REAL(wp) :: aj,i0,j0,i1,j1,i,j … … 545 544 ! 546 545 !convention: 547 ! non horizontal section: sens+ is toward left hand of section548 ! horizontal section: sens+ is toward north of section546 ! non horizontal section: direction + is toward left hand of section 547 ! horizontal section: direction + is toward north of section 549 548 ! 550 549 ! … … 552 551 ! ---------------- ----------------- --------------- -------------- 553 552 ! 554 ! isgnv=1 sens+555 ! ______ _____ ______ (i1,j1) (i1,j1)556 ! (i0,j0) | //| | | sens+557 ! | isgnu=1 // | |isgnu=1 |isgnu=1 /|\558 ! |_______ // ______| \\ | ---\ |559 ! | | isgnv=-1 \\ | | ---/ sens + ____________(i1,j1)560 ! | | __\\| | (i0,j0)561 ! | (i1,j1) |(i0,j0) sens+ | isgnv=1562 ! (i0,j0)553 ! isgnv=1 direction + 554 ! ______ _____ ______ 555 ! | //| | | direction + 556 ! | isgnu=1 // | |isgnu=1 |isgnu=1 /|\ 557 ! |_______ // ______| \\ | ---\ | 558 ! | | isgnv=-1 \\ | | ---/ direction + ____________ 559 ! | | __\\| | 560 ! | | direction + | isgnv=1 561 ! 563 562 !---------------------------------------------------------------------------------------------------- 564 563 isgnu = 1 … … 576 575 !------------------------------------------ 577 576 ! Select good coordinate to have velocity of the segment 577 ! 578 ! CASE( sec%direction(jseg) : 578 579 ! 579 580 ! |U CASE(3) … … 647 648 648 649 649 !compute velocity with the good sens650 !compute velocity with the correct direction 650 651 SELECT CASE( sec%direction(jseg) ) 651 652 CASE(0,1) … … 657 658 END SELECT 658 659 659 !velocity* cell's leng ht* cell's thickness660 !velocity* cell's length * cell's thickness 660 661 zTnorm=zumid*e2u(k%I,k%J)* fse3u(k%I,k%J,jk)+ & 661 662 zvmid*e1v(k%I,k%J)* fse3v(k%I,k%J,jk) … … 769 770 770 771 771 ELSE !if nb_inmesh=0772 ELSE !if sec%nb_point =0 772 773 zsum(1:2,:)=0. 773 774 IF (sec%llstrpond) zsum(3:10,:)=0. 774 775 zsum( 11:14,:)=0. 775 ENDIF !end of nb_inmesh=0 case776 ENDIF !end of sec%nb_point =0 case 776 777 777 778 !-------------------------------|
Note: See TracChangeset
for help on using the changeset viewer.