Changeset 743 for LMDZ4/branches


Ignore:
Timestamp:
Dec 11, 2006, 4:55:31 PM (18 years ago)
Author:
Laurent Fairhead
Message:

Nettoyage sur la version YM
LF

Location:
LMDZ4/branches/V3_test/libf
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/branches/V3_test/libf/dyn3dpar/leapfrog_p.F

    r734 r743  
    418418     &                                jj_Nb_caldyn,0,0,TestRequest)
    419419 
     420        do j=1,nqmx
     421         call Register_SwapFieldHallo(q(1,1,j),q(1,1,j),ip1jmp1,llm,
     422     &                                jj_nb_caldyn,0,0,TestRequest)
     423        enddo
     424
    420425         call SetDistrib(jj_nb_caldyn)
    421426         call SendRequest(TestRequest)
    … …  
    478483        call WriteField_p('pkf',reshape(pkf,(/iip1,jmp1,llm/)))
    479484        call WriteField_p('phis',reshape(phis,(/iip1,jmp1/)))
    480 c        do j=1,nqmx
    481 c          call WriteField_p('q'//trim(int2str(j)),
    482 c     .                reshape(q(:,:,j),(/iip1,jmp1,llm/)))
    483 c        enddo       
     485        do j=1,nqmx
     486          call WriteField_p('q'//trim(int2str(j)),
     487     .                reshape(q(:,:,j),(/iip1,jmp1,llm/)))
     488        enddo       
    484489      endif
    485490 

  • LMDZ4/branches/V3_test/libf/dyn3dpar/times.F90

    r709 r743  
    126126        V2=timer_table_sqr(jj_nb,no_timer,mpi_rank)
    127127        V=timer_table(jj_nb,no_timer,mpi_rank)
    128         timer_delta(jj_nb,no_timer,mpi_rank)=sqrt((V2-V*V/N)/(N-1))
     128        timer_delta(jj_nb,no_timer,mpi_rank)=sqrt(abs(V2-V*V/N)/(N-1))
    129129      else
    130130        timer_delta(jj_nb,no_timer,mpi_rank)=0

  • LMDZ4/branches/V3_test/libf/phylmd/ini_histrac.h

    r726 r743  
    33!
    44#ifndef INCA
     5c$OMP MASTER
    56         CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
    67c
    … …  
    125126         ndex3d = 0
    126127         ndex = 0
     128c$OMP END MASTER
    127129#endif

  • LMDZ4/branches/V3_test/libf/phylmd/interface_surf.F90

    r735 r743  
    20652065! sorties NETCDF des champs recus
    20662066!
    2067 !ym       ndexcs(:)=0
    2068 !ym       itau_w = itau_phy + itime
    2069 !ym       CALL histwrite(nidcs,cl_read(1),itau_w,read_sst,iim*(jjm+1),ndexcs)
    2070 !ym       CALL histwrite(nidcs,cl_read(2),itau_w,read_sic,iim*(jjm+1),ndexcs)
    2071 !ym       CALL histwrite(nidcs,cl_read(3),itau_w,read_alb_sic,iim*(jjm+1),ndexcs)
    2072 !ym       CALL histwrite(nidcs,cl_read(4),itau_w,read_sit,iim*(jjm+1),ndexcs)
    2073 !ym       CALL histsync(nidcs)
     2067     if (monocpu) THEN
     2068       ndexcs(:)=0
     2069       itau_w = itau_phy + itime
     2070       CALL histwrite(nidcs,cl_read(1),itau_w,read_sst,iim*(jjm+1),ndexcs)
     2071       CALL histwrite(nidcs,cl_read(2),itau_w,read_sic,iim*(jjm+1),ndexcs)
     2072       CALL histwrite(nidcs,cl_read(3),itau_w,read_alb_sic,iim*(jjm+1),ndexcs)
     2073       CALL histwrite(nidcs,cl_read(4),itau_w,read_sit,iim*(jjm+1),ndexcs)
     2074       CALL histsync(nidcs)
     2075     endif
    20742076! pas utile      IF (npas-itime.LT.nexca )CALL histclo(nidcs)
    20752077
    … …  
    24012403      tmp_field(ig) = 1.
    24022404    enddo   
    2403     call WriteField_phy('knindex',tmp_field,1)
     2405!ym    call WriteField_phy('knindex',tmp_field,1)
    24042406   
    24052407    tmp_field=0.
    … …  
    24082410      tmp_field(ig) = tsurf_new(i)
    24092411    enddo   
    2410     call WriteField_phy('tsurf_new',tmp_field,1)
     2412!ym    call WriteField_phy('tsurf_new',tmp_field,1)
    24112413   
    24122414    tmp_field=0.
    … …  
    24152417      tmp_field(ig) = alb_new(i)
    24162418    enddo   
    2417     call WriteField_phy('alb_new',tmp_field,1)
     2419!ym    call WriteField_phy('alb_new',tmp_field,1)
    24182420   
    24192421!    tmp_field=0.
    … …  
    24222424!      tmp_field(ig) = pctsrf_new(i,nisurf)
    24232425 !   enddo   
    2424     call WriteField_phy('pctsrf_new', pctsrf_new(:,nisurf),1)
     2426!ym    call WriteField_phy('pctsrf_new', pctsrf_new(:,nisurf),1)
    24252427  endif
    24262428!ym  do j=1,jjphy_nb

  • LMDZ4/branches/V3_test/libf/phylmd/isccp_cloud_types.F

    r703 r743  
    3232     &)
    3333       
    34 !$Id$
    3534
    3635! Copyright Steve Klein and Mark Webb 2002 - all rights reserved.

  • LMDZ4/branches/V3_test/libf/phylmd/oasis.dummy

    r675 r743  
    9595      cl_writ(18)='CORIVFLU'
    9696      cl_writ(19)='COCALVIN'
    97 c$$$      cl_writ(13)='COZOTAUX'
    98 c$$$      cl_writ(14)='COZOTAUV'
    99 c$$$      cl_writ(15)='COMETAUY'
    100 c$$$      cl_writ(16)='COMETAUU'
     97cxxx      cl_writ(13)='COZOTAUX'
     98cxxx      cl_writ(14)='COZOTAUV'
     99cxxx      cl_writ(15)='COMETAUY'
     100cxxx      cl_writ(16)='COMETAUU'
    101101      cl_writ(1)='COTAUXXU'
    102102      cl_writ(2)='COTAUYYU'
    … …  
    365365      REAL rivrunoff(imjm)
    366366      REAL calving(imjm)
    367 c$$$      REAL tauxu(imjm)
    368 c$$$      REAL tauxv(imjm)
    369 c$$$      REAL tauyu(imjm)
    370 c$$$      REAL tauyv(imjm)
     367cxxx      REAL tauxu(imjm)
     368cxxx      REAL tauxv(imjm)
     369cxxx      REAL tauyu(imjm)
     370cxxx      REAL tauyv(imjm)
    371371      REAL tauxx_u(imjm)
    372372      REAL tauxx_v(imjm)
    … …  
    496496     $          CALL locwrite(cl_writ(jf),calving, imjm,
    497497     $          file_unit_field(jf), ierror)
    498 c$$$            IF (jf.eq.13)
    499 c$$$     $          CALL locwrite(cl_writ(jf),tauxu, imjm,
    500 c$$$     $          file_unit_field(jf),ierror)
    501 c$$$            IF (jf.eq.1')
    502 c$$$     $          CALL locwrite(cl_writ(jf),tauxv, imjm,
    503 c$$$     $          file_unit_field(jf),ierror)
    504 c$$$            IF (jf.eq.15)
    505 c$$$     $          CALL locwrite(cl_writ(jf),tauyv, imjm,
    506 c$$$     $          file_unit_field(jf),ierror)
    507 c$$$            IF (jf.eq.16)
    508 c$$$     $          CALL locwrite(cl_writ(jf),tauyu, imjm,
    509 c$$$     $          file_unit_field(jf), ierror)
     498cxxx            IF (jf.eq.13)
     499cxxx     $          CALL locwrite(cl_writ(jf),tauxu, imjm,
     500cxxx     $          file_unit_field(jf),ierror)
     501cxxx            IF (jf.eq.1')
     502cxxx     $          CALL locwrite(cl_writ(jf),tauxv, imjm,
     503cxxx     $          file_unit_field(jf),ierror)
     504cxxx            IF (jf.eq.15)
     505cxxx     $          CALL locwrite(cl_writ(jf),tauyv, imjm,
     506cxxx     $          file_unit_field(jf),ierror)
     507cxxx            IF (jf.eq.16)
     508cxxx     $          CALL locwrite(cl_writ(jf),tauyu, imjm,
     509cxxx     $          file_unit_field(jf), ierror)
    510510            IF (jf.eq.1)
    511511     $          CALL locwrite(cl_writ(jf),tauxx_u, imjm,
    … …  
    574574          IF (jn.eq.18) CALL CLIM_Export(cl_writ(jn),kt,rivrunoff,info)
    575575          IF (jn.eq.19) CALL CLIM_Export(cl_writ(jn),kt,calving,info)
    576 c$$$          IF (jn.eq.13) CALL CLIM_Export(cl_writ(jn), kt, tauxu, info)
    577 c$$$          IF (jn.eq.14) CALL CLIM_Export(cl_writ(jn), kt, tauxv, info)
    578 c$$$          IF (jn.eq.15) CALL CLIM_Export(cl_writ(jn), kt, tauyv, info)
    579 c$$$          IF (jn.eq.16) CALL CLIM_Export(cl_writ(jn), kt, tauyu, info)
     576cxxx          IF (jn.eq.13) CALL CLIM_Export(cl_writ(jn), kt, tauxu, info)
     577cxxx          IF (jn.eq.14) CALL CLIM_Export(cl_writ(jn), kt, tauxv, info)
     578cxxx          IF (jn.eq.15) CALL CLIM_Export(cl_writ(jn), kt, tauyv, info)
     579cxxx          IF (jn.eq.16) CALL CLIM_Export(cl_writ(jn), kt, tauyu, info)
    580580          IF (jn.eq.1) CALL CLIM_Export(cl_writ(jn), kt, tauxx_u, info)
    581581          IF (jn.eq.2) CALL CLIM_Export(cl_writ(jn), kt, tauyy_u, info)
Note: See TracChangeset for help on using the changeset viewer.