Changeset 10057


Ignore:
Timestamp:
2018-08-20T17:56:23+02:00 (2 years ago)
Author:
dford
Message:

Fixes for physics-only compilation.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_GO6_package_asm_3d_bgc_v3/NEMOGCM/NEMO/OPA_SRC/ASM/asmbgc.F90

    r10055 r10057  
    990990         ! Select mixed layer 
    991991         IF ( ll_asmdin ) THEN 
     992#if defined key_top 
    992993            CALL ctl_warn( ' Doing direct initialisation with ocean colour assimilation', & 
    993994               &           ' Mixed layer depth taken to be background maximum mld_max_bkg' ) 
    994995            zmld(:,:) = mld_max_bkg(:,:) 
     996#else 
     997            CALL ctl_stop( ' Should not be doing biogeochemical assimilation without key_top' ) 
     998#endif 
    995999         ELSE 
    9961000            SELECT CASE( mld_choice_bgc ) 
     
    18541858 
    18551859         IF ( ln_pno3inc ) THEN 
     1860#if defined key_hadocc || ( defined key_medusa && defined key_foam_medusa ) 
    18561861#if defined key_hadocc 
    18571862            it = jp_had_nut 
    1858 #elif defined key_medusa && defined key_foam_medusa 
     1863#elif defined key_medusa 
    18591864            it = jpdin 
    1860 #else 
    1861             CALL ctl_stop ( ' bgc3d_asm_inc: no compatible BGC model defined' ) 
    18621865#endif 
    18631866            IF ( ln_phytobal ) THEN 
     
    18731876               pno3_bkginc(:,:,:) = pno3_bkginc(:,:,:) - ph_balinc(:,:,:,it) 
    18741877            ENDIF 
     1878#else 
     1879            CALL ctl_stop ( ' bgc3d_asm_inc: no compatible BGC model defined' ) 
     1880#endif 
    18751881         ENDIF 
    18761882 
     
    18781884#if defined key_medusa && defined key_foam_medusa 
    18791885            it = jpsil 
    1880 #else 
    1881             CALL ctl_stop ( ' bgc3d_asm_inc: no compatible BGC model defined' ) 
    1882 #endif 
    18831886            IF ( ln_phytobal ) THEN 
    18841887               psi4_bkginc(:,:,:) = psi4_bkginc(:,:,:) - phyto2d_balinc(:,:,:,it) 
     
    18931896               psi4_bkginc(:,:,:) = psi4_bkginc(:,:,:) - ph_balinc(:,:,:,it) 
    18941897            ENDIF 
     1898#else 
     1899            CALL ctl_stop ( ' bgc3d_asm_inc: no compatible BGC model defined' ) 
     1900#endif 
    18951901         ENDIF 
    18961902 
    18971903         IF ( ln_pdicinc ) THEN 
     1904#if defined key_hadocc || ( defined key_medusa && defined key_foam_medusa ) 
    18981905#if defined key_hadocc 
    18991906            it = jp_had_dic 
    1900 #elif defined key_medusa && defined key_foam_medusa 
     1907#elif defined key_medusa 
    19011908            it = jpdic 
    1902 #else 
    1903             CALL ctl_stop ( ' bgc3d_asm_inc: no compatible BGC model defined' ) 
    19041909#endif 
    19051910            IF ( ln_phytobal ) THEN 
     
    19151920               pdic_bkginc(:,:,:) = pdic_bkginc(:,:,:) - ph_balinc(:,:,:,it) 
    19161921            ENDIF 
     1922#else 
     1923            CALL ctl_stop ( ' bgc3d_asm_inc: no compatible BGC model defined' ) 
     1924#endif 
    19171925         ENDIF 
    19181926 
    19191927         IF ( ln_palkinc ) THEN 
     1928#if defined key_hadocc || ( defined key_medusa && defined key_foam_medusa ) 
    19201929#if defined key_hadocc 
    19211930            it = jp_had_alk 
    1922 #elif defined key_medusa && defined key_foam_medusa 
     1931#elif defined key_medusa 
    19231932            it = jpalk 
    1924 #else 
    1925             CALL ctl_stop ( ' bgc3d_asm_inc: no compatible BGC model defined' ) 
    19261933#endif 
    19271934            IF ( ln_phytobal ) THEN 
     
    19371944               palk_bkginc(:,:,:) = palk_bkginc(:,:,:) - ph_balinc(:,:,:,it) 
    19381945            ENDIF 
     1946#else 
     1947            CALL ctl_stop ( ' bgc3d_asm_inc: no compatible BGC model defined' ) 
     1948#endif 
    19391949         ENDIF 
    19401950 
     
    19421952#if defined key_medusa && defined key_foam_medusa 
    19431953            it = jpoxy 
    1944 #else 
    1945             CALL ctl_stop ( ' bgc3d_asm_inc: no compatible BGC model defined' ) 
    1946 #endif 
    19471954            IF ( ln_phytobal ) THEN 
    19481955               po2_bkginc(:,:,:) = po2_bkginc(:,:,:) - phyto2d_balinc(:,:,:,it) 
     
    19571964               po2_bkginc(:,:,:) = po2_bkginc(:,:,:) - ph_balinc(:,:,:,it) 
    19581965            ENDIF 
     1966#else 
     1967            CALL ctl_stop ( ' bgc3d_asm_inc: no compatible BGC model defined' ) 
     1968#endif 
    19591969         ENDIF 
    19601970 
     
    20932103#endif 
    20942104 
    2095             trn(:,:,:,jp_msa0:jp_msa1) = trn(:,:,:,jp_msa0:jp_msa1) + & 
    2096                &                         pco2_balinc(:,:,:,jp_msa0:jp_msa1) 
    2097             trb(:,:,:,jp_msa0:jp_msa1) = trn(:,:,:,jp_msa0:jp_msa1) 
    2098  
    20992105            IF ( ln_pno3inc ) THEN 
    21002106#if defined key_hadocc 
Note: See TracChangeset for help on using the changeset viewer.