Ignore:
Timestamp:
Oct 11, 2001, 11:50:51 AM (23 years ago)
Author:
lmdzadmin
Message:

Mauvais dimensionnement de rlon, rlat
LF

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ.3.3/branches/rel-LF/libf/bibio/writephys.F90

    r272 r279  
    8989  character*(*) :: nom_fichier
    9090  integer :: iim, jjm, llm, klon
    91   real, dimension(iim) :: rlon
    92   real, dimension(jjm) :: rlat
     91  real, dimension(klon) :: rlon
     92  real, dimension(klon) :: rlat
    9393  real, dimension(llm) :: zlev
    9494  real                 :: date0, dtime
     
    274274  integer              :: file_id, isize, jsize, lsize, phy_lon, zsize
    275275  integer              :: nhori, nvert, klon
    276   integer              :: nitau
     276  integer              :: itau
    277277  real, dimension(:,:,:), allocatable :: temp_data
    278278  character (len = 20) :: modname = 'writephy'
     
    289289      hist_files(file)%define = .false.
    290290      call histend(hist_files(file)%file_id)
    291 !      hist_files(file)%nitau = hist_files(file)%nitau + 1
    292291    else
    293292!
     
    309308        nvert = hist_files(file)%nvert
    310309      endif
     310!      if (var_name == 'phis') then
     311!        write(*,*)'Define: ',var_name
     312!        write(*,*)file_id, var_name, var_title, var_units,  &
     313! &                isize, jsize, nhori, lsize, 1, lsize, nvert, nbits,  &
     314! &                var_op, freq_op, freq_wri
     315!      endif
    311316      call histdef(file_id, var_name, var_title, var_units,  &
    312317 &                isize, jsize, nhori, lsize, 1, lsize, nvert, nbits,  &
     
    328333  freq_op = var_profs(iprof, file)%freq_op
    329334  freq_wri = var_profs(iprof, file)%freq_wri
    330   nitau = hist_files(file)%nitau
     335  itau = hist_files(file)%nitau
    331336
    332337  if (var_profs(iprof, file)%zsize == 0) then
     
    344349!
    345350  call gr_fi_ecrit(lsize, klon, isize, jsize, data, temp_data)
    346   call histwrite(file_id, var_name,nitau,temp_data, &
     351  if (var_name == 'tsol') then
     352    write(*,*)'writephys itau = ', file_id, var_name, itau
     353!    write(*,*)file_id, var_name,itau,temp_data, &
     354! &               var_profs(iprof, file)%reg_size, &
     355! &               var_profs(iprof, file)%reg_index
     356  endif
     357  call histwrite(file_id, var_name,itau,temp_data, &
    347358 &               var_profs(iprof, file)%reg_size, &
    348359 &               var_profs(iprof, file)%reg_index)
    349    
     360
     361  deallocate(temp_data)
     362 
    350363  return
    351364!
Note: See TracChangeset for help on using the changeset viewer.