Ignore:
Timestamp:
Apr 18, 2015, 5:06:36 PM (10 years ago)
Author:
jyg
Message:

Bug fixes concerning the number "nqo" of water
phases in "q" array: part II (sequel of fixes of
revision 2262).

Modified files:

carbon_cycle_mod.F90
ini_histrac.h
phyredem.F90
phyetat0.F90
traclmdz_mod.F90
write_histrac.h

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

Legend:

Unmodified
Added
Removed

  • LMDZ5/trunk/libf/phylmd/carbon_cycle_mod.F90

    r1907 r2265  
    157157    itc=0
    158158    DO it=1,nbtr
    159        iiq=niadv(it+2)
     159!!       iiq=niadv(it+2)                                                            ! jyg
     160       iiq=niadv(it+nqo)                                                            ! jyg
    160161       
    161162       SELECT CASE(tname(iiq))

  • LMDZ5/trunk/libf/phylmd/ini_histrac.h

    r1907 r2265  
    3030!----------------
    3131     DO it = 1,nbtr
    32         iiq = niadv(it+2)
     32!!        iiq = niadv(it+2)                                                         ! jyg
     33        iiq = niadv(it+nqo)                                                         ! jyg
    3334
    3435! CONCENTRATIONS

  • LMDZ5/trunk/libf/phylmd/phyetat0.F90

    r2252 r2265  
    2121       zstd, zthe, zval, ale_bl, ale_bl_trig, alp_bl
    2222  USE iostart, ONLY : close_startphy, get_field, get_var, open_startphy
    23   USE infotrac, only: nbtr, type_trac, tname, niadv
     23  USE infotrac, only: nbtr, nqo, type_trac, tname, niadv
    2424  USE traclmdz_mod,    ONLY : traclmdz_from_restart
    2525  USE carbon_cycle_mod, ONLY : carbon_cycle_tr, carbon_cycle_cpl, co2_send
    … …  
    349349
    350350  IF (type_trac == 'lmdz') THEN
    351      DO it=1, nbtr
    352         iiq=niadv(it+2)
     351     DO it=1, nbtr                                                                 
     352!!        iiq=niadv(it+2)                                                           ! jyg
     353        iiq=niadv(it+nqo)                                                           ! jyg
    353354        found=phyetat0_get(1,trs(:,it),"trs_"//tname(iiq), &
    354355              "Surf trac"//tname(iiq),0.)

  • LMDZ5/trunk/libf/phylmd/phyredem.F90

    r2255 r2265  
    337337     CALL traclmdz_to_restart(trs)
    338338     DO it=1, nbtr
    339         iiq=niadv(it+2)
     339!!        iiq=niadv(it+2)                                                           ! jyg
     340        iiq=niadv(it+nqo)                                                           ! jyg
    340341        CALL put_field("trs_"//tname(iiq), "", trs(:, it))
    341342     END DO

  • LMDZ5/trunk/libf/phylmd/traclmdz_mod.F90

    r2180 r2265  
    174174    id_pcsat=0; id_pcocsat=0; id_pcq=0; id_pcs0=0; id_pcos0=0; id_pcq0=0
    175175    DO it=1,nbtr
    176        iiq=niadv(it+2)
    177        IF ( tname(iiq) == "RN" ) THEN
     176!!       iiq=niadv(it+2)                                                            ! jyg
     177       iiq=niadv(it+nqo)                                                            ! jyg
     178       IF ( tname(iiq) == "RN" ) THEN                                               
    178179          id_rn=it ! radon
    179180       ELSE IF ( tname(iiq) == "PB") THEN
    … …  
    293294! ----------------------------------------------
    294295    DO it=1,nbtr
    295        iiq=niadv(it+2)
     296!!       iiq=niadv(it+2)                                                            ! jyg
     297       iiq=niadv(it+nqo)                                                            ! jyg
    296298       ! Test if tracer is zero everywhere.
    297299       ! Done by master process MPI and master thread OpenMP

  • LMDZ5/trunk/libf/phylmd/write_histrac.h

    r1907 r2265  
    1717!----------------
    1818     DO it=1,nbtr
    19         iiq=niadv(it+2)
     19!!        iiq=niadv(it+2)                                                           ! jyg
     20        iiq=niadv(it+nqo)                                                           ! jyg
    2021
    2122! CONCENTRATIONS
Note: See TracChangeset for help on using the changeset viewer.