Changeset 1045 for LMDZ4


Ignore:
Timestamp:
Nov 6, 2008, 5:31:39 PM (16 years ago)
Author:
lmdzadmin
Message:

On sort meantaucld dans histISCCP, car sortie simulateur (enleve du histday)
Lecture fichiers input ISCCP "debut" physiq
IM

Location:
LMDZ4/trunk/libf/phylmd
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/trunk/libf/phylmd/calcul_simulISCCP.h

    r827 r1045  
    6060      CALL haut2bas(klon, klev, dem_c, dem_cH2B)
    6161      CALL haut2bas(klon, klevp1, paprs, phalf)
    62 c
    63 cIM lecture invtau, tautab des fichiers formattes
    64 c
    65       IF (debut) THEN
    66 c$OMP MASTER
    67 c     open(99,file='tautab.bin',access='sequential',
    68 c    $     form='unformatted',status='old')
    69 c     read(99) tautab_omp
    70 c
    71       open(99,file='tautab.formatted', FORM='FORMATTED')
    72       read(99,'(f30.20)') tautab_omp
    73       close(99)
    74 c
    75       open(99,file='invtau.formatted',form='FORMATTED')
    76       read(99,'(i10)') invtau_omp
    77       close(99)
    78 c$OMP END MASTER
    79 c$OMP BARRIER
    80       tautab=tautab_omp
    81       invtau=invtau_omp
    82 c
    83       ENDIF !debut
    8462c
    8563cIM: initialisation de seed

  • LMDZ4/trunk/libf/phylmd/ini_histISCCP.h

    r827 r1045  
    146146     .                "ave(X)", zcals(n),zout_isccp(n))
    147147c
     148         CALL histdef(nid_isccp, "meantaucld"//verticaxe(n),
     149     .                "ISCCP mean cloud optical thickness", "1",
     150     .                iim, jj_nb,nhori,1,1,1,-99,32,
     151     .                "ave(X)", zcals(n),zout_isccp(n))
     152c
    148153         ENDDO
    149154c
    … …  
    168173          CALL histdef(nid_isccp, "nsunlit"//verticaxe(n),
    169174     .                "Nb of calls with sunlit ", "%",
     175     .                iim, jj_nb,nhori,1,1,1,-99,32,
     176     .                "ave(X)", zcals(n),zout_isccp(n))
     177c
     178         CALL histdef(nid_isccp, "meantaucld"//verticaxe(n),
     179     .                "ISCCP mean cloud optical thickness", "1",
    170180     .                iim, jj_nb,nhori,1,1,1,-99,32,
    171181     .                "ave(X)", zcals(n),zout_isccp(n))

  • LMDZ4/trunk/libf/phylmd/phys_output_mod.F90

    r1043 r1045  
    194194  integer, dimension(nfiles) , save :: flag_u10_sol      = (/ 10, 4, 10, 10 /)
    195195  integer, dimension(nfiles) , save :: flag_v10_sol      = (/ 10, 4, 10, 10 /)
    196   integer, dimension(nfiles) , save :: flag_meantaucld   = (/ 10, 5, 10, 10 /)
    197196  integer, dimension(nfiles) , save :: flag_cldtau       = (/ 10, 5, 10, 10 /)                     
    198197  integer, dimension(nfiles) , save :: flag_cldemi       = (/ 10, 5, 10, 10 /)
    … …  
    614613 CALL histdef2d(iff,flag_weakinv, "weakinv","Weak inversion", "-")
    615614 CALL histdef2d(iff,flag_dthmin,"dthmin","dTheta mini", "K/m")
    616  if (ok_isccp) then
    617    CALL histdef2d(iff,flag_meantaucld,"meantaucld","ISCCP mean cloud optical thickness","1")
    618  endif
    619615 CALL histdef2d(iff,flag_rh2m,"rh2m","Relative humidity at 2m", "%" )
    620616 CALL histdef2d(iff,flag_qsat2m,"qsat2m","Saturant humidity at 2m", "%")

  • LMDZ4/trunk/libf/phylmd/phys_output_write.h

    r1043 r1045  
    664664      CALL histwrite_phy(nid_files(iff),"dthmin",itau_w,dthmin)
    665665      ENDIF
    666 
    667 #ifdef histISCCP
    668       IF(ok_isccp) THEN
    669        IF (MOD(itap,NINT(freq_ISCCP/dtime)).EQ.0) THEN
    670         IF (flag_meantaucld(iff)<=lev_files(iff)) THEN
    671        CALL histwrite_phy(nid_files(iff),"meantaucld",itau_w,meantaucld)
    672         ENDIF
    673        ENDIF
    674       ENDIF !(ok_isccp)
    675 #endif
    676666
    677667       IF (flag_cldtau(iff)<=lev_files(iff)) THEN

  • LMDZ4/trunk/libf/phylmd/physiq.F

    r1043 r1045  
    17531753c
    17541754      IF (MOD(itap-1,lmt_pas) .EQ. 0) THEN
    1755          WRITE(lunout,*)' PHYS cond  julien ',julien
     1755         if(prt_level.ge.1) WRITE(lunout,*)' PHYS cond  julien ',julien
    17561756         CALL ozonecm( FLOAT(julien), rlat, paprs, wo)
    17571757      ENDIF
    … …  
    20822082
    20832083cIM begin
    2084         print*,'physiq: cin pbase dnwd0 ftd fqd ',cin(1),pbase(1),
    2085      .dnwd0(1,1),ftd(1,1),fqd(1,1)
     2084c       print*,'physiq: cin pbase dnwd0 ftd fqd ',cin(1),pbase(1),
     2085c    .dnwd0(1,1),ftd(1,1),fqd(1,1)
    20862086cIM end
    20872087cIM cf. FH
    … …  
    26972697      IF (ok_isccp) THEN
    26982698c
     2699cIM lecture invtau, tautab des fichiers formattes
     2700c
     2701      IF (debut) THEN
     2702c$OMP MASTER
     2703c
     2704      open(99,file='tautab.formatted', FORM='FORMATTED')
     2705      read(99,'(f30.20)') tautab_omp
     2706      close(99)
     2707c
     2708      open(99,file='invtau.formatted',form='FORMATTED')
     2709      read(99,'(i10)') invtau_omp
     2710
     2711c     print*,'calcul_simulISCCP invtau_omp',invtau_omp
     2712c     write(6,'(a,8i10)') 'invtau_omp',(invtau_omp(i),i=1,100)
     2713
     2714      close(99)
     2715c$OMP END MASTER
     2716c$OMP BARRIER
     2717      tautab=tautab_omp
     2718      invtau=invtau_omp
     2719c
     2720      ENDIF !debut
    26992721c
    27002722cIM appel simulateur toutes les  NINT(freq_ISCCP/dtime) heures

  • LMDZ4/trunk/libf/phylmd/write_histISCCP.h

    r1038 r1045  
    2929        CALL histwrite_phy(nid_isccp,"nsunlit"//verticaxe(n),itau_w,
    3030     .                 nbsunlit(1,:,n))
     31c
     32       CALL histwrite_phy(nid_isccp,"meantaucld"//verticaxe(n),itau_w,
     33     .                 meantaucld(:,n))
    3134c
    3235        ENDDO ! n=1, napisccp
    … …  
    7174     .                 itau_w,nbsunlit(1,:,n))
    7275c
     76       CALL histwrite_phy(nid_isccp,"meantaucld"//verticaxe(n),itau_w,
     77     .                 meantaucld(:,n))
     78c
    7379        zx_tmp_fi2d(1:klon)=float(seed(1:klon,n))
    7480c
Note: See TracChangeset for help on using the changeset viewer.