Changeset 2790 for LMDZ5/trunk


Ignore:
Timestamp:
Feb 8, 2017, 4:05:38 PM (8 years ago)
Author:
acozic
Message:

Modification for INCA
If there is no value for O3 concentration in start.nc, we decide to read value of OX (OX=O3 + O1D + O) instead of initialized to zero.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/dyn3dmem/dynetat0_loc.f90

    r2622 r2790  
    153153  DO iq=1,nqtot
    154154    var=tname(iq)
     155#ifdef INCA
     156    IF (var .eq. "O3" ) THEN
     157       IF(NF90_INQ_VARID(fID,var,vID) == NF90_NoErr) THEN
     158          CALL get_var2(var,q_glo); q(ijb_u:ije_u,:,iq)=q_glo(ijb_u:ije_u,:); CYCLE
     159       ELSE
     160          WRITE(lunout,*) 'Tracer O3 is missing - it is initialized to OX'
     161          IF(NF90_INQ_VARID(fID,"OX",vID) == NF90_NoErr) THEN
     162             CALL get_var2("OX",q_glo); q(ijb_u:ije_u,:,iq)=q_glo(ijb_u:ije_u,:); CYCLE
     163          ENDIF
     164       ENDIF
     165    ENDIF
     166#endif
    155167    IF(NF90_INQ_VARID(fID,var,vID)==NF90_NoErr) THEN
    156168      CALL get_var2(var,q_glo); q(ijb_u:ije_u,:,iq)=q_glo(ijb_u:ije_u,:); CYCLE
Note: See TracChangeset for help on using the changeset viewer.