Changeset 6 in lmdz_wrf


Ignore:
Timestamp:
Jul 18, 2014, 11:20:11 AM (10 years ago)
Author:
lfita
Message:

Fixing problem with the solar radiation:

1.- Correct values for the julian and the hour (0,1) to the LMDZ physiq

Adding NaN checking for 'PSFC' and 'T_2'

Location:
WRFV3
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • WRFV3/dyn_em/solve_em.F

    r1 r6  
    186186   CHARACTER(LEN=256)                         :: mminlu, mminsl
    187187   CHARACTER(LEN=50)                          :: errmsg
     188   INTEGER                                    :: hr, minute, sec, ms, julyr, &
     189     julday
     190   REAL                                       :: gmt
     191
    188192
    189193   errmsg = 'ERROR -- error -- ERROR -- error'
    … …  
    845849
    846850       grid%qv_2 = moist(:,:,:,P_QV)
     851
     852! L. Fita, LMD. July 2014. Getting hour of the day
     853       CALL domain_clock_get( grid, current_time=CurrTime )
     854       CALL WRFU_TimeGet( CurrTime, YY= julyr, dayOfYear=julday, H=hr, M=minute, S=sec, MS=ms, rc=rc)
     855! Julian day hour (0, 1) !!
     856       gmt=(hr+real(minute)/60.+real(sec)/3600.+real(ms)/(1000*3600))/24.
     857
     858! Checking for NaNs (should not be necessary but....)
     859   im2 = ims + (ime - ims) / 2
     860   jm2 = jms + (jme - jms) / 2
     861
     862   IF (grid%t_2(im2,iz,jm2) /= grid%t_2(im2,iz,jm2) .OR. ABS(grid%t_2(im2,iz,jm2)) > 10000. ) THEN
     863     PRINT errmsg
     864     WRITE(wrF_err_message,*)'solve_em: wrong T value=',                 &
     865       grid%t_2(im2,iz,jm2),' at: ', im2,', ', iz,', ', jm2,' !!!'
     866#ifdef DM_PARALLEL
     867     CALL wrf_error_fatal(TRIM(wrf_err_message))
     868#else
     869     PRINT *,TRIM(wrf_err_message)
     870     STOP
     871#endif
     872   END IF
     873! Checking for NaNs (should not be necessary but....)
     874   IF (grid%psfc(im2,jm2) /= grid%psfc(im2,jm2) .OR. ABS(grid%psfc(im2,jm2)) > 1000000. ) THEN
     875     PRINT errmsg
     876     WRITE(wrF_err_message,*)'solve_em: wrong PSFC value=',               &
     877       grid%psfc(im2,jm2),' at: ', im2 ,', ', jm2, ' !!!'
     878#ifdef DM_PARALLEL
     879     CALL wrf_error_fatal(TRIM(wrf_err_message))
     880#else
     881     PRINT *,TRIM(wrf_err_message)
     882     STOP
     883#endif
     884   END IF
     885
    847886       IF (config_flags%lmdz_physics) THEN
    848887
    … …  
    865904        &             ,WNUM3DM=num_3d_m, WPARFIRSTSCAL=PARAM_FIRST_SCALAR,           &
    866905        &        WNX=config_flags%e_we, WNY=config_flags%e_sn,                       &
    867         &        WNZ=config_flags%e_vert, WJULDAY=grid%julday, WGMT=grid%gmt,        &
     906        &        WNZ=config_flags%e_vert, WJULDAY=FLOAT(julday), WGMT=gmt,           &
    868907        &        WTIME_STEP=REAL(config_flags%time_step),                            &
    869908        &        WRF_FULLETA=grid%znw, WRF_HALFETA=grid%znu, WRF_DFULLETA=grid%dnw,  &

  • WRFV3/lmdz/physiq.F90

    r1 r6  
    13071307!
    13081308      CALL phys_cal_update(jD_cur,jH_cur)
     1309      PRINT *,'  Lluis physiq: jD_cur: ',jD_cur, ' jH_cur: ',jH_cur, &
     1310        ' days_elapsed: ',days_elapsed
    13091311
    13101312!c======================================================================
    … …  
    13551357      print*, '================================================='
    13561358      print*, 'Allocation des variables locales et sauvegardees'
    1357 !      PRINT *,'  Lluis before phys_local_var_init ftsol: ',ftsol(llp,:)
    13581359
    13591360      call phys_local_var_init
    1360       PRINT *,'  Lluis after phys_local_var_init ftsol: ',ftsol(llp,:)
    13611361
    13621362!c
    … …  
    13781378
    13791379
    1380       PRINT *,'  Lluis after conf_phys qsol: ',qsol(llp),      &
    1381         ' ftsol: ',ftsol(llp,:)
    1382 
    13831380!!      call phys_state_var_init(read_climoz)
    13841381      call phys_output_var_init
    13851382
    1386       PRINT *,'  Lluis after phys_output_var_init qsol: ',qsol(llp),      &
    1387         ' ftsol: ',ftsol(llp,:)
    13881383
    13891384!L. Fita, LMD. November 2013
    … …  
    13941389!!      CALL limit_initialize()
    13951390!!      CALL vars_limit_init(klon)
    1396       PRINT *,'  Lluis before all rst:', ALLOCATED(qsurf_rst), ALLOCATED(qsol_rst)
    1397 !      PRINT *,'  Lluis before all :', ALLOCATED(qsurf), ALLOCATED(qsol)
     1391
    13981392! L. Fita, LMD January 2014. Initializing output variables
    13991393      CALL init_ovars_lmdz_NOmodule(klon,klev,nbsrf)
    1400 
    1401       PRINT *,'  Lluis after init_ovas_lmdz_NOmodule qsol: ',qsol(llp),      &
    1402         ' ftsol: ',ftsol(llp,:)
    14031394
    14041395! Allocating restart variables
    … …  
    14131404!!      ALLOCATE(restart_runoff(klon))
    14141405
    1415 !!      PRINT *,'  Lluis physiq before pbl_surface_init 315 qsurf: ', qsurf(llp,:),    &
    1416 !!        ' dims: ',UBOUND(qsurf),' allocated: ', ALLOCATED(qsurf),' qsol: ',qsol(llp)
    14171406!      CALL pbl_surface_init(qsol_rst, fder_rst, snow_rst, qsurf_rst, evap_rst,       &
    14181407!        rugos_rst, agesno_rst, ftsoil_rst)
    1419 !!      PRINT *,'  Lluis physiq after pbl_surface_init 315 qsurf: ', qsurf(llp,:),     &
    1420 !!        ' dims: ',UBOUND(qsurf),' allocated: ', ALLOCATED(qsurf),' qsol: ',qsol(llp)
    14211408!      CALL fonte_neige_init(restart_runoff)
    14221409!      qsol=qsol_rst
    … …  
    14641451
    14651452       IF (debut) THEN
    1466          PRINT *,'Lluis debut!!!!!'
    14671453!rv
    14681454!cCRinitialisation de wght_th et lalim_conv pour la definition de la couche alimentation
    … …  
    15371523         radpas = NINT( 86400./dtime/nbapp_rad)
    15381524!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    1539          PRINT *,'Lluis before phyetat0'
    15401525
    15411526         CALL phyetat0 ("startphy.nc",clesphy0,tabcntr0)
    1542         PRINT *,'  Lluis after phyetat0 qsol: ',qsol(llp),      &
    1543         ' ftsol: ',ftsol(llp,:)
    1544 
    1545          PRINT *,'Lluis after phyetat0'
     1527
    15461528         IF (klon_glo==1) THEN
    15471529         coefh=0. ; coefm=0. ; pbl_tke=0.
    … …  
    16901672             rugoro(i) = f_rugoro * MAX(1.0e-05, zstd(i)*zsig(i)/2.0)
    16911673           ENDDO
    1692 
    1693       PRINT *,'  Lluis after LMDZ_CFMIP_point_locations qsol: ',qsol(llp),      &
    1694         ' ftsol: ',ftsol(llp,:)
    1695 
    16961674
    16971675!c34EK
    … …  
    17431721#ifdef CPP_IOIPSL
    17441722
    1745       PRINT *,'  Lluis before phys_output_open qsol: ',qsol(llp),      &
    1746         ' ftsol: ',ftsol(llp,:)
    1747 
    17481723!$OMP MASTER
    17491724       call phys_output_open(rlon,rlat,nCFMIP,tabijGCM,                              &
    … …  
    17581733!$OMP END MASTER
    17591734!$OMP BARRIER
    1760 
    1761       PRINT *,'  Lluis after phys_output_open qsol: ',qsol(llp),      &
    1762         ' ftsol: ',ftsol(llp,:)
    1763 
    17641735
    17651736#ifdef histISCCP
    … …  
    18311802
    18321803      call iniradia(klon,klev,paprs(1,1:klev+1))
    1833 
    1834       PRINT *,'  Lluis after iniradia qsol: ',qsol(llp),      &
    1835         ' ftsol: ',ftsol(llp,:)
    18361804
    18371805!$OMP single
    … …  
    18771845!
    18781846!
    1879       PRINT *,'  Lluis end of -debut- !!!!'
    18801847! Incrementer le compteur de la physique
    18811848!
    18821849      itap   = itap + 1
    1883       PRINT *,'  Lluis 5 11 itap: ',itap,' evap_srf: ',evap_srf(llp,:)
    18841850!c
    18851851!
    … …  
    18881854! on the surface fraction.
    18891855!
    1890       PRINT *,'  Lluis before change_srf_frac qsol: ',qsol(llp),    &
    1891         ' ftsol: ',ftsol(llp,:)
    18921856      CALL change_srf_frac(itap, dtime, days_elapsed+1,                              &
    18931857       &     pctsrf, falb1, falb2, ftsol, ustar, u10m, v10m, pbl_tke)
    1894       PRINT *,'  Lluis after change_srf_frac qsol: ',qsol(llp),    &
    1895         ' ftsol: ',ftsol(llp,:)
    18961858
    18971859! Update time and other variables in Reprobus
    … …  
    19491911!c Ne pas affecter les valeurs entrees de u, v, h, et q
    19501912!c
    1951       PRINT *,'  Lluis 550 k t_seri t_____'
    19521913      DO k = 1, klev
    19531914      DO i = 1, klon
    … …  
    19591920         qs_seri(i,k) = 0.
    19601921      ENDDO
    1961          PRINT *,k,t_seri(550,k), t(550,k)
    19621922      ENDDO
    19631923      tke0(:,:)=pbl_tke(:,:,is_ave)
    … …  
    20021962       &      , fs_bound, fq_bound )
    20031963      END IF
    2004 
    2005       PRINT *,'  Lluis after diageptq/diagphy qsol: ',qsol(llp),      &
    2006         ' ftsol: ',ftsol(llp,:)
    20071964
    20081965!c Diagnostiquer la tendance dynamique
    … …  
    20752032       PRINT*,' debut avant hgardfou min max ftsol',itap,amn,amx
    20762033      ENDIF !(check) THEN
    2077       PRINT *,'  Lluis before devu hgardfou qsol: ',qsol(llp),      &
    2078         ' ftsol: ',ftsol(llp,:)
    2079       PRINT *,'   t_seri: ',t_seri(550,:)
    20802034
    20812035!IM END
    … …  
    20832037
    20842038      CALL hgardfou(t_seri,ftsol,'debutphy')
    2085       PRINT *,'  Lluis after devu hgardfou qsol: ',qsol(llp),      &
    2086         ' ftsol: ',ftsol(llp,:)
    20872039!c
    20882040!IM BEG
    … …  
    20972049       PRINT*,' debut apres hgardfou min max ftsol',itap,amn,amx
    20982050      ENDIF !(check) THEN
    2099       PRINT *,'  Lluis after check qsol: ',qsol(llp),      &
    2100         ' ftsol: ',ftsol(llp,:)
    21012051
    21022052!IM END
    … …  
    21392089         wo(:, :, 1) = ozonecm(rlat, paprs, rjour=real(days_elapsed+1))
    21402090      ENDIF
    2141       PRINT *,'  Lluis after initial conditions qsol: ',qsol(llp),      &
    2142         ' ftsol: ',ftsol(llp,:)
    21432091
    21442092!c
    … …  
    22002148      endif
    22012149      if(prt_level.ge.1)                                                &
    2202        &    write(lunout,*)'Longitude solaire ',zlongi,solarlong0,dist
    2203       PRINT *,'  Lluis after solar orbit qsol: ',qsol(llp),      &
    2204         ' ftsol: ',ftsol(llp,:)
     2150       &    write(lunout,*)'Longitude solaire ',zlongi,solarlong0,dist,&
     2151       ' jD_cur: ',jD_cur,' jH_cur: ',jH_cur,' days_elapsed: ', &
     2152       REAL(days_elapsed+1)
    22052153
    22062154
    … …  
    22132161! non nul aux poles.
    22142162      IF (abs(solarlong0-1000.)<1.e-4) then
     2163         PRINT *,' Lluis calling zenang_an'
    22152164         call zenang_an(cycle_diurne,jH_cur,rlat,rlon,rmu0,fract)
    22162165      ELSE
    … …  
    22192168           zdtime=dtime*REAL(radpas) ! pas de temps du rayonnement (s)
    22202169           CALL zenang(zlongi,jH_cur,zdtime,rlat,rlon,rmu0,fract)
     2170           PRINT *,'  Lluis: radpas: ',radpas,' fract: ',fract(llp),  &
     2171             ' zdtime: ',zdtime,' jH_cur: ',jH_cur,' rmu0: ',rmu0(llp),&
     2172             ' dtime: ',dtime,' nbapp_rad: ',nbapp_rad
     2173           PRINT *,' Lluis @rlat@: ',rlat
     2174           PRINT *,' Lluis @rlon@: ',rlon
     2175           PRINT *,' Lluis @rmu0@: ',rmu0
     2176           PRINT *,' Lluis @fract@: ',fract
    22212177         ELSE
    22222178           CALL angle(zlongi, rlat, fract, rmu0)
    22232179         ENDIF
    22242180      ENDIF
     2181      PRINT *,'  Lluis zenang_an rlat: ',rlat(llp), ' rlon:',rlon(llp), &
     2182        ' rmu0: ',rmu0(llp),' fract: ',fract(llp),' jH_cur: ',jH_cur,  &
     2183        ' cycle_diurne: ',cycle_diurne
    22252184
    22262185      if (mydebug) then
    … …  
    22582217
    22592218      PRINT *,'  Lluis before pbl_surface qsol: ',qsol(llp),    &
    2260         ' ftsol: ',ftsol(llp,:)
     2219        ' rmu0: ',rmu0(llp),' jH_cur: ',jH_cur
    22612220
    22622221      CALL pbl_surface(                                                              &
    … …  
    22842243       &     zxfluxt,   zxfluxq,   q2m,     fluxq, pbl_tke )
    22852244
    2286       PRINT *,'  Lluis after pbl_surface qsol: ',qsol(llp),      &
    2287         ' ftsol: ',ftsol(llp,:)
     2245      PRINT *,'  Lluis after pbl_surface qsol: ',qsol(llp),    &
     2246        ' rmu0: ',rmu0(llp),' jH_cur: ',jH_cur
    22882247
    22892248!-----------------------------------------------------------------------------------------
    … …  
    22922251       &     (d_u_vdf,d_v_vdf,d_t_vdf+d_t_diss,d_q_vdf,dql0,'vdf')
    22932252!-----------------------------------------------------------------------------------------
    2294       PRINT *,'  Lluis after phys_add_tend qsol: ',qsol(llp),      &
    2295         ' ftsol: ',ftsol(llp,:)
    2296 
    22972253      if (mydebug) then
    22982254        call writefield_phy('u_seri',u_seri,llm)
    … …  
    23762332
    23772333!c Calcule de vitesse verticale a partir de flux de masse verticale
    2378       PRINT *,'    Lluis llp= ',llp,' k flxmass_w airephy omega _______'
    23792334      DO k = 1, klev
    23802335         DO i = 1, klon
    23812336            omega(i,k) = RG*flxmass_w(i,k) / airephy(i)
    2382             IF (i == llp) PRINT *,k,flxmass_w(i,k), airephy(i), omega(i,k)
    23832337         END DO
    23842338      END DO
    … …  
    35603514      ENDIF
    35613515!c
    3562       PRINT *,'  Lluis before calling radlwsw'
     3516      PRINT *,'  Lluis before calling radlwsw rmu0: ',rmu0(llp),   &
     3517        ' dist: ',dist,' frac: ',fract(llp),' jH_cur: ',jH_cur
    35633518         CALL radlwsw                                                                &
    35643519       &        (dist, rmu0, fract,                                                  &
    … …  
    35833538       &        solsw_aero, solsw0_aero,                                             &
    35843539       &        topswcf_aero, solswcf_aero)
     3540      PRINT *,'  Lluis after calling radlwsw rmu0: ',rmu0(llp),   &
     3541        ' dist: ',dist,' frac: ',fract(llp),' jH_cur: ',jH_cur
    35853542         
    35863543!c

  • WRFV3/phys/module_lmdz_phys.F

    r1 r6  
    9191         & wi_start,wi_end,wj_start,wj_end
    9292    INTEGER, INTENT(IN)                                  :: wparfirstscal, wnum3dm
    93     INTEGER, INTENT(IN)                                  :: wjulday
     93    REAL, INTENT(IN)                                     :: wjulday
    9494    REAL, INTENT(IN)                                     :: wtime_step
    9595    INTEGER, INTENT(IN)                                  :: wrf_qvid, wrf_qcid,      &
    … …  
    235235! wnum_tiles: number of tiles
    236236! wjulday: WRF julian day
    237 ! wgmt: WRF gmt hour
     237! wgmt: WRF gmt hour (0, 1)
    238238! wtime_step: WRF time-step (in seconds)
    239239! wrf_dbg: level of debug from WRF
    … …  
    375375      ' wrf_xtime: ',wrf_xtime,' wrf_rst: ', wrf_isrestart,' wrf_lwbdy: ',           &
    376376      wrf_islowbdyin
    377     PRINT *,'     '//TRIM(Spt)//': k p ph ___________'
    378     DO iz=1,dimz
    379       PRINT  *,iz,wrf_perP(ip,iz,jp),wrf_perGeopot(ip,iz,jp)
    380     END DO
    381     PRINT *,'dimx: ',dimx,' dimy: ',dimy,' dimz: ',dimz,' ddimx: ',ddimx,' ddimy: ',ddimy
    382377
    383378    phykeys=0.
    … …  
    390385    lmdz_lafin = .FALSE.
    391386
    392 !    PRINT *,'Lluis: i j ij lon lat lonvals latvals_______'
    393387    DO iy=1,ddimy
    394388      DO ix=1,ddimx
    … …  
    400394        LatValues(ddimx*(iy-1)+ix) = wrf_Lat(ix,iy)*Pi/180.
    401395        PsfctendValues(ddimx*(iy-1)+ix) = wrf_psfctend(ix,iy)
    402 !        PRINT *,ix,iy,(ddimx)*(iy-1)+ix,wrf_Lon(ix,iy),LonValues(ddimx*(iy-1)+ix),   &
    403 !          wrf_Lat(ix,iy),LatValues(ddimx*(iy-1)+ix)
    404396      END DO
    405397    END DO
    … …  
    533525          wrf_grid%loce,wrf_grid%lsic)
    534526
    535         PRINT *,'    Lluis ltksoil. '// Spt //' ter: ',wrf_grid%lter(ip,jp),' lic: ',&
    536           wrf_grid%llic(ip,jp),' loce: ', wrf_grid%loce(ip,jp),' sic: ',             &
    537           wrf_grid%lsic(ip,jp)
    538 
    539527! LANDMASK in LMDZ as fractions
    540528        wrf_grid%lmsk = wrf_grid%lter + wrf_grid%llic
    … …  
    565553! Value is given at the first time-step, not any more, since it will be continuously
    566554!   used by LMDZ
    567          PRINT *,'  Lluis before '//TRIM(Spt)//' SMOIS: ',wrf_grid%smois(ip,:,jp),   &
    568            ' dzs: ',wrf_grid%dzs,' lwsol: ',wrf_grid%lwsol(ip,jp)
    569555! l. Fita, LMD. January 2014
    570556! This should be in this way, but, because we can not use any scheme DZS=0., so
    … …  
    574560!            wrf_grid%soil_layers(:,iz,:)/100.
    575561!        END DO
    576         PRINT *,'  Lluis after: '//TRIM(Spt)//' SMOIS: ',wrf_grid%smois(ip,:,jp),    &
    577            ' lwsol: ',wrf_grid%lwsol(ip,jp)     
    578562        DO iz=1,4
    579563          wrf_grid%ltksoil(:,iz,:)=wrf_grid%st(:,1,:)
    … …  
    626610        ORCHIDEElevels = (/ 0.098, 0.1955, 0.5865, 1.3685, 2.9326, 6.0606, 12.3167,  &
    627611          24.8289, 49.8534, 99.9022, 200. /)
    628         PRINT *,'  Lluis interpolating:',SIZE(wrf_grid%soil_layers(ip,:,jp)),' : ',  &
    629           ' soil: ',                                                                 &
    630           wrf_grid%soil_layers(ip,SIZE(wrf_grid%soil_layers(ip,:,jp)):1:-1,jp),      &
    631           ' st: ',wrf_grid%st(ip,:,jp), 11, ' | ', ORCHIDEElevels,' soiltop: ',      &
    632           wrf_grid%soilctop(ip,:,jp),' soilbot: ',wrf_grid%soilcbot(ip,:,jp)
    633612
    634613        datasetname = soilcat
    … …  
    671650!L. Fita, LMD. January 2014 Soil layers are in cm!
    672651
    673             IF ( (ix == ip) .AND. (iy == jp)) PRINT *,' Lluis '//TRIM(Spt)//         &
    674               ' soiltop: ', wrf_grid%soilctop(ix,:,iy),' soilbot: ',                 &
    675               wrf_grid%soilcbot(ix,:,iy),' layers: ',                                &
    676               wrf_grid%soil_layers(ix,Nsoillayers:1:-1,iy)/100.,' sm: ',             &
    677               wrf_grid%sm(ix,:,iy)
    678652
    679653! NOTE: wrf_dx, wrf_dy should be matrix !!!
    … …  
    684658              soilTmass, wrf_grid%lwsol(ix,iy))
    685659
    686             IF ( (ix == ip) .AND. (iy == jp)) PRINT *,'TOTAL soil mass at '//        &
    687               TRIM(Spt)//': ', soilTmass,' [kg] water content: ',                    &
    688               wrf_grid%lwsol(ix,iy), ' [kg] equiv. water depth: ',                   &
    689               wrf_grid%lwsol(ix,iy)/(wrf_dx*wrf_dy*1000.)
    690 
    691660            wrf_grid%lwsol(ix,iy) =  wrf_grid%lwsol(ix,iy)/(wrf_dx*wrf_dy*1000.)
    692661
    … …  
    694663              wrf_grid%soil_layers(ix,SIZE(wrf_grid%soil_layers(ix,:,iy)):1:-1,iy),  &
    695664              wrf_grid%st(ix,:,iy),11, ORCHIDEElevels,wrf_grid%lotter(ix,:,iy))
    696             IF ( (ix == ip) .AND. (iy == jp)) PRINT *,' Lluis orchidee temperature'//&
    697               ' at '// TRIM(Spt)//': ',wrf_grid%lotter(ix,:,iy)
    698 !            IF ( (ix == ip) .AND. (iy == jp)) STOP
    699665          END DO
    700666        END DO
    … …  
    852818        wrf_grid%lzrugsrel = wrf_grid%lrugKsoil(:,1,:)
    853819
    854         PRINT *,'Lluis Topography values at',ip, jp, ' zmea: ',                      &
    855           wrf_grid%lzmea(ip,jp),' zstd: ',wrf_grid%lzstd(ip,jp),' zsig: ',           &
    856           wrf_grid%lzsig(ip,jp),' zgam: ',wrf_grid%lzgam(ip,jp),' zthe: ',           &
    857           wrf_grid%lzthe(ip,jp),' zpic: ',wrf_grid%lzpic(ip,jp),' zval: ',           &
    858           wrf_grid%lzval(ip,jp),' rugoro: ',wrf_grid%lzrugsrel(ip,jp)
    859 
    860820! R = 8.3144621 ([J mol/K] http://en.wikipedia.org/wiki/Gas_constant)
    861821! mu = 28.97 (Molecular Mass of Air [kg/kmol] from http://www.engineeringtoolbox.com/molecular-mass-air-d_679.html)
    … …  
    867827
    868828!  Grids do not start/end at the same index
    869 !!      PRINT *,'  Lluis: Creating area/cu/cv...'
    870829      wrf_area = 0.
    871830      wrf_cu = 0.
    … …  
    899858      CALL lmdz_vars_init(ddimxy)
    900859
    901 !      PRINT *,'  Lluis before loading 5 11 ltksoil: ',wrf_grid%ltksoil(5,:,11),      &
    902 !        ' ftsol: ',ftsol(315,:)
    903 !      PRINT *,'  Lluis before loading 5 11 qsurf: ',qsurf(315,:),' qsurf_rst: ',     &
    904 !        qsurf_rst(315,:),' smois :',wrf_grid%smois(5,:,11),' lqksoil: ',             &
    905 !        wrf_grid%lqksoil(5,:,11)
    906 
    907860      CALL load_lmdz(wrf_grid,ddimx,ddimy,dimz,wbdyw,ddimxy,4,11,rlon,rlat,zmasq,    &
    908861        pctsrf,ftsol,ftsoil_rst,qsurf_rst,qsol_rst,falb1,falb2,evap_rst,snow_rst,    &
    … …  
    916869        rugos_rst, agesno_rst, ftsoil_rst)
    917870      CALL fonte_neige_init(restart_runoff)
    918 
    919       PRINT *,'  Lluis on module restart_runoff: ',UBOUND(restart_runoff)
    920       PRINT *,'    ffonte_global: ',ALLOCATED(ffonte_global)
    921       PRINT *,'    fqfonte_global: ',ALLOCATED(fqfonte_global)
    922       PRINT *,'    fqcalving_global: ',ALLOCATED(fqcalving_global)
    923 
    924       PRINT *,'  Lluis after loading '//TRIM(Spt)//':  qsurf_rst: ', qsurf_rst(lp,:),&
    925         ' smois :',wrf_grid%smois(ip,:,jp),' lqksoil: ', wrf_grid%lqksoil(ip,:,jp)
    926 
    927       PRINT *,'  Lluis after loading '//TRIM(Spt)//' ltksoil: ',                     &
    928         wrf_grid%ltksoil(ip,:,jp),' ftsol: ',ftsol(lp,:)
    929 
    930       PRINT *,'  Lluis after loading '//TRIM(Spt)//' lwsol: ',wrf_grid%lwsol(ip,jp), &
    931         ' qsol_rst: ',qsol_rst(lp)
    932 
    933       PRINT *,'  Lluis lmdz pt= ',lp,' k qsurf_rst qsol_rst evap_rst snow_rst ' //   &
    934         'fder_rst agesno_rst rugos_rst'
    935       DO iz=1,4
    936         PRINT *,'    * ',iz,qsurf_rst(lp,iz),qsol_rst(lp), evap_rst(lp,iz),          &
    937           snow_rst(lp,iz),fder_rst(lp), agesno_rst(lp,iz),rugos_rst(lp,iz)
    938         PRINT *,'      ftsoil: ',ftsoil_rst(lp,:,iz)
    939       END DO
    940871
    941872      DO ix=1,ddimx
    … …  
    959890      CALL NOread_limit_init(ddimxy, klon_glo)
    960891
    961       PRINT *,'   Lluis klon_glo,', klon_glo,' klon: ',klon,'lp= ',lp,' sst: ',sst(lp)
    962892      sst_glo = sst(1:ddimxy)
    963893! In ocean_forced_mod is used as tsurf_limit = sst(knindex(1:knon))
    964894      tsurf_limit = ftsol(:,3)
    965       PRINT *,'   Lluis pt tsurf_limit_check_____'
    966       DO ix=-10,10
    967         PRINT *,lp+ix,tsurf_limit(lp+ix),sst(lp+ix)
    968       END DO
    969895
    970896! In surf_land_bucket_mod is used as z0_new=rugos(knindex(1:knon),1) [z0_new == z0_limit]
    971897      z0_limit = rugos_rst(:,1)
    972898      alb_limit = albedo
    973       PRINT *,'  Lluis pre-limit_values ftsol: ',ftsol(lp,:),' pctsrf: ',pctsrf(lp,:)
    974       PRINT *,'  Lluis limit_values lp= ',lp,' tsurf_limit: ',tsurf_limit(lp),       &
    975         ' z0_limit: ', z0_limit(lp),' alb_limit: ',alb_limit(lp)
    976899
    977900      rug_glo = rugos_rst(:,is_ter)
    … …  
    981904!      temps=wjulday
    982905
    983 !!      PRINT *,'  Lluis: calling iniphysiq'
    984906      CALL iniphysiq(ddimx*ddimy, dimz, 86400., wjulday, wtime_step,                    &
    985907        & LatValues, LonValues, wrf_area, CUValues, CVValues, 1./reradius, g, r_d, cp,  &
    … …  
    999921      PRINT *,'  Lluis: END of firststep !!!!!'
    1000922    END IF firststep
    1001 
    1002     ip = 12
    1003     jp = 7
    1004     PRINT *,'  ',wrf_grid%itimestep,' Lluis at 12 7: ',wrf_grid%lter(ip,jp), wrf_grid%llic(ip,jp), &
    1005       wrf_grid%loce(ip,jp), wrf_grid%lsic(ip,jp), ' rugksoil: ',                     &
    1006       wrf_grid%lrugksoil(ip,1,jp), wrf_grid%lrugksoil(ip,2,jp),                      &
    1007       wrf_grid%lrugksoil(ip,3,jp), wrf_grid%lrugksoil(ip,4,jp),' rug: ',             &
    1008       wrf_grid%lrug(ip,jp)
    1009     ip = 6
    1010     jp = 9
    1011     PRINT *,'  ',wrf_grid%itimestep,' Lluis at 6 9: ',wrf_grid%lter(ip,jp), wrf_grid%llic(ip,jp), &
    1012       wrf_grid%loce(ip,jp), wrf_grid%lsic(ip,jp), ' rugksoil: ',                     &
    1013       wrf_grid%lrugksoil(ip,1,jp), wrf_grid%lrugksoil(ip,2,jp),                      &
    1014       wrf_grid%lrugksoil(ip,3,jp), wrf_grid%lrugksoil(ip,4,jp),' rug: ',             &
    1015       wrf_grid%lrug(ip,jp)
    1016923
    1017924    ip=icheck_p
    … …  
    1067974      PRINT *,'  ' // subname // ': Problem allocating WRF moisture array!', ierr
    1068975    END IF
    1069 
    1070 !!    PRINT *,'  Lluis wrf MOIST:',LBOUND(wrf_MOIST),' x ',UBOUND(wrf_MOIST)
    1071 !!    PRINT *,'Mix ratio: ',LBOUND(MixingRatioValues),' x ', UBOUND(MixingRatioValues)
    1072976
    1073977    IF (ALLOCATED(MixingRatiotendValues)) DEALLOCATE(MixingRatiotendValues)
    … …  
    12301134    END IF
    12311135
    1232     PRINT *,'  Lluis before physiq '//TRIM(Spt)//' ltksoil: ',                       &
    1233       wrf_grid%ltksoil(ip,:,jp), ' ftsol: ',ftsol(lp,:)
    1234 
    1235     PRINT *,'  Lluis beofre physiq '//TRIM(Spt)//' qsurf: ',zxqsurf(lp),             &
    1236         ' qsurf_rst: ',qsurf_rst(lp,:),' smois :',wrf_grid%smois(ip,:,jp),           &
    1237         ' lqksoil: ', wrf_grid%lqksoil(ip,:,jp)
    1238 
    1239     PRINT *,'  Lluis before physiq: '//TRIM(Spt)//' FullPressValues HalfPressValues '&
    1240        // 'GeopotValues UValues VValues TValues: '
    1241     DO iz=1,dimz
    1242       PRINT *,FullPressValues(lp,iz), HalfPressValues(lp,iz), GeopotValues(lp,iz),   &
    1243         UValues(lp,iz), VValues(lp,iz), TValues(lp,iz)
    1244     END DO
    1245 
    1246     ip = 23
    1247     jp = 18
    1248     PRINT *,'  ',wrf_grid%itimestep,' Lluis before physiq at ',ip,', ', jp,' sst: ', &
    1249       wrf_grid%sst(ip,jp), sst((jp-1)*ddimx+ip),' ftsol 3: ',ftsol((jp-1)*ddimx+ip,is_oce)
    1250 
    12511136    ip=icheck_p
    12521137    jp=jcheck_p
    … …  
    12811166       &  lmdz_dudyn, lmdz_PVtheta)
    12821167
    1283     PRINT *,'  Lluis after physiq '//TRIM(Spt)//' ltksoil: ',                        &
    1284       wrf_grid%ltksoil(ip,:,jp), ' ftsol: ',ftsol(315,:)
    1285     PRINT *,'  Lluis after physiq '//TRIM(Spt)//' qsurf: ',zxqsurf(lp),              &
    1286         ' qsurf_rst: ', qsurf_rst(lp,:),' smois :',wrf_grid%smois(ip,:,jp),          &
    1287         ' lqksoil: ', wrf_grid%lqksoil(ip,:,jp)
    1288     PRINT *,'  Lluis after physiq: '//TRIM(Spt)//' FullPressValues HalfPressValues ' &
    1289        // 'GeopotValuesUValues VValues TValues: '
    1290     DO iz=1,dimz
    1291       PRINT *,FullPressValues(lp,iz), HalfPressValues(lp,iz), GeopotValues(lp,iz),   &
    1292         UValues(lp,iz), VValues(lp,iz), TValues(lp,iz)
    1293     END DO
    1294 
    1295     ip = 23
    1296     jp = 18
    1297     PRINT *,'  ',wrf_grid%itimestep,' Lluis after physiq at ',ip,', ', jp,' sst: ', &
    1298       wrf_grid%sst(ip,jp), sst((jp-1)*ddimx+ip),' ftsol 3: ',ftsol((jp-1)*ddimx+ip,is_oce)
    1299 
    13001168    ip=icheck_p
    13011169    jp=jcheck_p
    … …  
    13701238    wrf_Vtend(wime,:,wjme) = wrfVtenddestagg(wime-1,:,wjme-1)
    13711239
    1372     PRINT *,'  Lluis ' // TRIM(Spt) // ' wrf_Utend: ',wrf_Utend(ip,1,jp),            &
    1373       ' wrfUtenddestagg a: ', wrfUtenddestagg(ip-1,1,jp),' wrfUtenddestagg b: ',     &
    1374       wrfUtenddestagg(ip,1,jp), ' wrfUtenddestagg c: ',wrfUtenddestagg(ip+1,1,jp),   &
    1375       ' UtendValues: ', UtendValues(lp,1),' lp: ',lp
    1376     PRINT *,'dim UtendValues: ',1-wbdyw,':',dimx-1-wbdyw+1,', ',dimz+1,', ',1-wbdyw,':',dimy-1-wbdyw+1
    13771240!!    wrf_Ttend = wrf_Temptend
    13781241
    … …  
    14671330    LMDZvarmethod = 'prod'
    14681331
    1469     PRINT *,'  AFTER physic and before recomputing!!!!!'
    1470     ip = 12
    1471     jp = 7
    1472     PRINT *,'  ',wrf_grid%itimestep,' Lluis at 12 7: ',wrf_grid%lter(ip,jp), wrf_grid%llic(ip,jp), &
    1473       wrf_grid%loce(ip,jp), wrf_grid%lsic(ip,jp), ' rugksoil: ',                   &
    1474       wrf_grid%lrugksoil(ip,1,jp), wrf_grid%lrugksoil(ip,2,jp),                    &
    1475       wrf_grid%lrugksoil(ip,3,jp), wrf_grid%lrugksoil(ip,4,jp),' rug: ',           &
    1476       wrf_grid%lrug(ip,jp)
    1477     ip = 6
    1478     jp = 9
    1479     PRINT *,'  ',wrf_grid%itimestep,' Lluis at 6 9: ',wrf_grid%lter(ip,jp), wrf_grid%llic(ip,jp), &
    1480       wrf_grid%loce(ip,jp), wrf_grid%lsic(ip,jp), ' rugksoil: ',                   &
    1481       wrf_grid%lrugksoil(ip,1,jp), wrf_grid%lrugksoil(ip,2,jp),                    &
    1482       wrf_grid%lrugksoil(ip,3,jp), wrf_grid%lrugksoil(ip,4,jp),' rug: ',           &
    1483       wrf_grid%lrug(ip,jp)
    1484 
    14851332    CALL wrf_varKsoil(wims,wime,wjms,wjme,dimx,dimy,wbdyw,LMDZvarmethod,             &
    14861333      wrf_grid%lter,wrf_grid%llic,wrf_grid%loce,wrf_grid%lsic,                       &
    … …  
    14881335      wrf_grid%lrugksoil(:,4,:),wrf_grid%lrug)
    14891336
    1490     PRINT *,'  AFTER physic and recomputing!!!!!'
    1491     ip = 12
    1492     jp = 7
    1493     PRINT *,'  ',wrf_grid%itimestep,' Lluis at 12 7: ',wrf_grid%lter(ip,jp), wrf_grid%llic(ip,jp), &
    1494       wrf_grid%loce(ip,jp), wrf_grid%lsic(ip,jp), ' rugksoil: ',                   &
    1495       wrf_grid%lrugksoil(ip,1,jp), wrf_grid%lrugksoil(ip,2,jp),                    &
    1496       wrf_grid%lrugksoil(ip,3,jp), wrf_grid%lrugksoil(ip,4,jp),' rug: ',           &
    1497       wrf_grid%lrug(ip,jp)
    1498     ip = 6
    1499     jp = 9
    1500     PRINT *,'  ',wrf_grid%itimestep,' Lluis at 6 9: ',wrf_grid%lter(ip,jp), wrf_grid%llic(ip,jp), &
    1501       wrf_grid%loce(ip,jp), wrf_grid%lsic(ip,jp), ' rugksoil: ',                   &
    1502       wrf_grid%lrugksoil(ip,1,jp), wrf_grid%lrugksoil(ip,2,jp),                    &
    1503       wrf_grid%lrugksoil(ip,3,jp), wrf_grid%lrugksoil(ip,4,jp),' rug: ',           &
    1504       wrf_grid%lrug(ip,jp)
    1505 
    1506     ip = 23
    1507     jp = 18
    1508     PRINT *,'  ',wrf_grid%itimestep,' Lluis at ',ip,', ', jp,' sst : ',              &
    1509       wrf_grid%sst(ip,jp), sst((jp-1)*ddimx+ip),' ftsol 3: ',ftsol((jp-1)*ddimx+ip,is_oce)
    1510  
    15111337    CALL wrf_varKsoil(wims,wime,wjms,wjme,dimx,dimy,wbdyw,LMDZvarmethod,             &
    15121338      wrf_grid%lter,wrf_grid%llic,wrf_grid%loce,wrf_grid%lsic,                       &
    … …  
    15531379!       & wrf_grid%seed2)
    15541380
    1555     ip = 23
    1556     jp = 18
    1557     PRINT *,'  ',wrf_grid%itimestep,' Lluis after put at ',ip,', ', jp,' sst: ',     &
    1558       wrf_grid%sst(ip,jp), sst((jp-1)*ddimx+ip),' ftsol 3: ',ftsol((jp-1)*ddimx+ip,is_oce)
    1559 
    15601381    ip=icheck_p
    15611382    jp=jcheck_p
    … …  
    16541475!       & wrf_grid%lweakinv, wrf_grid%lwind10m, wrf_grid%lwind10max, wrf_grid%lzmax_th)
    16551476
    1656     PRINT *, '   Lluis: wrf_grid%la_th', UBOUND(wrf_grid%la_th),' fraca: ',UBOUND(fraca)
    16571477!    DO iz=wkms, wkme-2
    16581478!      CALL get_lmdz_out3D_z_i(wims, wime, wjms, wjme, iz, ddimx, ddimy, dimz, wbdyw, &
    … …  
    18201640!    END DO
    18211641
    1822     PRINT *,'  Lluis Loading done! iz ldtlwr _______________'
    1823     DO iz=1,dimz
    1824       PRINT *,iz,wrf_grid%ldtlwr(ip,iz,jp)
    1825     END DO
    1826 
    18271642!    CALL get_lmdz_out3D_ii(wims, wime, wjms, wjme, wkms, wkme, ddimx, ddimy, dimz,   &
    18281643!      & wbdyw, ddimxy, nbsrf, nsoilmx, lmdzmixingratios, 1, 2, wrf_l_thermals,       &
Note: See TracChangeset for help on using the changeset viewer.